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Abstract 

In  this  thesis,  I  present  high-resolution  stable-isotope  and  planktonic-fauna  records  from 
Bering  Sea  sediment  cores,  spanning  the  time  period  from  50,000  years  ago  to  the  present. 
During  Marine  Isotope  Stage  3  (MIS3)  at  30-20  ky  BP  (kiloyears  before  present)  in  a  core 
from  1467  m  water  depth  near  Umnak  Plateau,  there  were  episodic  occurrences  of  diagenetic 
carbonate  minerals  with  very  low  <513C  (— 22.4  %o),  high  £180  (6.5  %o),  and  high  [Mg]/[Ca], 
which  seem  associated  with  sulfate  reduction  of  organic  matter  and  possibly  anaerobic 
oxidation  of  methane.  The  episodes  lasted  less  than  1000  years  and  were  spaced  about  1000 
years  apart.  During  MIS3  at  55-20  ky  BP  in  a  core  from  2209  m  water  depth  on  Bowers 
Ridge,  N.  pachyderma  (s.)  and  Uvigerina  <5180  and  513C  show  no  coherent  variability  on 
millennial  time  scales. 

Bering  Sea  sediments  are  dysoxic  or  laminated  during  the  deglaciation.  A  high  sedi¬ 
mentation  rate  core  (200cm/ky)  from  1132  m  on  the  Bering  Slope  is  laminated  during  the 
B0lling  warm  phase,  Allerpd  warm  phase,  and  early  Holocene,  where  the  ages  of  lithological 
transitions  agree  with  the  ages  of  those  climate  events  in  Greenland  (GISP2)  to  well  within 
the  uncertainty  of  the  age  models.  The  subsurface  distribution  of  radiocarbon  was  esti¬ 
mated  from  a  compilation  of  published  and  unpublished  North  Pacific  benthic-planktonic 
14C  measurements  (475-2700  m  water  depth).  There  was  no  consistent  change  in  14C  pro¬ 
files  between  the  present  and  the  Last  Glacial  Maximum,  Bplling-Allerpd,  or  the  Younger 
Dryas  cold  phase.  N.  pachyderma  (s.)  5180  in  the  Bering  Slope  core  decreases  rapidly  (in 
less  than  220  y)  by  0. 7-0.8  %o  at  the  onset  of  the  Bplling  and  the  end  of  the  Younger  Dryas. 
These  isotopic  shifts  are  accompanied  by  transient  decreases  in  the  relative  abundance  of 
N.  pachyderma  (s.),  suggesting  that  the  isotopic  events  are  transient  warmings  and  sus¬ 
tained  freshenings. 


3 


Thesis  Advisor: 


Dr.  Lloyd  Keigwin 


Thesis  Committee: 

Dr.  William  Curry 
Dr.  Olivier  Marchal 
Dr.  Daniel  McCorkle 
Dr.  Jerry  McManus  (chair) 
Prof.  Carl  Wunsch 


5 


Biographical  Note 


The  author  was  born  on  25  January,  1977,  in  Englewood,  Colorado.  She  grew  up  in  the 
Washington,  D.  C.,  area  and  graduated  from  the  Thomas  Jefferson  High  School  for  Science 
and  Technology  in  Alexandria,  Virginia,  in  1995.  She  then  attended  Princeton  University, 
and  in  1999,  received  an  A.  B.  in  Geosciences,  magna  cum  laude ,  and  the  Edward  Sampson, 
Class  of  1914,  Prize  in  Economic  Geology.  Her  senior  thesis,  Late  Maastrichtian  Popula¬ 
tions  of  Planktonic  Foraminifera  of  the  South  Atlantic ,  was  advised  by  Prof.  Gerta  Keller. 
She  also  received  a  Certificate  in  Violoncello  and  Viola  da  Gamba  Performance  from  the 
Department  of  Music,  where  she  studied  with  Vivian  Barton  Dozor.  In  Fall,  1999,  she 
enrolled  in  the  MIT/WHOI  Joint  Program  in  Oceanography,  where  she  was  a  recipient  of 
the  WHOI  Stanley  Watson  Fellowship.  She  won  the  IMAGES  Student  Poster  Prize  at  the 
VHth  and  VHIth  International  Conferences  on  Paleoceanography  in  2001  and  2004.  At 
MIT,  she  was  awarded  a  Cherry  L.  Emerson  Music  Fellowship,  and  studied  violoncello  with 
Dr.  Carol  Ou  at  the  New  England  Conservatory  of  Music. 


7 


Acknowledgments 

The  work  in  this  thesis  was  supported  by  the  National  Science  Foundation  award  OPP- 
9912122  to  Lloyd  Keigwin,  the  Oak  Foundation  of  Boston,  Massachusetts,  the  Stanley 
Watson  Fellowship,  the  Paul  Fye  Fellowship,  and  the  Academic  Programs  Office  at  WHOI. 

My  thesis  advisor,  Lloyd  Keigwin,  gave  me  a  wonderful  opportunity  in  working  on  these 
Healy  02-02  cores.  I’d  like  to  thank  Lloyd  and  the  rest  of  my  thesis  committee,  Bill  Curry, 
Olivier  Marchal,  Dan  McCorkle,  Jerry  McManus,  and  Carl  Wunsch  for  their  guidance, 
especially  during  the  writing  of  this  thesis. 

This  project  would  not  have  been  possible  without  the  support  of  many  people  who 
helped  me  with  the  analytical  work.  Sue  Trimarchi  and  Luping  Zou  taught  me  how  to 
process  and  pick  foraminifer  samples.  Ellen  Roosen,  Mary  Carman,  Marti  Jeglinski,  and 
Kathy  Elder  helped  with  sample  preparation  and  picking.  A1  Gagnon,  Rindy  Ostermann, 
and  Marti  Jeglinski  taught  me  to  how  to  run  the  mass  spectrometers.  Kristy  Dahl  and  Rose 
Came  helped  me  with  preparing  and  running  ICP-MS  samples.  Joan  Bernhard  provided 
the  microscope  I  used  to  photograph  foraminifer  specimens.  All  the  14 C  analysis  were 
conducted  at  the  NOSAMS  facility. 

It  was  a  privilege  to  be  a  student  in  the  MIT/ WHOI  Joint  Program.  The  staff  of 
the  Academic  Programs  Office  and  JCMG&G  were  very  helpful  and  supportive.  I  was 
fortunate  to  enter  MIT  and  the  Joint  Program  with  a  cohort  of  paleoceanographers,  Rose 
Came,  Kristy  Dahl,  Peter  Huybers,  Dave  Lund,  and  Matt  Makou,  from  whom  I  learned 
much  and  who  all  became  good  friends.  Pd  like  to  thank  my  classmates  and  friends  in 
Woods  Hole  who  made  this  a  wonderful  place  to  live  and  work,  especially  Hugh,  Margaret, 
Linda,  Heidi,  Jonathan,  John,  Phoebe,  Adam,  Jeff  F.,  Tetjana,  Joe,  Marsha,  Julia,  Jayne, 
Mark,  Liz,  Cara,  Jeff  S.,  Rhea,  Boa,  Patty,  Luc,  Jessica,  Carlos,  and  Nathalie.  Pm  grateful 
to  Luping  Zou  and  Xin  Huang  for  their  friendship  and  for  being  a  second  family  to  me  here 
in  Woods  Hole. 

Finally,  Pd  like  to  thank  my  parents,  Joung  and  Jungnim  Cook,  and  my  sister,  Jean. 


9 


Contents 


1  Introduction  13 

2  The  deglacial  history  of  surface  and  intermediate  water  of  the  Bering  Sea  19 


3  The  deglaciation  in  the  Bering  Sea  31 

3.1  Introduction .  31 

3.2  Methods .  35 

3.3  Results . 43 

3.4  Discussion .  51 

3.5  Conclusions . 60 

4  Marine  Isotope  Stage  3  in  the  Bering  Sea  61 

4.1  Introduction .  61 

4.2  Methods .  63 

4.3  Results .  64 

4.4  Discussion .  64 

4.5  Conclusions .  68 


5  Evidence  of  Methane  in  Bering  Sea  Sediments  during  Marine  Isotope 


Stage  3  69 

5.1  Introduction .  69 

5.2  Methods .  75 

5.3  Results .  76 

5.4  Discussion .  81 

5.5  Conclusions  and  Future  Work .  88 


11 


6  Summary 


91 


A  Apparent  ventilation  age  of  the  North  Pacific  during  the  last  deglacial 
transition  97 

A.l  Benthic-planktonic  14C  age  difference .  97 

A. 2  A14C  of  subsurface  water .  99 

A. 3  Apparent  ventilation  age .  104 

A. 4  Discussion .  107 


12 


Chapter  1 


Introduction 


The  ocean  circulation  is  an  integral  part  of  Earth’s  climate  system.  Seawater  holds  large 
amounts  of  heat  and  dissolved  carbon  dioxide,  a  strong  greenhouse  gas.  As  the  ocean 
circulates,  these  quantities  are  moved  around  the  globe  and  exchanged  with  the  atmosphere. 
In  the  modern  ocean,  the  physical  properties  of  seawater,  such  as  salinity  and  temperature, 
and  the  distribution  of  dissolved  components,  such  as  nutrients,  oxygen,  and  carbon  dioxide, 
are  a  product  of  the  complex  interaction  of  circulation,  mixing,  and  biogeochemical  processes 
in  the  ocean 

Since  the  1960s,  paleoceanographers  have  been  developing  geochemical  proxies  of  the 
properties  of  seawater  in  the  past.  Now  we  have  a  wide  array  of  tools  with  which  to  study 
what  was  the  state  of  the  ocean  as  far  back  in  time  as  the  age  of  the  oldest  sediments  on 
the  seafloor.  When  these  marine  records  are  combined  with  continental  records  of  changes 
in  flora  and  ice  sheet  and  glacier  extent,  and  with  ice  core  records  of  the  composition  of  the 
ancient  atmosphere  and  precipitation,  we  can  get  a  sense  of  how  the  Earth’s  climate  system 
has  changed  in  the  past,  and  even  formulate  tentative  theories  of  why  it  changed. 

The  Pleistocene  (the  last  2  million  years)  is  the  most  recent  episode  of  Earth’s  history 
in  which  global  climate  was  cold  enough  for  there  to  be  ice  year-round  at  the  poles,  and 
where  the  Earth  cycles  into  and  out  of  ice  ages.  The  focus  of  this  thesis  is  the  time  from 
50,000  years  ago  (the  middle  of  the  last  glacial  period)  to  the  present.  During  this  time 
period,  superimposed  upon  the  large-scale  shift  from  a  glacial  state  to  an  interglacial  state 
(which  took  place  18,000-10,000  years  ago)  are  millennial-scale  climate  fluctuations  (events 
of  duration  and  spacing  on  the  order  of  1000  y)  which  are  about  half  the  amplitude  of  the 
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glacial-interglacial  change,  and  appear  to  have  had  the  greatest  impact  in  the  North  Atlantic 
region.  In  this  thesis,  I  study  the  manifestation  of  these  shorter  climate  events  in  the  North 
Pacific  in  order  to  learn  what  role  the  physical  and  biological  processes  in  this  region  play 
in  the  global  climate  system. 

Study  of  how  the  ocean,  atmosphere,  and  cryosphere  change  under  a  variety  of  forcings 
and  boundary  conditions  in  the  past  helps  us  to  understand  Earth’s  climate  system  better 
than  studying  just  the  state  of  the  modern  ocean.  This  is  vitally  important  today,  as  we 
humans  are  rapidly  altering  the  concentration  of  atmospheric  greenhouse  gases,  we  can 
apply  our  understanding  of  the  past  in  predicting  what  effects  human  activity  will  have  on 
climate,  and  how  this  will  impact  human  civilization. 

One  difficulty  in  studying  the  ancient  ocean  is  that  there  are  relatively  few  places  where 
marine  sediments  suitable  for  paleoceanographic  study  can  be  recovered.  So  the  geological 
record  has  rather  poor  spatial  coverage,  limited  to  the  ocean  margins  and  bathymetric  highs, 
such  as  seamounts,  hotspots  chains,  and  mid-ocean  ridges.  Another  difficulty  is  that  each  of 
the  geochemical  proxies  we  use  to  reconstruct  ocean  water  properties  is  imperfect,  so  that 
the  interpretation  of  any  one  proxy  record  is  fraught  with  enough  caveats  to  not  make  a 
very  compelling  story  on  its  own.  Yet  another  difficulty  is  that  the  ocean  is  heterogeneous, 
and  local  processes  can  influence  a  geological  record  as  much  as  the  basin- wide  processes  of 
interest.  So  reconstructing  what  the  ocean  was  like  in  the  past  is  a  thorny  problem! 

The  great  strength  of  our  arsenal  of  paleoceanographic  tools  is  that  there  is  often  more 
than  one  proxy  we  can  use  to  reconstruct  a  particular  property  of  the  ancient  ocean.  So  our 
strategy  is  to  build  a  database  of  sediment  records  from  various  locations,  using  more  than 
one  tool.  This  way,  if  several  records  of  independent  proxies  tell  a  consistent  story,  we  can 
believe  the  story  much  more  than  if  we  just  had  one  record  in  one  place.  The  process  of 
collecting  and  analyzing  sediment  cores  is  time-consuming,  labor-intensive,  and  expensive, 
so  it  can  take  years  of  work  by  multiple  specialists  to  find  out  what  story  one  sediment  core 
has  to  tell. 

Foraminifera,  single-celled  organisms  that  live  near  the  sea  surface  and  on  the  seafloor, 
are  the  workhorse  of  paleoceanography.  They  make  their  shells  out  of  calcium  carbonate, 
which  records  the  chemistry  of  the  environment  in  which  they  lived.  Calcium  carbonate  is 
well  preserved  in  the  Atlantic,  but  poorly  preserved  in  the  Pacific,  where  bottom  waters  are 
more  corrosive  with  respect  to  calcium  carbonate.  As  a  consequence,  there  are  relatively  few 
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Figure  1-1:  The  Bering  Sea.  The  scale  for  color  bathymetry  and  topography  is  in  meters. 
Coring  sites  from  HLY02-02  are  labeled  with  black  dots. 

high-resolution  paleoceanographic  records  from  the  Pacific,  and  it  is  not  well  known  what 
happened  in  this  largest  of  ocean  basins  during  the  Pleistocene.  Arguably  the  best  of  these 
records  come  from  the  Santa  Barbara  Basin.  However,  as  a  coastal  site,  the  Santa  Barbara 
Basin  is  influenced  by  coastal  upwelling,  and  it  is  ambiguous  how  much  the  geological  records 
from  this  site  reflect  local  processes,  and  how  much  they  reflect  what  is  occurring  in  the 
open  North  Pacific.  In  fact,  there  are  too  few  high-quality  North  Pacific  records  for  there 
to  be  much  of  any  consensus  as  to  what  are  the  most  important  physical  and  biological 
processes  in  this  basin  with  respect  to  climate,  and  how  they  changed  in  the  past. 

For  this  thesis,  I  studied  sediment  cores  that  were  collected  from  the  Bering  Sea  by  the 
science  party  lead  by  my  thesis  advisor,  Lloyd  Keigwin,  Neal  Driscoll  of  the  Scripps  Insti¬ 
tute  of  Oceanography,  and  Julie  Brigham-Grette  of  University  of  Massachusettes.  Amherst, 
aboard  the  USCGC  Healy  in  summer,  2002.  The  Bering  Sea  is  the  northern-most  marginal 
basin  of  the  Pacific.  The  sediment  cores  were  collected  from  the  Alaskan  continental  slope  in 
the  northeast  Bering  Sea,  the  Umnak  Plateau  in  the  southeast  Bering  Sea.  and  from  Bow¬ 
ers  Ridge,  an  aseismic  ridge  in  the  south-central  Bering  Sea  (Figure  1-1).  From  foraminifer 
shells  preserved  in  these  cores,  I  have  constructed  records  of  £l80,  a  temperature  and  salin¬ 
ity  proxy.  6[S C,  a  nutrient  proxy,  and  14C,  which  can  be  used  to  translate  depth  in  the 
sediment  core  to  calendar  age.  When  1  [C  is  measured  in  benthic  foraminifera.  which  live 
on  the  seafloor,  it  can  be  used  to  deduce  how  long  the  water  mass  bathing  the  seafloor 
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has  been  out  of  contact  with  the  atmosphere,  providing  insight  on  ocean  circulation.  I  also 
produced  records  of  the  relative  abundance  of  a  planktonic  (living  near  the  surface  ocean) 
foraminifer,  N.  pachy derma  (s.),  which  is  a  qualitative  proxy  for  sea  surface  temperature. 
The  temporal  resolution  of  the  records  presented  in  this  thesis  is  as  good  as  or  better  than 
most  published  North  Pacific  records,  good  enough  to  resolve  millennial-scale  events  if  they 
occured  in  this  basin. 

The  preliminary  stratigraphic  data  of  the  Healy  cores  are  included  in  a  paper  in  Deep 
Sea  Research  II  (Cook  et  ah,  2005),  where  I  report  the  discovery  of  cores  with  high  sediment 
accumulation  rate  and  good  carbonate  preservation.  In  this  paper  (Chapter  2  of  the  thesis), 
I  reviewed  the  hydrography  of  the  Bering  Sea,  why  the  region  is  of  paleoceanographic 
interest,  and  previous  studies  that  were  conducted  in  the  area. 

In  Chapter  3,  I  present  the  lithology  of  Bering  Sea  sediments,  and  high-resolution  data 
sets  of  N.  pachyderma  (s.)  relative  abundance  and  618 O  from  Bering  Slope  and  Bowers 
Ridge  that  span  the  Last  Glacial  Maximum  (LGM,  ^21,000  years  ago)  to  the  present. 
I  also  compile  published  and  unpublished  measurements  of  14C  in  coeval  planktonic  and 
benthic  foraminifera  from  the  North  Pacific  and  estimate  the  14C  of  subsurface  water  in  the 
North  Pacific  during  the  LGM  and  deglaciation  (14C  calculations  are  detailed  in  Appendix 
A).  The  goals  of  this  chapter  are  (1)  to  understand  the  history  of  temperature  and  salinity 
of  Bering  Sea  surface  water,  (2)  to  look  for  evidence  of  whether  the  circulation  in  the  North 
Pacific  changed  between  the  LGM  and  the  present,  and  during  the  B0lling-Aller0d  and 
Younger  Dryas,  two  climate  events  during  the  deglaciation.  A  change  in  circulation  of  the 
North  Pacific  could  affect  both  the  oxygen  content  (oxygen  affects  lithology),  and  the  14C 
content  of  subsurface  water. 

In  Chapter  4,  using  planktonic  and  benthic  <5lsO  and  <513C  as  well  as  relative  abundance 
of  N.  pachyderma  (s.)  in  a  core  from  Bowers  Ridge,  I  test  the  hypothesis  of  previous  investi¬ 
gators  that  the  sea  surface  temperature  and  the  <513C  of  intermediate-depth  water  changed 
on  millennial  time  scales  during  the  last  glacial  period  (55,000-28,000  years  ago).  In  Chap¬ 
ter  5,  I  present  planktonic  and  benthic  Sls O  and  5l3C  from  a  core  from  the  Umnak  Plateau 
region  from  the  same  time  period.  In  this  core,  there  are  large  excursions  in  benthic  and 
planktonic  S13 C  which  I  hypothesize  are  caused  by  the  precipitation  of  carbonate  minerals 
on  the  foraminifer  shells  after  they  were  buried  in  the  sediment  column.  I  formulate  a  sim¬ 
ple  mathematical  model  to  calculate  the  isotopic  composition  of  the  carbonate  overgrowths, 
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in  order  to  understand  why  these  minerals  formed.  In  Chapter  6,  I  summarize  the  main 
findings  of  this  thesis  and  discuss  future  research  directions. 
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Chapter  2 


The  deglacial  history  of  surface 
and  intermediate  water  of  the 
Bering  Sea 


This  chapter  was  published  as  Cook,  M.  S.,  Keigwin,  L.  D.,  Sancetta,  C.  S.,  2005.  The 
deglacial  history  of  surface  and  intermediate  water  of  the  Bering  Sea.  Deep-Sea  Research 
II  52,  2163-2173. 
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Abstract 

The  lithology  of  deglacial  sediments  from  the  Bering  Sea  includes  intervals  of  laminated  or  dysaerobic  sediments.  These 
intervals  are  contemporaneous  with  the  occurrence  of  laminated  sediments  from  the  California  margin  and  Gulf  of 
California,  which  suggests  widespread  low-oxygen  conditions  at  intermediate  depths  in  the  North  Pacific  Ocean.  The  cause 
could  be  reduced  intermediate  water  ventilation,  increased  organic  carbon  flux,  or  a  combination  of  the  two.  We  infer 
abrupt  decreases  of  planktonic  foraminifer  <5lsO  at  14,400  and  11,650  yr  BP,  which  may  be  a  combination  of  both 
freshening  and  wanning.  On  the  Shirshov  Ridge,  the  abundance  of  sea-ice  diatoms  of  the  genus  Nitzschia  reach  local 
maxima  twice  during  the  deglaciation,  the  latter  of  which  may  be  an  expression  of  the  Younger  Dryas.  These  findings 
expand  the  extent  of  the  expression  of  deglacial  millennial-scale  climate  events  to  include  the  northernmost  Pacific. 

©  2005  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

Water-mass  tracer  proxies  measured  from  the 
most  intense  cold  events  of  the  last  deglaciation, 
Heinrich  Event  1  and  the  Younger  Dryas  (YD),  are 
consistent  with  a  reduction  in  the  strength  of  the 
Atlantic  meridional  overturning  circulation  (MOC) 
(Boyle  and  Keigwin,  1987;  Zahn  et  al.,  1997;  Vidal 
et  al.,  1997;  McManus  et  al.,  2004).  Because  of  the 
large  transport  of  heat  associated  with  the  MOC, 
60%  of  the  Northern  Hemisphere  midlatitude 
poleward  ocean  heat  transport  takes  place  in  the 
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Atlantic  (Talley,  2003),  which  is  half  the  size  of  the 
Pacific.  There  is  no  analogous  circulation  like  the 
MOC  in  the  Pacific  today,  where  heat  transport  is 
dominated  by  that  associated  with  the  shallow 
wind-driven  gyre  circulation.  The  lack  of  deep  or 
intermediate  convection  is  attributed  to  the  low 
salinity  of  North  Pacific  surface  waters,  which  even 
at  their  freezing  point,  are  not  dense  enough  to  sink 
(Warren,  1983;  Emile-Geay  et  al.,  2003). 

In  the  modern  North  Pacific,  freshening  and  an 
increase  in  apparent  oxygen  utilization  of  NPIW  in 
the  last  60  yr  are  linked  to  hydrographic  changes  in 
the  Western  Subarctic  Gyre  (WSAG)  and  the 
Okhotsk  Sea  (Andreev  and  Watanabe,  2002;  Joyce 
and  Dunworth-Baker,  2003;  Hill  et  al.,  2003), 
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demonstrating  a  link  between  surface  and  inter¬ 
mediate  water  physical  properties.  From  the  pre¬ 
sence  of  CFCs  in  the  deepest  Bering  Sea  Basin 
(Warner  and  Roden,  1995),  there  is  a  very  small 
amount  of  deep-water  formation  occurring  there 
today.  Deep  water  in  the  rest  of  the  North  Pacific 
originates  from  the  Southern  Ocean  (Roden,  1995). 

From  geochemical  proxies  of  nutrients,  dissolved 
oxygen,  temperature,  and  salinity  within  intermedi¬ 
ate  (Keigwin  and  Jones,  1990;  Behl  and  Kennett, 
1996;  Zheng  et  ah,  2000;  Keigwin,  1998;  Ahagon 
et  al.,  2003)  and  deep  water  (Herguera  et  al.,  1992; 
Lund  and  Mix,  1998;  Matsumoto  et  al.,  2002),  it  is 
hypothesized  that  the  ventilation  of  the  intermediate 
and  the  deep  Pacific  water  masses  was  higher  in  the 
past.  This  implies  that  the  balance  of  freshwater  and 
salt  in  the  ocean  basins  would  have  been  signifi¬ 
cantly  different  than  today,  possibly  with  greater 
oceanic  heat  transport  in  the  Pacific.  Unfortunately, 
detailed  interpretation  of  the  paleo-geochemical 
evidence  at  most  open-ocean  sites  is  not  possible 
because  of  low  sediment  accumulation  rates. 

The  high-resolution  records  in  laminated  sedi¬ 
ments  from  anoxic  basins  from  the  California  coast 
and  the  Gulf  of  California  have  been  extensively 
studied  (Bchl  and  Kennett,  1996;  Zheng  et  al.,  2000; 
van  Geen  et  al.,  2003;  Ortiz  et  al.,  2004).  But 
because  productivity  in  California  margin  anoxic 


basins  is  a  function  of  coastal  upwelling,  it  is 
difficult  to  distinguish  between  reduced  ventilation 
and  local  increase  in  export  production  as  the  cause 
of  laminations. 

The  Bering  Sea’s  proximity  to  the  Okhotsk  Sea, 
which  is  the  source  of  the  precursor  to  modern 
NPIW,  and  its  position  as  the  gateway  between  the 
North  Pacific  and  the  Arctic  make  it  a  desirable 
place  to  study  Pacific  paleoclimatology.  In  addition, 
coring  of  the  thickly  sedimented  continental  margin 
(Carlson  and  Karl,  1988)  and  aseismic  ridges  allows 
for  high-resolution  studies. 

On  the  USCGC  Healy  in  summer,  2002,  we 
collected  new  sediment  cores  from  three  regions: 
Bowers  Ridge,  an  aseismic  ridge  in  south-central 
Bering  Sea;  Umnak  Plateau,  on  the  southeast  end  of 
the  continental  slope;  and  between  the  Pervenets 
and  Navarin  Canyons,  on  the  northwest  end  of  the 
slope  (Fig.  1,  Table  1).  Our  preliminary  data 
indicate  many  of  these  cores  have  very  high 
sediment  accumulation  rates  and  good  calcium 
carbonate  (CaCOj)  preservation  and  will  provide 
an  excellent  record  of  rapid  climate  variability  over 
the  last  complete  glacial  cycle.  Most  coring  loca¬ 
tions  are  far  enough  away  from  continental  land 
masses  to  be  less  influenced  by  coastal  upwelling. 

Due  to  the  low  CaC03  accumulation  rates 
endemic  to  the  Pacific,  most  of  the  previous  marine 


180° 


Fig.  I .  Bering  Sea  core  locations.  Bathymetric  contours  are  at  1,2,  and  3  km.  Details  of  the  numbered  core  locations  are  in  Table  I .  Mean 
surface  circulation  is  indicated  (modified  from  Stabeno  et  al.  (1999)).  Currents  (rectangles)  and  geographic  features  (circles)  referred  to  in 
the  text  are  labeled:  Kamchatka  Current  (KC),  Bering  Slope  Current  (BSC).  Alaskan  Stream  (AS),  Bering  Strait  (bs).  Kamchatka  Strait 
(ks).  Near  Strait  (ns),  Buldir  Pass  (bp),  Amchitka  Strait  (as),  Amutka  Pass  (ap),  and  Unimak  Pass  (up). 
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Table  I 

The  cores  used  in  this  study 


Location 

U 

Core 

Depth  (m) 

Latilude/Longitude 

Bering  Slope 

1 

HLY02-02-3JPC 

1132 

N60.12797EI  80.5582° 

Umnak  Plateau 

2 

H  LY02-02-51JPC 

1467 

N54.55327E  19 1.3331° 

Umnak  Plateau 

3 

RC 14-121 

2532 

N54.85007E1 89.3233° 

Bowers  Ridge 

4 

H  LY02-02-17JPC 

2209 

N53.9330°/EI  78.6988° 

Bowers  Ridge 

5 

VINOI9-I2GGC 

2745 

N53.73677EI  78.7083° 

Shirshov  Ridge 

6 

A  NT  P-4  2  PC 

2393 

N56.38707E17I.07700 

paleoclimatological  work  in  the  Bering  Sea  region  is 
based  on  siliceous  microfossil  abundances  (Zhuze. 
1962;  Sanectta  and  Robinson.  1983;  Sancetta  et  al.. 
1985;  Morley  and  Robinson.  1986;  Starratt,  1993). 
In  particular,  VoNitzschia  indicates  greater  sea-ice 
coverage  in  the  last  glacial  period  along  the  eastern 
shore  of  the  Bering  Sea,  including  the  Pervenets 
Canyon  region  and  the  Umnak  Plateau.  Modern 
distribution  of  Nitzschia  is  consistent  with  observed 
sea-ice  extent,  where  the  southern  slope  including 
the  Umnak  Plateau  is  relatively  ice-free. 

From  the  abundance  of  low-oxygen  benthic 
foraminifer  species,  near  Pervenets  Canyon,  Starratt 
(1995)  inferred  that  the  oxygen  minimum  zone 
(OMZ)  fluctuated  in  the  past.  Gorbarenko  (1996), 
in  the  Okhotsk  and  Bering  Seas,  found  two  episodes 
of  deglacial  freshening  of  surface  water  associated 
with  increases  in  %CaC03,  %opaI,  and  %organic 
carbon  (Corg).  However,  these  data  are  difficult  to 
correlate  with  the  global  observations  from  the  last 
glacial  maximum  (LGM)  and  deglaciation  because 
of  poor  age  control. 

In  this  paper,  we  will  discuss  our  observations 
of  the  deglaciation  in  the  Bering  Sea,  incor¬ 
porating  previously  unpublished  diatom  and  for¬ 
aminifer  <5i80  from  three  older  cores:  RC14-121 
(R.V.  Conrad ,  Umnak  Plateau),  VIN019-12GGC 
(R.V.  Vinogradov ,  Bowers  Ridge),  and  ANTP-42PC 
(R.V.  Melville ,  Shirshov  Ridge).  Our  synthesis  will 
show  that  the  changes  we  observe  were  widespread, 
not  only  within  the  Bering  Sea,  but  probably 
throughout  the  open  North  Pacific. 

2.  Bering  Sea  hydrography 

The  Bering  Sea  lies  between  the  Aleutian  Islands 
to  the  south,  Siberia  to  the  northwest,  and  Alaska  to 
the  northeast,  and  consists  of  a  broad  continental 
shelf  (50-1  50  m  deep)  and  a  deep  basin  (< 4000  m 


deep)  (Fig.  I).  The  principal  surface  circulation  in 
the  basin  is  a  cyclonic  gyre,  whose  western 
boundary  current  is  the  southward  flowing  Kam¬ 
chatka  Current.  Alaskan  Stream  water  flows  north¬ 
ward  through  passes  between  the  Aleutian  Islands 
and  is  incorporated  into  this  gyre  circulation.  The 
Bering  Slope  Current  (BSC),  which  flows  along  the 
edge  of  the  continental  shelf,  marks  the  eastern  edge 
of  the  gyre. 

On  the  shelf,  the  principal  surface  flow  is  north¬ 
ward,  with  net  northward  flow  through  the  Bering 
Strait,  which  is  85  km  wide  and  50  m  deep  (Schu¬ 
macher  and  Stabeno,  1998).  Water  from  the  Bering 
Sea  is  relatively  low  salinity  and  rich  in  nutrients, 
and  is  an  important  component  of  the  upper 
halocline  in  the  Arctic  (Cooper  et  al.,  1997).  With 
reduced  sea  level  (^130m  lower)  during  the  last 
glacial  period,  the  continental  shelf  was  subaerially 
exposed,  isolating  the  North  Pacific  from  the  Arctic 
and  North  Atlantic.  The  Bering  Strait  reopened 
during  the  deglaciation,  no  later  than  1 1,000  ,4Cyr 
before  present  (BP)  (latest  ~1 3,000  yr  BP)  (Elias 
et  al.,  1997). 

To  the  south,  the  Bering  Sea  is  connected  to  the 
North  Pacific  through  various  straits  and  passes,  the 
largest  of  which  are  (from  west  to  east)  Kamchatka 
Strait  (4420  m),  Near  Strait  (2000  m),  Buldir  Pass 
(640  m),  Amchitka  Strait  (1155  m),  and  Amutka 
Pass  (430  m)  (Fig.  1).  Physical  properties  of  water 
across  the  Aleutian  Islands  diverge  deeper  than 
2000  m,  indicating  relatively  unimpeded  exchange 
between  the  Bering  Sea  and  North  Pacific  at  and 
above  the  depth  of  Near  Strait.  Constricted 
exchange  of  deep  water  occurs  through  Kamchatka 
Strait.  In  the  eastern  Bering  Sea,  northward  flow 
through  the  Unimak  Pass  ( <  80  m  deep)  is  the  major 
conduit  between  the  North  Pacific  and  the  shelf 
(Stabeno  ct  al..  1999).  At  the  LGM,  this  pass  was 
closed,  but  exchange  through  the  remainder  of  the 
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Aleutian  Arc  was  probably  not  significantly  af¬ 
fected. 

NPIW  is  formed  by  cabbeling  of  the  northward 
flowing  warm  and  salty  Kuroshio  Current  and  the 
southward  flowing  cold  and  fresh  Oyashio  Current 
waters  where  they  meet  east  of  Hokkaido,  Japan 
(Yasuda,  1997;  You,  2003).  In  the  North  Pacific,  the 
core  of  NPIW  approximately  follows  the  26.8  oe 
potential  density  surface,  and  is  associated  with  a 
salinity  minimum  and  oxygen  maximum.  The 
density  of  the  core  of  NPIW  is  determined  by 
Okhotsk  Sea  Mode  Water  (OSMW),  the  product  of 
brine  rejection  during  winter  sea-ice  formation. 
OSMW  flows  out  of  the  Okhotsk  Sea  through  the 
Kuril  Islands,  and  mixes  with  WSAG  water, 
forming  Oyashio  Current  water.  Today,  NPIW  is 
found  at  200-400  m  in  the  Bering  Sea  (Macdonald 
et  al.,  2001).  The  water  mass  at  400-2000 m 
originates  from  mixing  of  NPIW  with  intermedi¬ 
ate-depth  water  that  originates  from  the  south 
(Roden,  2000).  Deeper  than  2000m,  compared  to 
the  North  Pacific,  Bering  Sea  water  is  slightly 
fresher,  warmer,  and  enriched  in  nutrients,  particu¬ 
larly  silicate,  which  dissolves  from  opal-rich  seafloor 
sediments  (Roden,  2000).  Consistent  with  the 
restricted  exchange  of  the  deepest  Bering  Sea  water 
with  the  North  Pacific,  A,4C  is  slightly  lower  than 
North  Pacific  water  of  equivalent  depth  (GEOSECS 
Executive  Committee  and  Ostlund,  1987). 

There  is  a  band  of  very  high  productivity  along 
the  BSC.  Low  stratification  along  the  slope  break  is 
caused  by  tidal  mixing  and  transverse  circulation, 
bringing  nutrients  from  300-800  m  depths  to  the 
surface  and  supporting  production  to  high  tropic 
levels  (Springer  et  al.,  1 996).  The  OMZ  in  the  Bering 
Sea  is  the  most  intense  at  ~900m  water  depth,  with 
concentrations  of  dissolved  oxygen  ([02])  as  low  as 
15pmol/kg  (Roden,  2000). 

3.  Methods 

For  stable  isotope  analysis,  0.5-cm-thick  (in 
laminated  intervals)  or  1-cm-thick  (in  massive 
intervals)  samples  with  dry  mass  of  1-5  g  were 
taken  every  12cm  in  the  HLY  cores.  Two-centi¬ 
meter-thick  samples  were  taken  every  4  cm  in 
VIN019-12GGC  and  every  5  cm  in  ANTP-42PC. 
Sediment  samples  were  oven  dried  at  50  °C,  then 
washed  with  tap  water  through  a  63-pm  sieve. 
Benthic  foraminifera  were  picked  from  the  >250- 
pm  fraction,  planktonic  from  the  150-250-pm 
fraction.  Single  tests  of  a  hispid  species  of  the 


Table  2 

Uncorrected  AMS  l4C  ages  measured  on  N.  pachyderma  (s.)  and 
corrected,  calibrated  ages 


Core 

Depth  (cm) 

l4C  yr  BP±<t 

yr  BP  ± 

H  LY02-02-3JPC 

140.5-142.5 

10,050  +  60 

10  "140 10,620 
^hj10  300 

HLY02-02-3JPC 

352.5-353.5 

10,850  +  65 

1I.650!J$S 

HLY02-02-3JPC 

600.5-601.5 

12,400  +  65 

13,500|»i» 

HLY02-02-3JPC 

906.5-909.5 

13,350  +  80 

14  4001 5,360 

^>WI4,330 

H  L  Y 02-02-5 1 J  PC 

134.5-135.5 

10,600  +  60 

1 1  250m,60° 

1  I ,  ^->^10.850 

HLY02-02  51JPC 

177.5-178.5 

12,500  +  60 

13  S4013,830 

i  j,  j*tu13  490 

HLY02-02-51 JPC 

241.5-242.5 

14,050  +  85 

1  5  Q2016,170 

13,  VZU15  680 

HLY02-02-17JPC 

155.5-156.5 

10,000  +  60 

10  32010,600 
iu,  ->^u1(U90 

HLY02-02-17JPC 

185.5-186.5 

12,550  +  65 

13  59Q 14-000 
i  j,  3^u]3  490 

HLY02-02-17JPC 

215.5-216.5 

13,400  +  60 

14  4  501 5,420 
Mi350 

Calibrations  performed  with  Calib4.4html  (Sluivcr  et  al.,  1998) 
using  A R  =  300 yr,  corresponding  to  a  ~700-yr  reservoir  correc¬ 
tion.  We  report  the  maximum  probability  calendar  ages  with  the 
upper  and  lower  limits  of  the  \a  range. 

benthic  genus  Uvigerina ,  and  eight  tests  of  the 
planktonic  species  Neogloboquadrina  pachyderma 
(sinistral)  were  analyzed  from  each  sample  on  a 
VG  PRISM-III  mass  spectrometer  with  an  Isocarb 
system  (90  °C  common  acid  bath)  for  <5lsO  with  a 
precision  of  ±0.07%o.  Exceptions  are  HLY02-02- 
51JPC,  where  a  costate  Uvigerina  species  was  used, 
and  in  VIN019-12GGC  and  ANTP-42PC  where 
15  N.  pachyderma  (s.)  comprised  each  analysis. 
Uvigerina  is  not  present  in  all  samples.  Percent 
CaC03  was  measured  as  in  Keigwin  et  al.  (1992). 
Methods  for  counting  diatom  species  relative 
abundance  are  described  in  Sancetta  et  al.  (1985). 
Diatom  <5lsO  data  and  a  lower  resolution  version  of 
the  diatom  relative  abundance  data  from  RC14-121 
were  previously  published  in  Sancetta  et  al.  (1985). 

Radiocarbon  dating  was  performed  at  the  Na¬ 
tional  Ocean  Sciences  Accelerator  Mass  Spectro¬ 
metry  (NOSAMS)  facility  at  WHOI  on  N. 
pachyderma  (s.)  >  150  pm  (Table  2).  In  calibrating 
our  dates,  we  assumed  the  local  anomaly  from  the 
mean  global  reservoir  correction  (A R)  was  constant 
at  300  yr.  The  modem  A R  in  the  North  Pacific  is  not 
well  known;  estimates  include  80-360  yr  in  the 
northern  and  eastern  Bering  Sea  (Dumond  and 
Griffn,  2002),  320±40yr  at  the  Kuril  Islands  and 
40  +  40  yr  in  the  Sea  of  Japan  (Kuzmin  et  al..  2002), 
and  240  +  50yr  near  Kamchatka  Strait  (Robinson 
and  Thompson.  1981).  It  is  very  likely  that  A R  is 
spatially  and  temporally  variable  as  a  function  of 
local  upwelling  strength  and  circulation,  so  our 
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Fig.  2.  Positive  X-radiograph  of  a  u-channel  sample  taken  from  HLY02-02-1GGC  (255-267cm).  a  companion  core  to  HLY02-02-3JPC. 
This  12  x  2.5cm  sample  is  from  a  depth  interval  that  is  the  equivalent  of  the  uppermost  laminated  interval  HLY02-02-3JPC.  From  our 
estimates  of  sediment  accumulation  rate,  each  light-dark  pair  probably  represents  one  year  of  accumulation.  Note  the  cyclic  variability  in 
lamina  thickness  indicating  sub-decadal  variability  in  the  seasonal  deposition  of  sediment. 


assumed  value  of  300  yr  should  only  be  regarded  as 
a  reasonable  guess.  A R  may  have  varied  by  up  to  a 
few  hundred  years  through  time  and  between  these 
sites,  but  this  would  not  affect  the  conclusions  of 
this  paper. 


4.  Results1 

4.1.  Lithology 

All  of  our  cores  contain  one  or  two  deglacial 
intervals  which  are  visually  identified  as  high 
%CaC03,  high  %opal,  and  high  %Corg  intervals 
with  a  distinctly  greenish  hue.  HLY02-02-3JPC 
contains  the  most  expanded  deglacial  sequence  of 
the  cores  in  this  study,  at  ~8m  long,  with  an 
average  sediment  accumulation  rate  of  190cm/kyr 
from  the  beginning  to  the  end  of  the  laminated 
sequences.  There  are  four  visibly  laminated  inter¬ 
vals,  which  are  composed  of  pairs  of  dark  and  light 
olive  laminae  ranging  in  thickness  from  <  1  to  2  mm 
(Fig.  2).  Intervening  massive  intervals  and  the 
deepest  6  m  of  sediment  are  composed  of  very  fine 
grained  sticky  homogeneous  dark  olive-gray  mud. 
The  topmost  1.4m  of  sediment  is  firm  dark-olive- 
gray  mud.  Holocene  sediments  from  the  Bowers  and 
Shirshov  Ridges  are  diatomaceous  and  contain 
moderately  high  CaC03.  The  sediment  accumula¬ 
tion  rates  in  HLY02-02-51 JPC  and  HLY02-02- 
17JPC  are  15-20cm/kyr  during  the  deglaciation. 
We  currently  have  no  age  control  for  the  Holocene 
portions  of  the  HLY  cores  because  of  the  rarity  of 
planktonic  foraminifera  and  the  unknown  amount 
of  sediment  lost  from  the  top  of  piston  cores  during 
the  coring  process. 

‘All  daia  reported  in  this  paper  are  posted  on  the  NOAA 
Paleoclimatology  Program  website  (\v\\w  .ncdc.noaa.gov  palco 
paleo.html). 


For  the  time  interval  corresponding  to  the 
deglaciation,  stratigraphic  correlation  between  the 
HLY  cores  is  based  on  ten  I4C  dates,  and  in  the 
other  cores  it  is  based  on  lithology.  The  upper  and 
lower  laminated  intervals  in  HLY02-02-51  JPC 
correspond  to  the  upper  one  and  lower  three 
laminated  intervals  in  HLY02-02-3JPC.  In  the  other 
four  cores,  the  analogs  to  the  laminated  sequence 
are  the  one  or  two  “green  layers,”  which  are 
indicated  with  shaded  bars  in  Figs.  3  and  4.  For 
the  cores  with  no  age  control,  this  correlation  is 
approximate,  assuming  that  green  and  laminated 
intervals  represent  contemporaneous  dysaerobic 
conditions  at  the  sediment-water  interface.  In 
RC 14-121,  three  uncorrected,  uncalibrated  14C 
dates  were  measured  on  20-cm-thick  sections  of 
sediment  and  were  interpreted  by  Sancetta  et  al. 
(1985)  as  maximum,  not  exact  ages,  establishing  the 
uppermost  lithological  unit  of  this  core  (0-1 75  cm) 
as  Holocene  (Table  2). 

4.2.  Oxygen  isotope  data 

There  are  two  rapid  decreases  of  planktonic  c5180 
in  HLY02-02-3JPC,  at  888-924  cm  and  348-372  cm 
(Fig.  3),  which  occur  at  14,400  and  ll,650yrBP, 
respectively.  Linearly  interpolating  between  the 
four  ,4C  dates,  the  transitions  occur  in  ^  1 90  yr 
and  ^  1 30  yr  (Fig.  3).  The  magnitude  of  the  <5180 
events  is  0.7-0. 8%o.  In  RC14-121  and  ANTP-42PC, 
the  <5,80  measurements  begin  probably  after  the 
start  of  the  deglaciation  since  the  highest  benthic 
<5i80  measurements  in  those  cores  are  lower  than 
the  LGM  values  in  similar-depth  cores  HLY02-02- 
17 JPC  and  VINO! 9-1 2GGC,  which  are  at  -5%o.  In 
RC  14-121,  there  is  a  ~l%o  decrease  in  benthic  <$,80 
in  adjacent  samples  at  170-190cm  (Fig.  4).  In  the 
rest  of  the  Bering  Sea  cores,  there  are  no  robust 
features  in  foraminifer  <5,80  in  the  data  currently 
available. 
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N.pachy  derma  (s.)  5180 
— Uvigerina  510O 


0  100  200  300  400 


depth  (cm) 

Fig.  3.  Shaded  bars  represent  “green”  layers,  which  are  characterized  by  high  %CaC03,  %Corg,  and  %opal.  In  (A)  and  (B),  these  layers 
are  laminated,  where  the  laminae  are  probably  annual  pairs  (see  text).  <$lsO  data  is  of  Uvigerina  with  costate  test  in  (B),  hispid  test  in  (C). 
The  l4C  dates  in  (A)  are  calibrated  to  calendar  years  (Table  2).  The  two  abrupt  ~0.7-0.8%o  decreases  in  <5I80  in  panel  (A)  occur  at  14,400 
and  ll,650yrBP,  each  event  lasting  at  most  ~200yr.  Assuming  direct  correlation  of  the  green  layers  between  cores,  the  sediment 
accumulation  rates  range  from  >  12cm/kyr  in  (B)  and  (C)  to  190cm/kyr  during  the  deglaciation  in  (A).  Note  that  the  x-axis  scales  are  not 
the  same  between  the  panels. 
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Fig.  4.  Shaded  bars  represent  “green’'  layers,  which  are  characterized  by  high  %CaC03,  %Corg,  and  %opal.  £180  data  is  of  Uvigerina 
with  hispid  test  in  (A),  (B),  and  (C).  The  thin  line  in  (C)  is  %CaC03,  in  (A)  it  is  #  planktonic  foraminifera  (PF)  per  gram  dry  sediment,  a 
proxy  for  %CaC03.  Presence  of  Nitzschia  is  associated  with  sea  ice,  R.  hebetata  is  associated  with  high  nutrient  availability  and  Corg  flux. 
Nitzschia  abundance  patterns  show  that  maximum  extent  of  sea  ice  may  not  have  reached  the  Bowers  or  Shirshov  Ridge  core  sites  during 
the  LGM  or  YD.  Note  that  the  x-axis  scales  are  not  the  same  between  the  panels.  The  >>-axis  scales  for  a  property  on  each  panel  are 
identical.  Previously  reported  in  Sancetta  et  al.  (1985)  were  the  diatom  <518 O  data  and  the  three  uncorrected,  uncalibrated  14C  dates  in  (B), 
which  were  measured  on  the  organic  fraction  of  20cm-thick  sections  of  sediment  and  interpreted  as  maximum,  not  exact  ages,  establishing 
the  uppermost  lithologic  unit  of  this  core  (0-1 75  cm)  as  Holocene. 


4.3.  Diatom  relative  abundance  data 

In  Umnak  Plateau  core  RC14-121  (Fig.  4),  the 
new  data  at  closer  sample  spacing  support  the 


earlier  observations  of  Sancetta  et  al.  (1985).  The 
sum  of  the  abundances  of  the  diatoms  Nitzschia 
grunowii  and  Nitzschia  cylindra  (both  sometimes 
assigned  to  Fragilariopsis)  indicate  the  presence  of 
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sea  ice.  During  the  late  glacial  or  early  deglacial 
transition,  these  sea-ice  species  occur  at  a  relative 
abundance  of  >20%.  There  is  a  subsequent 
decrease  in  Nitzschia  that  is  contemporaneous  with 
a  ~-9%o  decrease  in  bulk  opal  <5180.  This  is  followed 
by  a  newly  resolved  increase  in  Nitzschia ,  which 
occurs  as  the  Uvigerina  <5,80  decreases.  Nitzschia 
abundance  from  ANTP-42PC  of  Shirshov  Ridge  is 
low  throughout  the  deglaciation,  and  is  zero  in  each 
of  the  samples  at  VIN019-12GGC.  The  species 
Rhizosolenia  hebetata ,  which  indicates  high  nutrient 
availability  and  Corg  flux,  shows  highest  relative 
abundance  during  the  green  layers  in  ANTP-42PC. 
In  RC14-121,  R  hebetata  shows  no  significant 
changes,  with  an  average  relative  abundance  of 
~3%  through  the  deglaciation  and  Holocene. 

5.  Discussion 

5.1.  Chronology  and  correlations 

The  sequence  of  Bering  Sea  laminated  sediment  is 
similar  to  the  California  Borderland  Basins  and 
Gulf  of  California,  where  laminations  occur  during 
the  Holocene  and  the  Bolling-Allerod  (Keigwin  and 
Jones,  1990;  Behl  and  Kennett,  1996;  Zheng  et  ah, 
2000;  van  Geen  et  al.,  2003).  Based  on  a  simple  age 
model  where  we  linearly  interpolate  between  the 
four  14C  dates,  the  four  laminated  intervals  in 
HLY02-02-3JPC  may  correspond  to  (1)  the  Bolling, 
(2)  the  early  Allerod,  (3)  the  Allerod  following  the 
Inter-Allerod  Cold  Period,  and  (4)  the  early 
Holocene.  This  correlation  will  need  confirmation 
from  more  detailed  14C  dating  and  a  more  robust 
estimate  of  the  reservoir  correction.  From  our 
estimate  of  sediment  accumulation  rate  and  the 
lamina  thickness,  each  light-dark  pair  probably 
represents  one  year  of  accumulation. 

The  lithology  of  the  deglacial  sequence  in  the 
Bering  Sea  is  similar  to  that  observed  in  cores  from 
the  Okhotsk  Sea  and  Bering  Sea  (Gorbarenko, 
1996)  and  the  open  North  Pacific  (Keigwin  et  al., 
1992;  Gorbarenko,  1996;  Keigwin,  1998)  where 
during  the  deglaciation,  rapid  decreases  in  plank¬ 
tonic  foraminifer  <5180  precede  intervals  of  high 
%CaC03.  Percent  Corg  and  %opal  are  generally 
high  during  these  intervals  as  well.  From  dry  bulk 
density  measurements  and  an  estimate  of  sediment 
accumulation  rate,  Keigwin  et  al.  (1992)  determined 
that  not  only  the  percentage,  but  the  fluxes  of  the 
biogenic  components  increased.  This  conclusion  has 
received  further  support  from  accumulation  rates 


calculated  using  the  excess  Th  method  (Crusius 
ct  al..  2004). 

In  cores  from  the  far  Northwest  Pacific,  Okhotsk 
Sea,  and  Bering  Sea  Gorbarenko  (1996)  determined 
that  the  negative  planktonic  <5,sO  events  occurred  at 
12,500  and  9300  14Cyr  BP  using  a  1 100-yr  reservoir 
correction.  Calibrated  to  calendar  years  using 
A R  =  300  yr,  these  become  1 5,400  and  1 1 , 1 00  yr  BP, 
respectively.  Considering  the  difficulty  in  identifying 
the  location  of  the  <5lsO  events  in  cores  with  low 
sediment  accumulation  rates,  these  ages  are  reason¬ 
ably  close  to  our  dating  of  the  events  in  HLY02-02- 
3JPC  as  14,400  and  ll,650yrBP. 

The  deglacial  <5lsO  decrease  of  ~9%o  measured  in 
diatom  silica  in  RC14-121  is  much  larger  than  the 
foraminifer  <5,80  decreases  measured  in  any  of  the 
other  Bering  Sea  cores.  As  discussed  in  Sancetta 
et  al.  (1985),  diatom  <5lsO  is  more  sensitive  to 
changes  in  surface  hydrography  than  foraminifer 
<$180,  and  because  of  the  magnitude  of  the  event,  it 
probably  represents  a  freshening  of  surface  water.  If 
it  were  entirely  due  to  warming,  the  diatom  <5lgO 
signal  would  represent  a  temperature  increase  of 
~18°C  (calibration  of  Shemesh  et  al.  (1992)),  which 
is  not  probable.  Sancetta  et  al.  (1985)  hypothesized 
that  the  source  of  the  freshening  was  a  large  flux  of 
meltwater  from  retreating  alpine  glaciers  in  Alaska. 

Of  the  three  study  regions  with  diatom  data,  only 
the  Umnak  Plateau  (RC14-121)  has  indications  of 
significant  change  in  sea  ice  in  the  past,  with  two 
maxima  in  Nitzschia  at  what  appear  to  be  the  LGM 
and  the  YD.  The  Nitzschia  maximum  at  2.5  m  is 
interpreted  to  represent  extreme  glacial  conditions 
with  sea-ice  cover  throughout  the  year  (Sancetta 
et  al.,  1985).  Today,  the  maximum  winter  ice  edge  in 
the  Bering  Sea  does  not  reach  the  Bowers  Ridge, 
and  just  reaches  the  Umnak  Plateau  (Niebauer  et 
al.,  1999).  No  sea-ice  diatoms  are  present  in  any  of 
the  samples  analyzed  from  the  VIN019-12GGC 
(Bowers  Ridge).  Katsuki  et  al.  (2003)  report  that 
abundance  of  Thalassiosira  gravida  suggests  sea  ice 
was  proximate  at  the  crest  but  not  the  western  side 
of  the  Bowers  Ridge  at  the  LGM.  Nitzschia  is 
present  at  Shirshov  Ridge  core  ANTP-42PC  but  in 
low  abundance  throughout  the  record.  Percent  R. 
hebetata  is  relatively  high  in  the  older  green  layer 
and  at  the  base  of  the  younger  green  layer, 
suggesting  high  Corg  flux  to  the  sediments  when 
these  layers  were  deposited.  Numbers  of  planktonic 
foraminifera  per  gram  dry  sediment  (a  proxy  for 
%CaC03)  reach  a  peak  in  both  the  green  layers  at 
1-2  orders  of  magnitude  higher  than  in  intervening 
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sediments,  so  probably  represent  high  CaC03  flux 
at  those  times  as  well. 

5.2.  What  produces  the  laminated  sediment? 

Laminated  sediments  occur  when  the  dissolved 
oxygen  at  the  sediment-water  interface  is  too  low  for 
burrowing  macrofauna  ([02]<5jiM).  This  can  be 
caused  by  the  intersection  of  an  intense  OMZ  with 
the  seafloor.  It  also  can  occur  when  the  flux  of  Corg  is 
so  high  that  the  oxygen  at  the  bottom  of  the  water 
column  is  consumed  by  respiration.  A  third  mechan¬ 
ism  is  seasonal  fall-out  of  mat-forming  diatoms  or 
bacteria  which  can  form  a  thick  enough  layer  to 
suffocate  benthic  macrofauna  (see  Kemp,  1996). 

In  cores  from  the  low-latitude  eastern  North 
Pacific  during  the  deglaciation,  there  is  an  increase 
in  the  extent  of  laminations  from  north  to  south 
(Zheng  et  ai.,  2000),  toward  the  low-oxygen  inter¬ 
mediate  water  end  member  from  the  equatorial 
upwelling  region.  Stott  et  al.  (2000)  and  van  Geen 
et  al.  (2003)  document  asynchrony  of  laminations  in 
cores  from  the  Pacific  coast  of  California  and  Baja, 
highlighting  the  influence  of  local  export  production 
in  the  occurrence  of  laminated  sediment.  In  the 
Bering  Sea,  despite  its  proximity  to  the  source  of 
well-oxygenated  NPIW,  the  belt  of  high  productiv¬ 
ity  in  the  BSC  appears  to  have  drawn  [OJ  low 
enough  to  bring  it  near  the  5  pM  threshold  for 
laminae  preservation. 

A  simple  explanation  for  the  contemporaneous 
sequences  of  laminated  and  dysaerobic  sediments  in 
widely  separated  regions  of  the  North  Pacific  is  that 
the  primary  control  is  the  [02]  of  the  large-scale 
intermediate  water  mass  and  that  flux  of  NPIW  was 
reduced  in  the  YD  and  other  cold  deglacial  climate 
events.  But  anoxia  caused  by  an  increase  in  export 
production  leaves  a  similar  geochemical  imprint  in 
the  sedimentary  record,  and  there  is  evidence  of 
higher  fluxes  of  biogenic  sedimentary  components 
during  the  deglacial  events.  Based  on  data  available 
at  present,  we  cannot  distinguish  between  the  two 
hypotheses.  It  is  likely  that  both  processes  are 
involved. 

6.  Conclusions 

(1)  The  deglacial  lithology  in  the  Bering  Sea  can  be 
correlated  to  other  parts  of  the  North  Pacific. 
The  deglaciation  is  characterized  by  one  to  two 
“green  layers,”  which  are  occasionally  lami¬ 
nated. 


(2)  There  is  a  widespread  occurrence  of  laminated 
sediments  in  the  North  Pacific,  from  the  well- 
documented  California  margin  and  Gulf  of 
California  anoxic  basins  to  the  continental  slope 
of  the  Bering  Sea.  In  one  of  our  cores  from  the 
continental  slope  near  the  Pervenets  Canyon, 
the  laminations  occur  during  the  Bolling-Aller- 
od  and  the  early  Holocene,  which  is  similar  to 
Behl  and  Kcnnett's  (1996)  findings  in  Santa 
Barbara  Basin. 

(3)  The  primary  control  on  the  presence  of  lami¬ 
nated  sediments  in  the  North  Pacific  could  be 
the  oxygen  content  of  NPIW.  But  both  ventila¬ 
tion  and  export  production  probably  contribute 
to  the  dysaerobic  conditions. 

(4)  Two  abrupt  freshening  and/or  warming  events 
are  recorded  in  planktonic  foraminifer  <5I80  at 
14,400  and  1 1,650  yr  BP,  which  directly  precede 
the  “green  layers.” 

(5)  There  is  a  YD-like  cold  event  in  the  Bering  Sea 
region,  which  is  expressed  as  a  maximum  of  sea- 
ice  diatoms  Nitzschia  in  a  core  from  the  Umnak 
Plateau. 
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Chapter  3 


The  deglaciation  in  the  Bering  Sea 


Abstract 

Cores  from  within  the  oxygen  minimum  zone  in  the  Bering  Sea  were  laminated  during 
the  B0lling-Aller0d  (BA,  14.7-12.9  ky  BP),  a  phenomenon  which  has  been  observed  in 
cores  from  around  much  of  the  rim  of  the  North  Pacific,  from  the  Japanese  margin  to  the 
California  Margin  and  Gulf  of  California.  A  compilation  of  published  and  unpublished 
14C  measured  in  pairs  of  benthic  and  planktonic  foraminifera  from  the  North  Pacific  shows 
no  consistent  change  in  the  14 C  of  subsurface  water  between  the  BA  and  the  Younger 
Dryas  (YD,  12.9-11.7  ky  BP)  or  Last  Glacial  Maximum  (LGM,  23-21  ky  BP).  The  relative 
timing  of  the  lithological  transitions  in  a  partially-laminated,  high  sediment-accumulation- 
rate  (200cm/ky)  core  from  1132m  in  the  Northern  Bering  Sea  is  very  close  to  the  timing 
of  the  centennial-  and  millennial-scale  climate  events  in  Greenland  ice  core  (GISP2) 
to  well  within  the  uncertainty  of  the  age  model.  I  do  not  know  the  absolute  timing  of  these 
events  between  the  Bering  Sea  and  GISP2  because  the  surface  reservoir  correction  (A/?)is 
not  well  constrained.  However,  the  match  of  the  relative  timing  suggests  that  A R  was  fairly 
constant  through  the  deglaciation.  The  two  rapid  decreases  in  N.  pachyderma  (s.)  5180  at 
the  beginning  of  the  B0lling  and  the  end  of  the  YD  in  this  core  (Cook  et  a!.,  2005)  are 
accompanied  by  large,  but  brief  decreases  in  the  relative  abundance  of  N.  pachyderma  (s.). 
This  suggests  that  the  6lsO  reflects  transient  warmings,  but  sustained  freshening  of  surface 
water. 


3.1  Introduction 

The  heat  transport  associated  with  the  Atlantic  meridional  overturning  circulation  (AMOC) 
is  considered  one  of  the  most  important  components  of  the  climate  system.  Because  of  the 
sheer  breadth  of  the  Pacific,  the  poleward  heat  transport  associated  with  the  wind-driven 
shallow  gyre  circulation  at  24°N  is  two-thirds  the  magnitude  of  the  heat  transport  associated 
with  the  AMOC  at  the  same  latitude  (Talley,  2003).  The  Pacific  contains  the  ocean’s  largest 
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Figure  3-1:  The  scale  for  color  bathymetry  and  topography  is  in  meters.  Coring  sites 
from  HLY02-02  are  labeled  with  black  dots  and  the  site  number.  See  Table  3.1  for  details. 
WOCE  stations  are  in  yellow.  Hydrocasts  from  the  HLY02-02  cruise  are  in  blue  stars. 
Sediment  traps  discussed  in  text  are  purple  triangles.  On  the  inset  map,  red  circles/ triangles 
denote  North  Pacific  coring  sites  where  dysoxic/laminated  sediments  are  found  during  the 
deglaciation  (Keigwin  and  Jones  (1990);  Behl  and  Kennett  (1996):  Zheng  et  ah  (2000): 
Marret  et  al.  (2001):  van  Geen  et  al.  (2003);  Ortiz  et  ah  (2004):  Crusius  et  ah  (2004): 
Narita  et  al.  (2004):  and  this  study). 
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reservoir  of  deep  water  and  dissolved  inorganic  carbon,  so  a  significant  change  in  circulation 
in  the  Pacific  could  potentially  have  a  large  impact  on  global  climate.  Since  relatively  few 
high-quality  Pleistocene  paleoceanographic  records  exist  from  the  open  Pacific,  it  is  not 
well  understood  what  changes  occurred  in  the  Pacific  during  the  deglaciation,  when  large 
reorganizations  occurred  in  the  North  Atlantic.  In  the  western  North  Atlantic,  sedimen¬ 
tary  [231Pa]/[230Th]  approaches  the  production  ratio  of  0.093  during  Heinrich  Event  1  (HI, 
~17.5ky  BP)  and  the  Younger  Dryas  (YD,  12.9-11.6  ky  BP),  suggesting  that  the  AMOC 
was  very  weak  or  shut  down  during  these  cold  climate  events  (McManus  et  ah,  2004).  For 
these  two  periods,  passive  water  mass  tracers  (513C  and  Cd/Ca)  suggest  that  there  was  a 
greater  contribution  of  southern  source  intermediate  and  deep  waters  in  the  North  Atlantic 
(Boyle  and  Keigwin,  1987;  Zahn  et  al.,  1997;  Vidal  et  al.,  1997).  Some  investigators  have 
suggested  that  the  circulations  in  the  North  Atlantic  and  North  Pacific  are  out  of  phase, 
where  the  AMOC  is  strong  when  the  rate  of  formation  of  North  Pacific  Intermediate  Water 
(NPIW)  is  weak,  and  vice  versa  (Kiefer  et  ah,  2001).  The  geological  records  from  the  North 
Pacific  are  suggestive,  but  not  conclusive. 

Ventilation  of  the  North  Pacific  during  the  deglaciation 

In  the  North  Pacific,  there  was  a  widespread  occurrence  of  dysoxic  and  laminated  sediment 
during  the  deglaciation  at  oxygen-minimum-zone  (OMZ)  depths  (~800-1200  m),  in  partic¬ 
ular,  during  the  Bplling-Allerpd  (BA,  14,89-12.89  ky  BP,  citetStuiM95-f ).  Observations  of 
dysoxic  and  laminated  sediments  have  been  made  across  most  of  the  rim  of  the  Pacific, 
including  the  Mexican  margin  (Hendy  et  al.,  2003),  the  Gulf  of  California  (van  Geen  et  al., 
2003;  Keigwin  and  Jones,  1990),  the  California  borderland  basins  (Behl  and  Kennett,  1996), 
off  the  Pacific  northwest  (McKay  et  al.,  2004),  the  Gulf  of  Alaska  (Marret  et  al.,  2001),  the 
Bering  Sea  (Cook  et  al.,  2005),  and  the  Japanese  margin  (Narita  et  al.,  2004)  (see  Figure  3- 
1).  There  is  yet  no  consensus  on  the  primary  mechanism  for  this  phenomenon,  although 
the  apparent  basin-wide  extent  of  reduction  of  the  [O2]  in  the  OMZ  suggests  a  basin-wide 
change  in  chemistry. 

The  position  of  the  minimum  in  oxygen  in  the  water  column  is  a  function  of  circulation 
and  the  consumption  of  oxygen  by  remineralization  of  the  organic  matter  which  originates 
from  the  surface  ocean  (Wyrtki,  1962).  There  are  several  mechanisms  which  can  lower  [O2] 
of  intermediate-depth  Pacific  water  including:  (1)  decrease  in  the  formation  rate  of  02-rich 
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North  Pacific  intermediate  water  (NPIW),  (2)  decrease  in  the  [O2]  of  newly- formed  NPIW 
without  a  change  in  ventilation  rate  (Crusius  et  ah,  2004),  (3)  intensification  of  the  OMZ 
from  greater  C org  export  from  the  surface  ocean  and  subsequent  remineralization  at  depth 
(McKay  et  ah,  2004). 

If  the  formation  rate  of  NPIW  is  entirely  responsible  for  the  changes  in  lithology,  then 
this  would  support  the  hypothesis  that  the  anti-phase  relationship  of  the  North  Atlantic 
and  North  Pacific  circulations.  One  way  to  test  this  hypothesis  is  by  reconstructing  the 
14C  of  the  subsurface  North  Pacific  using  pairs  of  benthic  and  planktonic  foraminifer  14C 
measurements.  All  else  being  equal,  if  NPIW  grew  weaker,  then  I  would  expect  a  decrease 
in  14C  content  of  intermediate-depth  water-that  is,  a  decrease  of  benthic  14C  relative  to 
planktonic  14C.  There  are  only  a  few  time  series  of  pairs  of  benthic-planktonic  14C  from 
the  LGM  through  the  deglaciation:  from  the  Japan  margin  (Ahagon  et  ah,  2003;  Duplessy 
et  ah,  1989),  the  Santa  Barbara  Basin  (Ingram  and  Kennett,  1995),  and  the  northeast 
Pacific  near  the  California-Oregon  border  (Mix  et  ah,  1999). 

The  core  of  NPIW  follows  the  26.8  a$  potential  density  surface  (Yasuda,  1997;  You, 
2003),  which  is  shallow  (200-400  m)  in  the  modern  Bering  Sea  and  the  Gulf  of  Alaska 
(Macdonald  et  ah,  2001).  Macdonald  et  ah  (2001)  note  that  in  the  Gulf  of  Alaska,  the 
26.8  oq  isopycnal  surface  may  outcrop  for  short  periods  in  the  winter  and  ventilate  the 
intermediate  depths.  This  is  supported  by  hydrographic  and  tritium  data  and  potential 
vorticity  calculations  of  Van  Scoy  et  ah  (1991).  The  26.8  oq  isopycnal  surface  is  also  shallow 
ing  the  Bering  Sea.  The  presence  of  CFC’s  in  the  abyssal  Bering  Sea  Basin  (Warner  and 
Roden,  1995),  suggests  a  small  amount  of  deep  water  formation  within  the  last  40  years, 
probably  from  brine  rejection  during  sea-ice  formation.  Ohkushi  et  ah  (2003)  suggest  from 
radiolarian  data  that  the  Bering  Sea  was  a  significant  locus  of  intermediate-depth  ventilation 
during  the  glacial  period,  when  sea  surface  temperatures  were  colder. 

Deglacial  SST  change  in  the  North  Pacific 

There  is  sparse  marine  evidence  of  a  YD  cooling  in  the  open  Pacific,  an  event  with  a  small 
enough  duration  to  perhaps  be  obscured  by  the  relatively  low  sediment  accumulation  rates 
in  the  Pacific  basin  and  poor  preservation  of  calcite.  Kiefer  and  Kienast  (2005)  review 
published  Pacific  SST  proxy  records  and  conclude  that  the  spatial  pattern  of  deglacial 
warming  was  heterogeneous,  with  Younger-Dryas-like  deglacial  cold  intervals  observed  in 
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location 

core 

depth 

latitude/longitude 

Bering  Slope 

HLY02-02-3JPC 

1132m 

N60.1279VE180.55820 

Bowers  Ridge 

HLY02-02-17JPC 

2209  m 

N53.9330°/E178.6988° 

Bowers  Ridge 

HLY02-02-25GGC 

721  m 

N53.9963°/E179.5842° 

Bowers  Ridge 

HLY02-02-18GGC 

856  m 

N53.60997E178.7225” 

Umnak  Plateau 

HLY02-02-51JPC 

1467  m 

N54.55327E191.33310 

Umnak  Plateau 

HLY02-02-55JPC 

1005  m 

N54.3991°/E192.0861° 

Umnak  Plateau 

HLY02-02-48JPC 

1960  m 

N54. 791 1°/E189. 5997° 

Table  3.1:  The  cores  used  in  this  study.  See  Figure  3-1  for  core  locations. 


marginal  basins  but  not  in  the  open  Pacific,  and  commencement  of  deglacial  surface  warming 
beginning  after  HI.  Lacustrine  records  of  pollen  (Peteet  and  Mann,  1994;  Ager,  2003)  and 
<5180  of  biogenic  silica  (a  temperature  proxy,  Hu  and  Shemesh  (2003))  from  Alaska  suggest 
a  warm  and  moist  BA,  and  a  relatively  cool  and  dry  YD.  This  inference  is  supported  by  a 
glacier  advance  in  southwest  Alaska  dated  to  the  YD  (Briner  et  al.,  2002). 

The  purpose  of  this  study  is  to  address  some  of  the  outstanding  questions  about  the 
deglacial  North  Pacific  using  new  sediment  cores  collected  during  a  cruise  to  the  Bering  Sea 
in  June,  2002,  and  with  a  compilation  of  published  and  unpublished  North  Pacific  data. 
In  particular,  how  did  the  14 C  of  subsurface  water  change  during  deglacial  climate  events? 
And  how  did  the  surface  water  of  the  North  Pacific  change  during  the  deglaciation? 


3.2  Methods 

The  cores  in  this  study  were  collected  in  June,  2002,  from  the  Bering  Sea  on  HLY02-02  (see 
Figure  3-1  and  Table  3.1  for  core  locations  and  depths).  Samples  were  prepared  as  described 
in  Cook  et  al.  (2005).  Briefly,  sediment  samples  were  dried  at  50°C,  then  washed  with  tap 
water  through  a  63  (im  sieve.  N eogloboquadrina  pachyderma  (sinistral)  and  Uvigerina  sp. 
were  picked  from  the  212-250  pm  and  >250  pun  size  fractions,  respectively.  The  planktonic 
fauna  >150 /im  were  counted  and  identified,  and  I  calculated  the  relative  abundance  of 
N.  pachyderma  (s.)  as  the  number  of  that  species  divided  by  the  total  number  of  individuals 
counted  in  the  sample.  Confidence  limits  on  the  relative  abundance  of  N.  pachyderma  (s.) 
were  calculated  at  the  95%  level  as  in  Fatela  and  Taborda  (2002). 

Isotope  measurements  on  8-10  N.  pachyderma  (s.)  or  a  single  Uvigerina  sp.  were  made 
on  a  PRISM  mass  spectrometer  with  a  common  acid  bath  at  90  °C,  and  are  reported  relative 
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depth  (cm) 

species 

14C  y  BP±cr 

cal  y  BP  a 

±2  <j 

HLY02-02-3JPC 

140.5-142.5 

N.  pachy derma  (s.) 

10,050±60 

10,510±5q 

i  130 
±150 

352.5-353.5 

N.  pachyderma  (s.) 

10,850±65 

11,550±}|q 

.320 

^270 

600.5-601.5 

N.  pachyderma  (s.) 

12,400±65 

13,440±foo 

.180 

^150 

839.5-840.5 

N.  pachyderma  (s.) 

13,300±70 

14,600±250 

.370 

^390 

839.5-840.5 

G.  affinis 

14,250±65 

906.5-909.5 

N.  pachyderma  (s.) 

13,350±80 

14,700±2jo 

_l350 

±460 

1427.5-1428.5 

N.  pachyderma  (s.) 

18,100dbl30 

20,430±i9o 

i  410 
^340 

HLY02-02-51JPC 

134.5-135.5 

N .  pachyderma  (s.) 

10,600±60 

11,190±4q 

.110 

■^80 

177.5-178.5 

N.  pachyderma  (s.) 

12,500±60 

13,550±9o 

.160 

150 

241.5-242.5 

N.  pachyderma  (s.) 

14,050±85 

15,700±220 

.420 
^  390 

HLY02-02-17JPC 

155.5-156.5 

N.  pachyderma  (s.) 

10,000±60 

10,460±go 

_L  110 

^  170 

185.5-186.5 

N.  pachy deima  (s.) 

12,550±65 

13,600±ioo 

.150 

*170 

214.5-215.5 

N.  pachyderma  (s.) 

13,600±50 

15,110±i§8 

.260 

*210 

214.5-215.5 

N.  labradorica 

14,200±110 

215.5-216.5 

N.  pachyderma  (s.) 

13,400±60 

14,830±1650 

i_300 

*390 

383.5-384.5 

N.  pachyderma  (s.) 

25,200±130 

28,500± 

± 

Table  3.2:  Uncorrected  AMS  14C  ages  and  corrected,  calibrated  ages.  Calibrations  per¬ 
formed  with  CalibS.0.1  (Reimer  et  al.,  2004)  using  the  Marine04  calibration  (Hughen  et  al., 
2004b).  I  use  AR=400y,  corresponding  to  a  reservoir  correction  of  800  y.  I  report  the 
median  calendar  ages  with  the  upper  and  lower  limits  of  the  a  and  2 a  ranges. 


core 

species 

Wc 

5180 

HLY-02-02-MC02 

Cibicidoides  tenmmoargo 

-0.825 

2.143 

HLY-02-02-MC02 

Uvigerina  peregrina 

-0.902 

3.419 

HLY-02-02-MC05 

Cibicidoides  spp.A 

-0.390 

2.804 

HLY-02-02-MC05 

Cibicidoides  spp.A 

-0.388 

2.708 

HLY-02-02-MC05 

Uvigerina  canariensis 

-0.764 

3.545 

HLY-02-02-MC05 

Uvigerina  peregrina 

-0.770 

3.382 

HLY-02-02-MC16 

Cibicidoides  mundulus 

-0.394 

2.991 

HLY-02-02-MC16 

Cibicidoides  mundulus 

-0.414 

2.964 

HLY-02-02-MC19 

Uvigerina  peregrina 

-1.237 

4.108 

HLY-02-02-MC24 

Uvigerina  peregrina 

-1.134 

3.434 

HLY-02-02-MC31 

Cibicidoides  kullenberg 

-0.052 

2.867 

HLY-02-02-MC43 

Uvigerina  peregrina 

-1.288 

3.665 

HLY-02-02-MC50 

Uvigerina  peregrina 

-1.149 

3.501 

HLY-02-02-MC56 

Uvigerina  peregrina 

-0.991 

3.210 

HLY-02-02-MC56 

Uvigerina  peregrina 

-1.107 

3.219 

Table  3.3:  Stable  isotope  ratios  from  rose-Bengal  stained  Cibicidoides  and  Uvigerina  from 
sediments  from  the  upper-most  (0-2  cm)  sediments  from  multicores,  which  capture  the 
sediment- water  interface  (unpublished  data  courtesy  of  D.  McCorkle).  These  data  are 
plotted  in  Figure  3-2. 
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to  PDB  via  NBS-19  (Coplen,  1996).  Radiocarbon  ages  in  Table  3.2  were  measured  at  the 
National  Ocean  Sciences  Accelerator  Mass  Spectrometer  Facility  at  WHOI.  I  converted 
radiocarbon  ages  to  calendar  years  with  CalibS.0.1,  using  A R  =  400  y,  so  that  the  reservoir 
age  of  surface  water  is  assumed  to  be  constant  at  R  =  800  y.  I  report  the  median  of  the 
probability  density  function  of  the  calibrated  calendar  age  (Telford  et  al.,  2004).  A R  is  the 
local  anomaly  of  the  apparent  age  of  surface  water  from  the  global  mean,  and  is  poorly 
constrained  in  the  Bering  Sea  (see  discussion  in  Cook  et  al.  (2005)).  Reported  uncertainties 
in  Table  3.2  do  not  take  into  account  uncertainty  in  ( R  +  A R).  For  the  radiocarbon  ages 
older  than  25  ky,  the  reported  calendar  age  is  the  intersection  of  the  radiocarbon  age  with 
the  Globigerina  bulloides  radiocarbon  record  from  the  Cariaco  Basin  (Hughen  et  al.,  2004a) 
on  the  ss09sea  age  model  of  the  GRIP  ice  core  (Johnsen  et  al.,  2001).  Age  models  for  each 
core  are  simply  constructed  as  linear  interpolation  of  median  calibrated  ages. 


Equilibrium  <$18OcaiCjte  calculation 

In  experiments  performed  by  Kim  and  O’Neil  (1997),  where  calcite  precipitated  from  solu¬ 
tions  of  NaHCOa  and  CaCl2  in  deionized  water,  the  fractionation  factor,  o^(C8dc\te-H20)y  was 
determined  as  an  exponential  function  of  T~1, 


In  O(caicjte_jj2o)  —  AT  +  B 


where 

_  (  16Q  )calcite  _  1000  +  <5180caCQ3(PDB)  _(  ~re0)pDB 

calc.te  H20  (“g)^  1000  +  518C>h20(SMOW)  ({§§)SMOW 


and 


— i8q - =  1.03092. 

(tsq)smow 

In  this  paper,  all  518Ocaicite  is  reported  relative  to  the  PDB  standard  through  NBS-19 
(Coplen,  1996),  and  5180h2O  ami  5180sw  are  relative  to  the  SMOW  standard.  Kim  and 
O’Neil  (1997)  performed  experiments  with  a  range  of  T  and  [Ca]G,  the  initial  concentration 
of  calcium  in  the  solution,  where  aCaicite-H20  increased  with  greater  T  and  with  greater 
[Ca]Q.  The  ratio  of  the  concentrations  of  major  ions  in  seawater  is  fairly  uniform  in  the 
ocean  (variation  is  on  the  order  of  1%);  at  a  salinity  of  35,  [Ca]  =  10.5  mM  (Wilson,  1975).  I 
approximated  the  relationship  between  aca|cite_H2o  and  T  for  [Ca]G  =  10  mM  by  calculating 
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Figure  3-2:  Compilation  of  (<518Ocaicjte(PDB)  -  <518Osw(SMOW))  and  temperature  (T) 
(modified  from  Lynch-Stieglitz  et  al.  (1999)).  Planktonic  foraminifer  cultures  are  of  Glo- 
bigerinoides  sacculifer  (Erez  and  Luz,  1983)  and  Orbuhna  universa  (Bemis  et  al.,  1998). 
Mollusk  data  include  both  cultured  and  field  samples  of  a  variety  of  genera  and  species 
(Epstein  et  al.,  1953).  Cibicidoides  and  Planulina  species  are  multicore-top  sediments  from 
Little  Bahama  Bank  (Slowey  and  Curry,  1995;  Curry  et  ah,  1993).  HLY02-02  data  are 
Uvigerinci  and  Cibicidoides  individuals  stained  with  Rose  Bengal  from  multicore-top  sedi¬ 
ments  (unpublished  data  courtesy  of  D.  McCorkle,  see  Table  3.3).  All  data  are  corrected, 
if  necessary,  so  that  the  standards  for  £l80ca|Cite  and  <P8Osw  are  PDB  and  SMOW,  re¬ 
spectively.  The  thin  lines  are  linear  least-squares  best  fits  to  each  set  of  planktonic  culture 
data  and  the  Little  Bahama  Bank  data.  The  thick  grey  dashed  line  is  the  estimate  of  the 
theoretical  exponential  relationship  between  (518OcaiCite(PDB)  -  £180sw(SM0W))  and  T 
where  [Ca]0  =  10  mM  (based  on  the  inorganic  calcite  precipitation  experiments  of  Kim  and 
O’Neil  (1997))  (Equation  3.1). 
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S  (PSU) 

Figure  3-3:  HLY02-02  hydrocast  data.  See  Figure  3-1  for  hydrocast  locations.  Unpublished 
data  courtesy  of  D.  McCorkle.  The  thick  grey  dashed  line  is  a  linear  least-squares  fit  to  the 
HLY02-02  data  (Equation  3.3). 
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the  mean  parameters,  A  and  B ,  for  the  experiments  at  5mM  and  15  mM  (Figure  3-2).  The 
resulting  expression  for  (£18Ocaicite-£18OH2o)  as  a  function  of  T  (in  Kelvin)  is 


<518Oca]Cite  -  <5180h2o 


-1000  + 


1000  +  5180h2o 
1.03092 


exp 


/ 19.87 

V  T 


-  38.09  x  10 


(3.1) 


The  difference  is  relatively  insensitive  to  <5180h2o-  For  a  10 %o  increase  in  £18Oh2o>  which 
is  a  larger  range  than  most  measured  values  in  the  modern  ocean,  and  for  T  =  0-40°C, 
(<5180caicite — <518Oh2o)  increases  by  an  average  of  0.006  %o.  This  is  less  than  measurement 
precision  of  6180.  In  Figure  3-2,  the  <5180ca)cite  curve  is  calculated  for  3180h2o  =  0%o. 

It  is  generally  considered  sufficient  to  estimate  the  relationship  between  (518Ocaicite--518Osw) 
and  T  with  a  first-  or  second-order  polynomial  because  the  scatter  in  field  and  laboratory 
data  is  large  enough  that  use  of  an  exponential  relationship  does  not  significantly  improve 
goodness  of  fit.  However,  at  temperatures  lower  than  10&C,  beyond  the  range  of  most  in¬ 
organic  and  planktonic  foraminifer  culturing  experiments,  the  polynomial  and  exponential 
curves  diverge.  Benthic  foraminifera  Cibicidoides  and  Planuhna  (Lynch-Stieglitz  et  ah, 
1999)  from  the  Bahamas,  and  Cibicidoides  from  the  Bering  Sea  fall  on  the  extrapolation 
of  a  linear  fit  through  the  planktonic  foraminifer  culture  data  (Bemis  et  ah,  1998;  Erez 
and  Luz,  1983)  (Figure  3-2).  The  quantity  [<518Oca]cite— <5180sw]  of  field  samples  of  mol- 
lusks  (Epstein  et  ah,  1953)  and  Uvigerina  (McCorkle,  unpublished)  are  higher,  and  fall 
on  the  theoretical  exponential  curve.  I  chose  to  use  the  exponential  relationship  to  calcu¬ 
late  equilibrium  <hl8Ocaicite  in  this  paper  because  it  is  the  relationship  I  would  expect  from 
an  Arrhenius  behavior  of  c*(caicite-H20)*  Using  this  relationship,  Uvigerina  sp.  calcifies  in 
equilibrium,  as  postulated  by  Shackleton  (1974). 

When  estimating  a  change  in  T  from  two  measurements  of  £18Ocaicite>  both  values  of 
must  be  corrected  for  the  whole-ocean  change  in  £18Osw  due  to  the  volume  of  continental 
ice  (Waelbroeck  et  ah,  2002).  After  this  correction,  I  may  set  <5180h2O  to  0%o  because  the 
absolute  value  of  5180h2o  does  not  significantly  affect  AT.  The  expression  for  AT  becomes 


t2-tx  = 


19.87 


19.87 


38.09  x  lO-3  +  In ( ^ 1  ”  ooo/ll 03092°° )  38.09  x  10~3 


i  Inf  ^18^calcitel  +  1000  \ 
^  “H  1000/1  03092  ) 


(3.2) 


Berner  (1975)  conducted  inorganic  calcite  precipitation  experiments  in  seawater  and 
found  that  the  precipitates  contained  7-10  mol%  Mg,  however  he  did  not  measure  the 
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5IS O  of  these  precipitates.  Tarutani  et  al.  (1969)  determined  that  the  amount  of  Mg  in 
inorganically  precipitated  calcite  is  related  to  the  8IS O,  where  <5180  increases  by  0.078  %o  for 
an  increase  in  1  mol%  Mg  of  the  calcite.  The  [Mg]/[Ca]  of  a  typical  Quaternary  foraminifer 
is  1-10  mmol/mol  (Rosenthal  et  al.,  1999),  so  this  effect  is  probably  not  significant.  I  will 
make  the  assumption  that  Equation  3.2  is  also  applicable  to  calcite  precipitation  in  seawater, 
and  will  substitute  5180sw  f°r  ^18Oh20  in  Equation  3.2. 

I  calculated  a  linear  least-squares  fit  to  hydrocast  data  collected  during  the  HLY02-02 
cruise  (Figure  3-3)  (McCorkle,  unpublished)  to  find  this  relationship  between  5  and  £18Osw 

<S180Sw  =  -12.5  +  0.3575  (n  =  142;r2  =  0.9657;p  <  1(T6).  (3.3) 


This  relationship  is  almost  identical  to  the  one  calculated  for  Okhotsk  Sea  hydrocasts  and 
GEOSECS  data  in  Keigwin  (1998).  Combining  equations  3.1  and  3.3,  we  can  calculate 
equilibrium  6180caicite,  given  T  and  S: 


$  Ocaicite  =  —12.5  +  0.3575  -  1000  T 


1000  -  12.5  +  0.3575 
1.03092 


exp 


^  1U87  _  3g  og  x  10-3^ 

(3-4) 


In  Figure  3-4c  is  equilibrium  <518Ocaidte  calculated  with  Equation  3.4  and  T  and  5  from 
the  WOCE  dataset  (refer  to  Figure  3-1  for  WOCE  station  locations).  The  WOCE  transect 
was  conducted  in  Summer,  1993.  In  Figures  3-4a-c,  I  use  Levitus  et  al.  (1998)  climatology 
from  the  nearest  gridbox  to  calculate  the  equilibrium  518Ocakite  m  the  top  400  m  of  the 
water  column  throughout  the  year  for  three  Bering  Sea  coring  sites  to  see  how  much  the 
equilibrium  <5lsO  in  the  surface  ocean  changes  seasonally. 


Planktonic  fauna  and  SST 

The  two  predominant  species  in  the  Bering  Sea  during  the  last  glacial  cycle  are  N.  pachy- 
derma  (s.)  and  G.  bulloides .  N.  pachyderma  (s.)  is  the  dominant  species  in  the  polar  regions. 
Coretop  planktonic-foraminifer  faunal  assemblages  are  often  calibrated  against  SST  in  order 
to  reconstruct  SST  in  the  past  (Imbrie  and  Kipp,  1971).  The  key  assumption  of  this  method 
is  that  one  can  define  assemblages  of  foraminifera  whose  relative  abundances  are  a  linear 
function  of  T,  and  that  the  function  defined  with  coretop  sediments  and  modern  SST  is 
applicable  in  the  past.  The  most  complete  fauna-SST  calibration  dataset  (Prell  et  al.,  1999) 
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Figure  3-4:  (a)  Flux  of  N.  pachyderma  (s.)  and  G.  bulloides  from  sediment  trap  deployments 
from  the  Gulf  of  Alaska  (PAPA,  September  1982  to  August  1986)  (Sautter  and  Thunell, 
1989),  the  far  northwest  Pacific  (50N,  Dec,  1997  to  Nov,  1998)(Kuroyanagi  et  al..  2002). 
and  the  Sea  of  Okhotsk  (SHOYU,  August  1990  to  August  1991)  (Alderman,  1989).  (b)-(d) 
color  contours  are  equilibrium  £18OcaiCjte  calculated  for  three  coring  sites,  from  Levitus  et  al. 
(1998)  data,  (e)  equilibrum  £18Ocaicite  calculated  from  WOCE  data  collected  in  summer. 
1993  (see  Figure  3-1  for  station  locations).  In  (b)-(e),  black  contour  lines  are  T  (°C). 
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includes  few  data  from  the  Pacific,  and  no  data  from  the  high-latitude  Pacific.  Therefore, 

I  do  not  try  to  convert  measured  changes  in  the  relative  abundance  of  N.  pachyderma  (s.) 
to  SST,  but  use  these  data  as  a  qualitative  measure  of  temperature. 

3.3  Results 

Lithology 

The  lithology  of  the  cores  consist  of  massive,  dysoxic,  and  laminated  sediment  (Figure  3-5). 
There  are  radiocarbon  dates  in  one  core  from  each  of  the  three  coring  regions  (Table  3.2). 
On  Bowers  Ridge,  cores  from  the  ridge  crest  typically  had  low  sediment  accumulation 
rates,  where  Termination  I  (TI)  and  the  Holocene  are  contained  in  the  top  10-30  cm  ( e.g . 
cores  25GGC  and  18GGC).  Core  17JPC  was  collected  on  the  west  flank  of  Rowers  Ridge, 
where  sediment  accumulation  rates  are  an  order  of  magnitude  higher  than  on  the  crest. 
I  tentatively  correlate  17JPC  with  the  other  Bowers  Ridge  cores  at  the  base  of  the  lami- 
nated/dysoxic  layers  in  each  core.  This  correlation  is  supported  by  N.  pachyderma  (s.)  5180 
in  25GGC,  where  I  can  identify  the  deglaciation  and  LGM.  In  18GGC,  the  correlation  is 
supported  by  features  in  the  magnetic  suceptibility  and  7- attenuation  data  (not  shown).  I 
use  the  bases  of  two  thick  black  ash  layers  found  in  each  Urnnak  Plateau  core  for  correlation, 
assuming  that  they  were  deposited  synchronously  and  instantaneously.  Using  these  corre¬ 
lations,  and  assuming  a  coretop  age  of  0  y  BP,  the  lithology  of  the  cores  is  plotted  against 
age  in  Figure  3-6.  Though  the  coretops  are  almost  certainly  older  than  Oy,  the  lithology 
shown  at  the  tops  of  the  cores  in  Figure  3-6  are  corroborated  by  the  companion  multicores 
at  these  sites.  The  ages  of  the  lithological  transitions  are  in  Table  3.4. 

Apparent  ventilation  ages 

I  estimated  apparent  ventilation  age  in  two  cores  during  the  Bplling  warm  phase.  In  17JPC, 
I  measured  14C  of  benthic  foraminifer  Nonionella  labradorica  at  its  maximum  abundance 
and  N.  pachyderma  (s.)  in  the  same  sample  (Figure  3-7  and  Table  3.2).  I  measured  the 
benthic-planktonic  pair  at  a  benthic  abundance  peak  in  17JPC  in  order  to  minimize  bias  in 
the  A14C  due  to  bioturbation  of  non-coeval  foraminifera  into  the  sample  (Broecker  et  ah, 
1999).  In  3JPC,  I  measured  Globohuhminia  affims  and  N.  pachyderma  (s.).  The  Bplling  is 
laminated  in  3JPC,  so  I  was  not  constrained  to  a  benthic  abundance  peak  (Figure  3-7  and 
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Figure  3-5:  Lithology  of  HLY02-02  cores.  Note  the  x-axes  are  different  for  each  core.  The 
radiocarbon  dates  measured  on  N.  pachyderma  (s.)  and  calibrated  to  calendar  years  with  the 
Marine04  dataset  with  Calib5.0.L  using  AR  =  400 y.  so  that  R  +  A R  =  800 y  (Table  3.2). 
Correlation  between  51JPC  and  the  undated  Umnak  Plateau  cores  is  at  the  base  of  2  black 
ash  layers  that  are  assumed  to  be  synchronously  deposited.  Tentative  correlation  is  made 
between  17JPC  and  the  undated  Bowers  Ridge  cores  at  the  base  of  the  laminated-dysoxic 
layer.  In  25GGC.  this  correlation  is  supported  by  N.  pachyderma  (s.)  <5180;  in  18GGC  it  is 
supported  by  magnetic  suceptibility  and  bulk  density  measurements. 
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Figure  3-6:  Top:  Lithology  of  HLY02-02  cores  versus  calibrated  calendar  age.  The  cores  are 
arranged  according  to  water  depth  of  coring  site.  Ash  layers  in  48JPC,  51 JPC,  and  55JPC 
are  assumed  to  have  been  deposited  instantaneously,  and  are  indicated  with  black  dashed 
lines.  See  Figure  3-5  for  a  legend  describing  symbols  and  shading.  Bottom:  [O2]  from  the 
WOCE  dataset  is  contoured,  with  coring  site  depths  indicated.  Cores  collected  from  the 
Bowers  Ridge  are  plotted  on  the  left,  cores  collected  from  the  continental  slope  are  plotted 
on  the  right.  The  cores  within  the  modern-day  OMZ  are  laminated  during  the  deglaciation. 
The  bathymetry  drawn  is  interpolation  of  the  depth  of  the  bottom  for  each  WOCE  station. 


45 


Figure  3-7:  I  measured  paired  benthic-planktonic  radiocarbon  dates  in  2  cores.  In  17JPC, 
I  chose  to  date  the  maximum  in  abundance  of  Nonionella  labradorica  (orange  star),  which 
is  1cm  above  the  base  of  the  dysoxic  layer  equivalent  to  the  BA.  In  3JPC,  I  measured 
Globobulimina  affinis  in  a  sample  from  the  laminated  interval  equivalent  to  the  B0lling.  In 
both  cores,  I  measured  N.  pachyderma  (s.)  in  the  same  samples. 

Table  3.2). 

To  develop  a  context  for  these  two  measurements,  I  compiled  published  and  unpub¬ 
lished  pairs  of  benthic  and  planktonic  foraminifer  14 C  dates  from  the  North  Pacific  region 
from  the  Last  Glacial  Maximum  (LGM,  23-19  ky  BP)  through  the  present.  In  Appendix 
A,  I  describe  in  detail  three  methods  for  quantifying  the  depletion  of  14 C  in  subsurface 
water  relative  to  surface  water:  (1)  the  simple  difference  between  benthic  14 C  age  and 
planktonic  14C  age,  A (2)  estimate  of  the  A14C  of  subsurface  water  relative  to  surface 
water,  A14Co— A14Csiy,  and  (3)  apparent  ventilation  age  of  subsurface  water,  tv,  taking  into 
account  changes  in  atmospheric  A14C. 


The  d180  and  faunal  record 

In  Figure  3-8  are  the  raw  isotope  data  and  N .  pachyderma  (s.)  relative  abundance.  The 
highest  N.  pachyderma  (s.)  5180  at  the  bottom  of  3JPC  (3.7  %o)  does  not  reach  the  max¬ 
imum  value  reached  in  N.  pachyderma  (s.)  <518 0  in  17JPC  (3.9 %o),  but  according  to  the 
age  model  in  17JPC,  at  20  ky  BP  the  two  cores  have  the  same  planktonic  d180  to  within 
measurement  uncertainties  (Figure  3-9  and  Figure  3-10).  The  N .  pachyderma  (s.)  relative 
abundance  in  3JPC  (70%)  is  lower  than  17JPC  (90%)  at  this  time  as  well,  but  since  there 
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Figure  3-8:  Raw  Jl80  data  from  N.  pachyderma  (s.)  (black)  and  Uvigerina  sp .  (green)  for 
HLY02-02-17JPC  (left)  and  HLY02-02-3JPC  (right).  Each  dot  is  an  analysis,  the  lines  go 
through  the  mean  value  in  each  sample  (n=3  for  each  sample).  Triangles  denote  location  of 
measured  radiocarbon  dates.  Blue  lines  are  relative  abundance  of  N ’  pachyderma  (s.),  with 
95%  confidence  limits  calculated  as  in  Fatela  and  Taborda  (2002).  The  age  model  assumes 
constant  sediment  accumulation  rate  between  radiocarbon  dates.  The  grey  vertical  bars  on 
the  left  panel  denote  dysoxic  lithological  units  in  17JPC,  on  the  right  panel,  they  denote 
laminated  layers  in  3JPC. 
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Figure  3-9:  17JPC  data  plotted  against  age.  In  the  left-hand  panel,  N.  pachyderma  (s.) 
5180  and  <513C  are  in  black.  Grey  contours  are  the  95%  confidence  limits.  Uvigerina  sp. 
<5180  and  <!)13 C  are  in  green.  Relative  abundance  of  N .  pachyderma  (s.)  is  in  blue,  where  the 
blue  countours  are  the  95%  confidence  limit.  In  the  right-hand  panel,  in  red  is  the  GISP2 
Jl80  of  ice  (Stuiver  et  ah.  1995).  In  orange  is  the  mean  change  in  <5l80  of  seawater  due  to 
the  rnelt  of  continental  ice  (Waelbroeck  et  aL,  2002).  On  this  panel,  the  mean-ocean  change 
in  J18 O  is  subtracted  from  the  17JPC  5180. 
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Figure  3-10:  3JPC  data  plotted  against  age.  In  the  left-hand  panel,  N.  pachyderma  (s.)  <5180 
and  (513C  are  in  black.  Grey  contours  are  the  95%  confidence  limits.  Relative  abundance  of 
N.  pachyderma  (s.)  is  in  blue,  where  the  blue  countours  are  the  95%  confidence  limit.  In 
the  right-hand  panel,  in  red  is  the  GISP2  <5i80  of  ice  (Stuiver  et  al.,  1995).  In  orange  is 
the  mean  change  in  <5180  of  seawater  due  to  the  melt  of  continental  ice  (Waelbroeck  et  al., 
2002).  On  this  panel,  the  mean-ocean  change  in  £180  is  subtracted  from  the  3JPC  <5180. 
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Figure  3-11:  N.  pachyderma  (s.)  5lsO  from  17JPC  and  its  two  companion  cores  15GGC  and 
16MC-C.  The  1.5  %o  increase  in  the  top  30  cm  of  17JPC  is  not  seen  in  its  two  companion 
cores.  I  conclude  that  it  is  not  a  climatic  signal. 


are  few  faunal  data  from  the  bottom  3  meters  of  3JPC,  this  may  not  be  a  robust  observation. 

In  the  TV.  pachyderma  (s.)  61S O  record  in  17JPC,  5l80  begins  to  decrease  35  cm  below 
the  first  dysoxic  layer.  At  the  lithological  transition  into  the  dysoxic  layer  there  is  a  rapid 
decrease  in  TV.  pachyderma  (s.)  relative  abundance.  In  3JPC,  where  the  sediment  accu¬ 
mulation  rate  is  an  order  of  magnitude  higher,  two  decreases  in  TV.  pachyderma  (s.)  <5180 
occurred  very  rapidly  at  the  boundary  between  massive  and  laminated  sediment,  which  cor¬ 
responds  to  the  beginning  of  the  Bplling  and  the  end  of  the  Younger  Dryas.  The  amplitude 
of  these  8iS O  events  is  0.7-0.8%c.  These  changes  in  <5180  are  accompanied  by  decreases  in 
the  relative  abundance  of  TV.  pachyderma  (s.). 

During  the  Holocene,  the  TV.  pachyderma  (s.)  (5180  in  3JPC  decreases  continuously  to 
1.8  %o  at  the  top  of  the  core  (Figure  3-8).  In  17JPC,  TV.  pachyderma  (s.)  <5180  reaches  a 
minimum  of  1.3  %o  at  30  cm,  then  increases  sharply  by  1.8  %o  to  3.1  %o  at  the  top  of  the 
core.  This  core-top  value  is  not  seen  in  the  companion  gravity  core  (15GGC)  and  multicore 
(16MC-C)  for  this  site  (Figure  3-11),  calling  into  question  whether  it  is  a  real  climatic  signal. 
In  17JPC,  the  Uvigerina  sp.  <5180  curve  is  dashed  above  128  cm  because  the  5  analyses  above 
this  point  were  conducted  on  the  only  5  Uvigerina  sp.  found  in  the  core  0-128  cm,  and  I 
am  not  confident  they  are  in  place,  and  not  bioturbated  upward  from  deeper  in  the  core 
where  they  are  more  abundant  (see  Discussion). 

3.4  Discussion 

Chronology 

The  age  model  in  each  core  assumes  a  constant  sediment  accumulation  rate  between  14  C 
dates.  Though  I  do  not  expect  the  sediment  accumulation  rate  in  these  cores  to  not  be 
smoothly  varying,  this  age  model  is  very  simple.  The  greatest  uncertainty  in  the  age  models 
is  the  magnitude  of  the  reservoir  correction,  R  +  A R.  The  A R  is  the  local  anomaly  from 
the  global  mean  surface  water  apparent  age  (Bard,  1988)  and  is  very  poorly  constrained  in 
the  modern  North  Pacific.  Physical  processes  which  affect  the  magnitude  of  AT?  include 
circulation,  upwelling,  and  sea-ice  formation  and  melting,  which  affects  air-sea  gas  exchange 
and  upper-ocean  stratification. 

In  order  to  convert  measured  radiocarbon  ages  to  calendar  ages,  I  must  choose  a  rea¬ 
sonable  value  for  A R}  and  assume  that  it  has  remained  constant  through  the  deglaciation. 
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The  nearest  locations  where  there  are  modern  estimates  of  A R  are  the  Kuril  Islands  (320  y, 
Kuzmin  et  al.  (2002))  and  Kamchatka  Strait  (240 y,  Robinson  and  Thompson  (1981)).  I 
chose  A R  =  400  y  because  if  it  were  any  smaller,  there  are  events  in  3JPC.  I  believe  that 
these  layers  are  correlated  with  GISP2  would  precede  the  events  in  the  ice  core. 

The  N.  pachyderma  (s.)  UC  age  measured  at  215  cm  is  200  y  older  than  the  age  at 
216  cm,  identifying  an  age  reversal,  though  these  two  dates  are  indistinguishable  within  2 a 
uncertainty.  I  do  not  use  the  age  at  215  cm  in  constructing  the  age  model,  but  do  use  it 
in  calculating  the  apparent  ventilation  age  of  this  sample.  I  discuss  the  reliability  of  the 
measurements  in  17JPC  later  in  the  next  section. 

Lithology 

The  oxygen  minimum  zone  (OMZ)  in  the  modern  Bering  Sea  is  most  intense  at  about 
1000  m  (Figure  3-6b),  where  lowest  [O2]  measured  on  the  WOCE  profile  is  15/imol/kg.  No 
laminated  coretop  sediments  have  been  found  in  the  Bering  Sea.  During  the  deglaciation, 
it  appears  that  the  OMZ  intensified,  so  that  cores  from  within  the  OMZ  (18GGC,  55JPC, 
3JPC,  51JPC,  and  48JPC)  have  intervals  of  laminated  sediment.  Core  17JPC  (from  2209  m 
water  depth),  is  deep  enough  to  be  below  the  OMZ.  There  are  contemporaneous  dysoxic 
layers  with  sharp  lower  contacts  during  those  two  time  periods  in  17JPC  and  dysoxic  sed¬ 
iments  elsewhere  in  the  Pacific  (van  Geen  et  al.,  2003;  Zheng  et  al.,  2000)  which  I  believe 
are  analogous  to  the  laminated  layers  in  the  shallower  cores  because  they  occur  during  BA 
and  the  early  Holocene. 

If  the  cores  with  lowest  sediment-accumulation  rate,  25GGC  and  18GGC,  also  accumu¬ 
lated  laminated  sediments  during  these  two  time  periods,  the  evidence  of  laminae  would 
be  destroyed  with  eventual  increase  in  [O2]  during  the  Holocene,  as  benthic  macrofauna 
homogenized  the  uppermost  sediments  (the  activity  of  benthic  macrofauna  is  typically  de¬ 
scribed  mathematically  as  Fickian  diffusion  (Berger  and  Heath,  1968),  with  a  characteristic 
length  scale).  For  the  same  reason,  in  the  other  cores,  the  top  boundary  of  a  laminated 
interval  with  the  overlying  massive  interval  probably  appears  stratigraphically  lower  than 
when  it  was  originally  deposited  by  about  the  thickness  of  the  length  scale  of  diffusive  bio- 
turbation,  which  is  typically  on  the  order  of  10 cm  (Boudreau,  1994).  For  simplicity,  I  do 
not  attempt  to  correct  the  apparent  age  of  the  top  of  an  interval  of  laminated  sediment. 
So  the  interpolated  ages  of  these  lithological  transitions  in  Table  3.4  may  be  biased  toward 
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slightly  older  ages. 

Bioturbation  can  be  heterogeneous,  or  “advective,”  where  chunks  of  sediment  are  trans¬ 
ported  up  or  down  in  the  sediment  column  relative  to  the  layer  in  which  it  originally 
accumulated  (Frangois  et  al.,  2001).  This  is  suggested  in  the  deglacial  sediments  of  17JPC, 
where  the  two  dysoxic  layers  have  such  a  different  texture  and  color  than  the  intervening 
sediments,  that  centimeter-scale  mottling  is  visible.  It  is  likely  that  such  processes  occur 
in  the  rest  of  this  sediment  core,  even  in  the  absence  of  color  changes  in  the  sediment  that 
allow  the  mottling  to  be  easily  seen. 

Given  the  sedimentation  rate  in  this  core,  the  thickness  of  sediment  that  would  accumu¬ 
late  during  a  millennial-scale  climate  event  is  on  the  order  of  10  cm,  which  is  on  the  same 
spatial  scale  as  the  mottling  in  17JPC.  This  means  that  a  sample  analyzed  for  5180  could 
contain  non-coeval  foraminifera,  and  if  the  individuals  span  a  transition  in  5180,  the  tran¬ 
sition  in  the  measured  profile  could  be  offset  significantly  from  where  it  would  be  observed 
if  the  core  were  not  bioturbated.  A  radiocarbon  measurement  of  a  sample  of  foraminifera 
composed  of  individuals  which  are  not  coeval  could  be  significantly  biased  from  the  age  of 
the  depth  horizon  from  which  the  sample  was  collected.  The  sedimentation  rate  of  3JPC  is 
an  order  of  magnitude  higher  than  17JPC,  and  so  the  isotopic  and  14 C  record  in  this  core 
is  more  reliable. 

The  planktonic  and  benthic  5180  measured  in  the  uppermost  20  cm  of  17JPC  increases 
by  1.5  %o.  I  do  not  see  such  an  isotopic  feature  in  any  other  Bering  Sea  cores.  I  expect 
the  companion  gravity  core  at  this  site,  15GGC,  to  stratigraphically  overlap  with  the  top 
of  17JPC  and  contain  more  of  the  sediment  column  closest  to  the  sediment- water  inter¬ 
face.  Neither  the  gravity  core  nor  the  companion  multicore  show  such  an  isotopic  feature 
(Figure  3-11).  Therefore  I  conclude  that  the  planktonic  record  in  17JPC  was  somehow  dis¬ 
turbed  during  the  coring  procedure,  or  from  bioturbation  of  high-518 O  individuals  up  from 
deeper  within  the  core.  Both  benthic  and  planktonic  foraminifera  gradually  become  very 
rare  toward  the  top  of  the  core,  so  that  even  a  small  amount  of  contamination  of  sediment 
containing  foraminifera  with  different  5180  could  have  a  large  impact  on  measured  5180 
values  at  the  top  of  the  core  ( N .  pachyderma  (s.)  5180  measurements  are  made  on  a  sampe 
of  8  individuals).  In  the  case  of  Uvigerina  sp.  (which  are  analyzed  one  at  a  time),  the  sole 
five  individuals  found  and  analyzed  in  top  128  cm  of  17JPC  could  have  been  transported 
from  below. 
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GISP2 

(y  BP) 

3JPC 

(y  bp) 

17JPC 
(y  BP) 

A  ice  -  3JPC) 
(y) 

A(ice  -  17JPC) 

(y) 

Pre-Boreal 

11,650±230 

ll,540±i8o 

12,240±go 

no±38 

-590±250 

Younger  Dryas 

End-Aller0d 

12,890±260 

12,880±i2o 

13,600±100 

io±i§ 

-710±280 

iAcp 

13,070±260 

13,070±1°° 

0±29O 

Aller0d 

13,250±260 

13,440±i(j0 

-190±280 

14,010±280 

14,000±i5q 

10±320 

Older  Dryas 

14,090±280 

14,110±i!o 

-20±330 

Bplling 

14,670±290 

14,690±2io 

14,830±j5o 

-20±ig 

-160±330 

Oldest  Dryas 

mean  offset 

-10±110 

-530±160 

Table  3.4:  Age  in  y  BP  of  the  mid-point  of  transitions  in  GISP2  <5180;ce  (Stuiver  et  ah,  1995), 
and  the  age  of  lithological  transitions  between  massive  and  laminated  sediment  in  3JPC, 
and  massive  and  dysoxic  sediment  in  17JPC.  In  the  age  models,  I  calibrate  radiocarbon 
dates  with  A R  —  400  y,  so  that  the  total  reservoir  correction,  R  =  800  y. 

The  ages  of  transitions  between  massive  and  laminated  sediment  in  3JPC  and  17JPC 
are  compared  to  the  ages  of  transitions  between  climate  events  in  the  Greenland  ice  cores 
in  Table  3.4  and  Figure  3-12.  In  the  ice  cores,  the  age  of  a  transition  is  considered  the 
mid-point  of  the  change  in  <5180jce  (Stuiver  et  ah,  1995).  The  offsets  between  the  climate 
events  in  the  Greenland  ice  cores  and  3JPC  range  from  -190  to  110  years,  which  is  within  the 
uncertainty  of  the  14C  age  model.  The  offests  between  GISP2  and  17JPC  are  significantly 
larger,  ranging  from  -710  to  -160  y,  where  the  transitions  in  17JPC  are  older  than  Greenland. 

At  the  start  of  the  deglaciation,  the  decrease  in  planktonic  Sl80  in  17JPC  is  gradual 
and  begins  earlier  than  in  3JPC  (Figure  3-9  and  3-10).  This  is  similar  to  what  is  seen 
in  a  compilation  of  SST  proxy  datasets  from  the  Pacific  by  Kiefer  and  Kienast  (2005). 
The  greater  age  of  the  transitions  in  the  lithology  of  17JPC  as  well  as  the  decrease  in 
N.  pachy derma  (s.)  <P80  that  occurs  below  the  lithological  transition  to  dysoxic  sediment 
could  both  be  the  result  of  the  presence  of  N .  pachyderma  (s.)  in  sediments  older  than  when 
they  were  deposited.  If  this  is  the  case,  the  differences  between  the  records  in  the  17JPC 
and  3JPC  cores  may  not  be  climatic,  but  rather,  an  artifact  of  post-depositional  processes. 

If  we  assume  that  the  lithological  changes  in  3JPC  and  the  5180  changes  in  GISP2  are 
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age  (ky)  in  GISP2 


Figure  3-12:  Age  of  events  in  GISP2  plotted  against  age  of  lithological  transitions  in  3JPC 
(black)  and  17JPC  (green).  The  vertical  error  bars  are  the  la  uncertainty  in  the  age  model. 
The  horizontal  error  bars  are  the  estimated  uncertainty  in  the  absolute  age  of  the  layer- 
counted  age  model  for  GISP2  is  2%,  but  the  relative  uncertainty  is  much  smaller  (Meese 
et  ah,  1997). 


synchronous,  this  implies  that  A R  did  not  change  significantly  during  the  deglaciation.  I 
chose  A R  —  400  y  because  if  it  were  smaller,  the  lithological  changes  in  3JPC  would  preceed 
climate  transitions  in  GISP2.  Adjustment  of  A R  would  change  the  absolute  timing  of  the 
events,  but  not  the  relative  timing,  and  the  magnitude  of  the  mismatch  between  the  timing 
of  the  events  is  well  within  the  uncertainty  of  the  age  model.  Therefore,  if  A R  did  not 
change  much  during  the  deglaciation,  the  correlation  between  GISP2  and  the  lithology  of 
3JPC  is  not  sensitive  to  the  absolute  value  of  A R. 

A  relatively  constant  A R  is  consistent  with  no  change  in  sea  ice  extent  or  no  change  in 
stratification  of  the  upper  water  column,  processes  which  can  affect  A R.  This  is  a  simpler, 
and  in  my  opinion,  more  likely  scenario  than  if  A R  was  not  constant  during  the  deglaciation, 
because  this  would  require  a  separate  mechanism  which  would  make  the  lithological  events 
appear  to  have  the  same  relative  timing. 

Apparent  ventilation  ages 

To  test  whether  reduced  ventilation  of  intermediate  water  contributed  to  the  low  [O2]  con¬ 
ditions  revealed  by  the  lithology  of  these  cores,  I  compiled  all  the  published  pairs  of  benthic 
and  planktonic  14C  measurements  from  the  North  Pacific  and  added  unpublished  data  of 
L.  Keigwin,  and  the  two  new  data  points  from  this  study.  Analysis  of  the  compilation 
shows  that  there  was  no  consistent  change  in  the  14C  of  subsurface  water  shallower  than 
2700  m  between  the  BA  and  the  Younger  Dryas  (YD,  12.9-11.7  ky  BP)  or  Last  Glacial 
Maximum  (LGM,  23-19  ky  BP)  (Figure  3-13).  There  are  some  limitations  to  the  dataset 
used,  including  the  small  number  of  data  points,  especially  during  the  YD.  A  particular 
challenge  in  measuring  14C  on  foraminifer  samples  is  choosing  an  interval  of  sediment  in 
which  the  data  will  be  the  least  biased  by  bioturbation.  This  potential  bias,  combined  with 
the  uncertainty  in  the  transformation  of  a  14 C  age  measurement  to  calendar  years,  has  the 
potential  to  prevent  us  from  being  able  to  know  how  the  changes  in  the  benthic-planktonic 
14C  data  from  different  cores  are  related  to  one  another  in  time.  This  is  particularly  true 
for  the  deglaciation  since  the  BA  and  YD  are  brief  enough  that  they  are  not  always  pos¬ 
sible  to  resolve  well  in  sediment  cores.  If  I  had  observed  a  significant  decrease  in  14C  of 
intermediate-depth  water,  it  may  not  necessarily  mean  that  there  was  a  decrease  in  the 
rate  of  formation  of  that  water  mass.  Reduced  rate  of  formation  of  NPIW  would  result  in 
lower  oxygen  and  14 C  content  of  intermediate-depth  water  only  if  this  were  the  only  process 
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Figure  3-13:  This  is  a  summary  of  the  analysis  of  a  compilation  of  North  Pacific  pairs  of 
benthic  and  planktonic  measurements  of  14C  conducted  in  Appendix  1.  The  three  columns, 
from  left  to  right,  are  time  slices  from  the  Younger  Dryas  (YD,  panels  a,  d,  and  g),  the 
R0lling-Aller0d  (BA,  panels  b,  e,  and  h),  and  the  Last  Glacial  Maximum  (LGM,  panels 
c,  f,  and  i).  The  three  rows,  from  top  to  bottom,  are  A^_p,  the  difference  between  the 
14  C  age  of  benthic  and  planktonic  foraminifera  (panels  a-c);  tv,  the  apparent  ventilation 
age  of  subsurface  water,  taking  into  account  change  in  atmospheric  A14C  (panels  d-f); 
and  A14Co~ A14Cslt;,  the  difference  between  A14C  of  subsurface  water  and  surface  water 
(panels  g-i).  For  comparison  to  the  modern  ocean,  on  each  panel,  the  solid  grey  profile  is 
the  average  North  Pacific  A14C  from  the  WOCE  dataset  (Key  et  ah,  2002).  In  panels  g-i, 
the  dashed  grey  profile  is  the  modern  A14C  profile,  corrected  for  the  higher  amount  of 14 C  in 
the  ocean-atmosphere  system  for  each  time  slice.  See  Figure  A-4  for  the  calculation  of  mean 
A14C  for  each  time  slice.  Each  open  circle  is  a  data  point  which  falls  within  the  time  slice. 
The  closed  circles  are  data  points  which  fall  within  the  time  slice  with  2 a  certainty.  The 
error  bars  are  la.  See  the  Appendix  for  details  on  the  calculation  of  these  three  quantities, 
and  Table  A.l  for  data  and  references.  The  vertical  profiles  of  A 6_p,  tv,  and  A 14 Co  —  A14CSf(. 
are  not  significantly  different  from  the  modern  profile  in  any  of  the  3  time  slices. 
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that  was  changing.  For  example,  a  change  in  the  14C  of  end-member  components  of  the 
water  mass  at  this  depth  could  affect  14C  without  a  change  in  circulation  (Wunsch,  2003); 
a  change  in  the  oxygen  content  of  the  Okhotsk  Sea  water  mass  which  forms  NPIW  could 
affect  the  oxygen  conent  of  intermediate  water,  without  a  change  in  formation  rate  (Crusius 
et  al.,  2004). 

The  two  new  data  points  from  this  study  were  measured  in  samples  where  known  bias 
from  bioturbation  was  minimized.  In  17JPC,  the  sample  was  from  a  maximum  of  absolute 
abundance  of  the  benthic  foraminifer.  In  3JPC,  the  sample  was  from  a  laminated  interval 
of  sediment.  With  eventual  expansion  of  the  North  Pacific  14C  dataset,  this  question  could 
be  revisited  in  the  future,  using  only  data  where  efforts  have  been  made  to  minimize  bias 
from  post-depositional  processes. 

Deglacial  SST  and  salinity 

When  plotted  against  age,  the  shape  of  the  N.  pachyderma  (s.)  <5180  record  in  3JPC  is 
similar  to  the  <5180ice  record  from  GISP2  (Figure  3-10),  where  there  is  a  gradual  increase  in 
N.  pachyderma  (s.)  <5180  in  3JPC  (to  relatively  colder  temperatures)  in  3JPC  and  a  gradual 
decrease  in  Sxs Olce  in  GISP2  (to  relatively  colder  temperatures). 

There  are  large  changes  in  the  fraction  N.  pachyderma  (s.)  at  the  same  time  as  the 
changes  in  planktonic  5180  (Figure  3-8).  At  station  PAPA  in  the  Gulf  of  Alaska,  Sautter 
and  Thunell  (1989)  estimate  that  only  about  25%  of  the  foraminifera  from  the  upper  water 
column  are  preserved  in  the  sediment,  the  rest  dissolve  or  are  broken  up  in  the  descent  to 
seafloor  and  subsequent  burial  and  mixing.  A  change  in  the  relative  abundance  of  a  par¬ 
ticular  species  in  the  sediment  record  may  be  the  consequence  of  a  change  in  the  amount 
of  dissolution  of  the  CaC03  fraction  of  the  sediment  that  has  occurred,  if  the  species  is 
more  or  less  resistant  to  dissolution  (Berger  and  Piper,  1972).  Across  both  these  transi¬ 
tions,  the  sediment  becomes  less  oxic,  and  the  increase  in  pC02  in  porewater  from  organic 
matter  oxidation  would  undersaturate  the  porewater  with  respect  to  calcite.  If  N.  pachy¬ 
derma  (s.)  is  more  susceptible  to  dissolution  than  the  other  dominant  species  in  the  plank¬ 
tonic  foraminiferal  assemblage,  G.  bulloides ,  then  this  could  account  for  an  exaggeration 
in  the  change  in  its  relative  abundance.  However,  the  mean  mass  of  a  N.  pachyderma  (s.) 
test  from  the  212-250  ^m  size  fraction  does  not  change  systematically  through  the  deglacial 
transition  (data  not  shown),  which  suggests  the  relative  abundance  is  probably  not  affected 
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by  a  differential  preservation  across  these  events. 

The  rapid  decreases  in  N.  pachyderma  (s.)  6IS O  of  -0.7  and  -0.8  %o  at  the  beginning  of 
the  Bplling  and  the  end  of  the  YD  in  3JPC  are  probably  due  to  a  combination  of  warming 
(from  the  faunal  data  in  this  study)  and  freshening  (Cook  et  al.,  2005),  In  each  case, 
the  low  relative  abundance  of  N.  pachyderma  (s.)  is  transient,  but  the  5180  stays  low, 
suggesting  that  though  the  warming  itself  is  transient,  there  is  a  step  decrease  in  salinity 
at  this  transition. 

If  an  increase  in  precipitation  contributed  to  the  freshening,  this  would  be  consistent 
with  a  scenario  in  which  there  was  a  strong  East  Asian  Monsoon  during  warm  climate 
events.  Speleothem  <5180  records  dated  with  uranium-series  isotopes  from  southeast  China 
suggest  the  East  Asian  monsoon  (EAM)  was  strong  during  the  BA,  weak  during  the  YD, 
and  strong  in  the  early  Holocene  (Wang  et  al.,  2001).  Flux  of  freshwater  to  the  North  Pacific 
subpolar  gyre  from  the  EAM  is  thought  to  contribute  to  the  relatively  low  salinity  of  the 
subpolar  North  Pacific  (Emile-Geay  et  al.,  2003;  Warren,  1983).  In  addition,  the  geometry 
of  the  Pacific  basin  results  in  relatively  little  exchange  of  water  between  the  subtropical  and 
subpolar  gyres  (Emile-Geay  et  al.,  2003).  A  consequence  of  the  fresh  surface  water  of  the 
subpolar  North  Pacific  is  that  it  is  not  able  to  support  convective  overturning  such  as  the 
AMOC. 

We  thus  speculate  that  the  strength  of  the  EAM  may  influence  the  amount  of  sea-ice 
formation  and  the  density  of  the  brine  produced  in  the  Okhotsk  and  Bering  Seas,  and 
ultimately,  the  rate  of  formation  of  intermediate  water.  A  greater  flux  of  freshwater  to  the 
subpolar  North  Pacific  (as  precipitation  associated  with  the  EAM)  could  have  decreased 
the  salinity  of  surface  water.  If  fresh  Okhotsk  Sea  surface  water  inhibits  the  formation  of 
NPIW,  this  could  result  in  decrease  in  ventilation  and  oxygen  content  of  intermediate  water, 
resulting  in  laminated  and  dysoxic  sediments  throughout  the  North  Pacific  basin.  However, 
another  effect  of  the  freshening  of  surface  water  could  be  an  increase  in  production  of  winter 
sea  ice,  which  would  have  the  opposite  effect,  an  increase  in  production  of  brine,  and  NPIW. 

Another  source  of  freshwater  to  the  Bering  Sea  could  be  glacial  meltwater  from  Siberia 
and  Alaska.  Sancetta  et  al.  (1985)  suggested  that  the  low  <5180  event  in  bulk  opal  (— 9%o) 
during  the  deglaciation  from  a  core  in  the  southeast  Bering  Sea  could  be  due  to  a  flux 
of  freshwater  from  the  Yukon.  It  is  generally  believed  that  Siberia  and  Alaska  were  not 
heavily  glaciated  during  the  LGM,  based  on  cosmogenic  dating  of  terrestrial  glacial  deposits 
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(Brigham-Grette  et  al.,  2003).  However,  Hill  et  al.  (2006)  have  found  deep  incised  channels 
on  the  Beringian  shelf  which  were  completely  excavated  then  filled  with  sediment  since  the 
last  deglaciation,  implying  a  large  flux  of  meltwater  from  Alaska.  This  would  mean  that 
Alaska  was  more  heavily  glaciated  than  previously  thought. 

Mikolajewicz  et  al.  (1997)  perturbed  a  coupled  ocean-atmosphere  model  with  a  large  flux 
of  freshwater  to  the  North  Atlantic  in  order  to  study  the  far  field  effects  of  the  consequent 
YD-like  climate  event.  They  find  that  in  their  model,  a  change  in  atmospheric  circulation 
cools  surface  waters  in  the  North  Pacific,  increasing  ventilation  of  thermocline-depth  water. 
We  thus  conclude  that  both  cooling  of  North  Pacific  subpolar  surface  water  from  cold 
Siberian  winds,  and  freshwater  flux  to  the  subpolar  gyre  from  the  EAM  are  atmospheric 
mechanisms  which  could  explain  the  correlation  between  GISP2  and  the  lithology  at  3JPC. 

3.5  Conclusions 

Laminated  and  dysoxic  sediments  accumulated  around  the  rim  of  the  North  Pacific  at 
OMZ  depths  during  the  Bplling-Allerpd  and  early  Holocene.  Analysis  of  a  compilation 
of  published  benthic-planktonic  14  C  measurements  from  the  North  Pacific  does  not  clarify 
whether  the  lithological  changes  were  due  to  a  change  in  the  rate  of  formation  of  intermediate 
water.  There  is  a  strong  correlation  between  the  timing  of  changes  between  laminated  and 
massive  sediments  in  the  Bering  Sea  and  the  <5180  of  ice  at  GISP2,  which  is  robust,  even 
given  the  uncertainty  in  the  local  reservoir  correction  (A R)  at  this  site.  The  correlation 
suggests  that  there  are  consistent  and  rapid  variations  of  both  air  temperature  at  the  summit 
of  Greenland  and  the  oxygen  content  of  subsurface  water  in  the  North  Pacific  at  centennial 
and  millenial  time  scales  during  the  deglaciation,  even  though  the  climate  system  was 
undergoing  a  major  transition. 
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Chapter  4 


Marine  Isotope  Stage  3  in  the 
Bering  Sea 


Abstract 

In  HLY02-02-17JPC,  a  core  from  the  west  flank  of  the  Bowers  Ridge  in  the  south-central 
Bering  Sea,  I  constructed  records  of  benthic  and  planktonic  £180  and  <513 C  and  relative 
abundance  of  N.  pachyderma  (s.),  a  proxy  for  sea  surface  temperature.  There  is  statistically 
significant  but  low-amplitude  variability  in  the  faunal  record.  Episodes  of  apparent  warm 
surface  water  do  not  correlate  with  any  of  the  most  prominent  features  of  the  benthic 
or  planktonic  isotopic  records.  Though  interpretation  of  this  dataset  is  limited  by  the 
preliminary  age  model,  either  there  is  not  a  coherent  relationship  between  the  isotopic 
and  faunal  records  with  Heinrich  or  Daansgard-Oeschger  climate  events  of  the  last  glacial 
period,  or  the  effects  of  these  events  in  the  Bering  Sea  region  are  too  subtle  to  be  detected 
in  these  particular  proxy  records.  These  findings  do  not  corroborate  previous  investigations 
which  conclude  that  there  was  significant  variability  in  subsurface  Sl3C  and  sea  surface 
temperature  in  the  North  Pacific  coinciding  with  Dansgaard-Oeschger  and  Heinrich  events. 


4.1  Introduction 

During  Marine  Isotope  Stage  3  (MIS3,  ca.  50-28  ky  BP),  there  were  dramatic  millennial- 
scale  climate  events,  which  appear  to  have  been  global  in  impact  (see  Voelker  and  workshop 
participants  (2002)  for  a  review)  ,  though  the  most  extensive  geological  records  documenting 
the  events  are  from  the  North  Atlantic  region.  Warm  Dansgaard-Oeschger  events  (or  inter- 
stadials)  are  associated  with  high  air  temperature  in  Greenland  and  elevated  atmospheric 
CH4  (Brook  et  al.,  1996).  Kennett  et  al.  (2003)  proposed  that  warmer  upper-intermediate 
water  during  interstadials  destabilizes  methane  hydrates  from  continental  margin  sediments. 
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This  releases  CH4  gas  into  the  water  column,  and  if  the  flux  is  large  enough,  not  all  of  the 
CH4  is  oxidixed  in  the  water  column,  and  some  escapes  to  the  atmosphere.  This  CH4, 
a  strong  greenhouse  gas,  could  provide  a  positive  feedback  of  interstadial  warming.  Ken- 
nett  et  al.  (2003)  axgue  that  wetlands  could  not  respond  rapidly  enough  to  account  for  the 
observed  fluctuations  in  atmospheric  CH4.  There  is  sedimentary  evidence  of  the  episodic 
presence  of  methane  in  marine  sediments  during  MIS3  (see  review  in  the  Introduction  of 
Chapter  5). 

Our  knowledge  of  the  glacial  history  of  the  Pacific  is  sparse.  Arguably  the  best  sedi¬ 
mentary  climate  record  from  anywhere  in  the  Pacific  comes  from  Santa  Barbara  Basin,  an 
anoxic  basin  on  the  California  margin  (sill  depth  =  475  m),  where  there  is  a  good  corre¬ 
lation  between  interstadial  events  and  lithology;  laminated  sediments  occur  during  every 
interstadial  (Behl  and  Kennett,  1996).  However,  at  this  coastal  site,  it  is  uncertain  whether 
the  laminated  sediment  is  due  to  reduced  ventilation  of  upper  intermediate  water  in  the 
Pacific  or  the  effects  of  a  local  change  in  export  production.  Elsewhere  in  the  North  Pacific, 
reconstructions  of  glacial  oceanography  have  been  hampered  by  the  typically  low  sediment 
accumulation  rates  endemic  to  this  ocean  basin.  From  a  site  on  the  Emperor  seamounts, 
Kiefer  et  al.  (2001)  concluded  that  North  Pacific  and  North  Atlantic  sea  surface  tempera¬ 
tures  (SSTs)  are  out  of  phase,  and  that  warm  Pacific  SSTs  are  caused  by  reduced  upwelling 
of  cold  water.  However,  uncertainties  in  the  age  model  for  this  core  call  into  question 
whether  the  planktonic  SST  record  is  in  or  out  phase  with  GISP2.  In  a  record  from  a 
2712m-deep  core  in  the  northeast  Pacific,  Lund  and  Mix  (1998),  find  a  correlation  between 
high  benthic  <513  C  and  Heinrich  events.  Heinrich  events  were  the  coldest  millennial-scale 
climate  events  of  the  last  glacial  period.  In  the  North  Atlantic,  nutrient  proxies  measured 
from  during  Heinrich  events  are  consistent  with  a  reduction  of  the  strength  of  the  Atlantic 
meridional  overturning  circulation  (AMOC)  (Vidal  et  al.,  1997;  Keigwin  and  Lehman,  1994; 
Curry  and  Oppo,  1997;  Keigwin  and  Boyle,  1999).  Lund  and  Mix  (1998)  interpret  their 
data  as  indicating  greater  ventilation  of  deep  Pacific  water  when  the  AMOC  was  thought 
to  be  weak.  This  interpretation  is  similar  to  better  ventilated  intermediate  depths  inferred 
from  LGM  <513C  on  the  Emperor  Seamounts  (Keigwin,  1998). 

A  new  sediment  core,  collected  from  2200  m  water  depth,  from  Bowers  Ridge,  in  the 
south-central  Bering  Sea,  has  a  sediment  accumulation  rate  during  MIS3  that  is  comparable 
or  higher  than  those  of  existing  records  of  MIS3  variability,  and  could  help  resolve  some  of 
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the  ambiguity  of  published  records  from  relatively  low-sediment  accumulation  rate  cores. 
With  the  Bowers  Ridge  core,  I  will  test  the  interpretations  in  previous  studies  that  there 
was  significant  millennial-scale  variability  in  the  513C  of  intermediate/deep  water,  and  that 
sea  surface  temperature  in  the  far  North  Pacific  was  high  during  DO  events. 


4.2  Methods 

Sample  preparation 

HLY02-02-17JPC  was  collected  from  a  site  at  2209  m  water  depth  on  the  western  flank  of 
the  Bowers  Ridge  (Figure  3-1).  Sediment  samples  were  dried  at  50°C,  then  washed  with  tap 
water  through  a  63  pm  sieve.  The  sample  was  dried  again  at  50°C  and  the  dry  sample  was 
sieved  so  that  planktonic  and  benthic  foraminifera  could  be  picked  from  larger  size  fractions. 
N.  pachyderma  (s.)  and  Uvigerina  sp.  were  picked  from  the  212-250  pm  and  >250  pm  size 
fractions,  respectively.  5lsO  and  <513C  measurements  of  8-10  N.  pachyderma  (s.)  or  a  single 
Uvigerina  sp.  were  made  on  a  PRISM  mass  spectrometer  and  are  reported  relative  to 
PDB  via  NBS-19  (Coplen,  1996).  At  least  300  individual  planktonic  fauna  >150  pm  were 
counted,  and  95%  confidence  limits  on  the  relative  abundance  of  species  were  calculated  as 
in  Fatela  and  Taborda  (2002).  Relative  abundance  of  N.  pachyderma  (s.)  is  a  proxy  for  sea 
surface  temperature,  where  higher  relative  abundance  corresponds  to  lower  temperatures 
(see  Chapter  3). 


Age  model 

Samples  of  N.  pachyderma  (s.)  were  analyzed  for  14 C  at  the  National  Ocean  Sciences  Ac¬ 
celerator  Mass  Spectrometer  Facility  at  WHOI.  I  calibrated  measured  radiocarbon  dates 
with  Calib5.0.1,  using  a  reasonable  value  for  the  reservoir  correction,  AR=400y  (R  +  AR  = 
800  y).  The  surface  ocean  reservoir  correction  is  poorly  constrained  in  the  modern  North 
Pacific,  see  Cook  et  al.  (2005)  for  discussion.  The  age  model  is  simply  constructed  as  linear 
interpolation  of  median  calibrated  ages.  To  calibrate  radiocarbon  dates  beyond  25,000  y  to 
calendar  years,  I  use  a  record  of  14  C  age  of  Globigerina  bulloides  from  the  Cariaco  Basin 
on  the  GRIP  ss09sea  age  model.  The  average  sample  spacing  in  17JPC  is  230  y. 
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depth  (cm) 

14C  y  BP±ct 

cal  y  BP  t  cr 

1  2cr 

155.5- 156.5 

185.5- 186.5 

10,000±60 

12,550±65 

in  icn  10,400 
10»460  10,540 

13,600 

— iu,^yu 

10,570 

13,430 

13,750 

215.5- 216.5 

383.5- 384.5 

13,400±60 

25,200±130 

14,830  J}g 
28,500 

14,440 

15,130 

551.5-552.5 

41,500±370 

44,530 

668.5-669.5 

49,100±640 

52,310 

Table  4.1:  Uncorrected  AMS  14C  ages  and  corrected,  calibrated  ages  of  TV.  pachyderma  (s.) 
in  HLY02-02-17JPC.  For  the  three  ages  younger  than  2514C  ky  BP,  calibrations  are  per¬ 
formed  with  Calib5.0html  (Reimer  et  ah,  2004)  using  the  Marine04  calibration  (Hughen 
et  ah,  2004b).  I  use  AR=400y,  corresponding  to  a  reservoir  age  of  surface  water  of  800 y. 
I  report  the  median  calendar  ages  with  the  upper  and  lower  limits  of  the  a  and  2a  ranges. 
For  the  measured  ages  older  than  25,000  14 C  y,  I  report  the  intercept  of  this  age  with  the 
G.  bulloides  record  from  the  Cariaco  Basin  (Hughen  et  ah,  2004a),  which  is  tuned  to  the 
GRIP  ss09sea  age  model.  Uncertainty  in  these  calibrations  is  on  the  order  of  1-2  ky. 


4.3  Results 


The  age  model  for  17JPC  (Figure  4-la)  is  preliminary,  as  most  of  this  record  is  constrained 
by  planktonic  14C  dates  whose  calibration  to  calendar  years  is  fairly  uncertain  beyond 
25,000 y  BP  (Table  4.1).  As  a  result,  I  am  unable  to  correlate  events  in  this  record  with 
other,  well-dated  paleoclimate  records  with  any  confidence.  The  average  sample  spacing  in 
17JPC  core  (230  y)  is  theoretically  small  enough  to  resolve  millennial-scale  events  if  present. 

In  the  raw  isotope  data  (Figure  4-1),  there  is  broad  variability  in  TV.  pachyderma  (s.) 
61S0  with  an  amplitude  of  ~0.4  %o,  and  local  maxima  at  352  cm,  410  cm,  520  cm,  and  642  cm. 
There  is  one  notable  maximum  in  Uvigerina  sp .  Sxs 0,  defined  by  3  data  points  at  280  cm, 
with  an  amplitude  of  0.4  %o.  There  is  broad  variability  in  the  rest  of  the  Uvigerina  sp .  5180 
record,  with  an  amplitude  of  ~0.4%c.  Uvigerina  sp.  8l3C  is  variable  with  an  amplitude  of 
up  to  1  %o,  with  minima  at  328  cm,  492  cm,  and  594  cm.  There  are  statistically  significant 
but  small  changes  in  relative  abundance  of  TV.  pachyderma  (s.)  in  this  core. 


4.4  Discussion 


We  plot  the  HLY02-02-17JPC  foraminifer  data  sets  against  age  in  Figure  4-2.  Considering 
the  large  uncertainties  in  the  age  model,  I  plot  the  51S O  data  without  attempting  to  subtract 
whole-ocean  <S180  change,  which  is  on  the  order  of  1  %o  during  MIS3  (Waelbroeck  et  ah, 
2002).  I  will  focus  on  the  relationships  between  features  of  the  isotope  records  versus  depth 
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Figure  4-1:  (a)  Age  model  for  HLY02-02-17JPC,  constructed  assuming  constant  sediment 
accumulation  rate  between  measured  N.  pachyderma  (s.)  radiocarbon  ages,  which  are  cal¬ 
ibrated  to  calendar  years  with  Calib  5.0.1,  using  A R  —  400  y,  or  a  surface  reservoir  age  of 
R  +  AjR  —  800 y.  (b)  Sediment  accumulation  rate  (cm/ky)  inferred  from  the  age  model. 
(c.)-(d)  HLY02-02-17JPC  N.  pachyderma  (s.)  (black)  and  Uvigerina  sp.  (green)  6iS 0.  Dots 
are  individual  measurements,  the  solid  line  connects  the  mean  of  the  measurements  in  each 
sample  (n=3).  (e)-(f)  N.  pachyderma  (s.)  (black)  and  Uvigerina  sp.  (green)  <513C.  (g) 
Marked  with  triangles,  sample  depths  in  which  N.  pachyderma  (s.)  radiocarbon  was  mea¬ 
sured.  (h)  Fraction  relative  abundance  of  N.  pachyderma  (s.),  with  95%  confidence  limits. 
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age  (ky) 

Figure  4-2:  (a)  5180iCe,  of  GISP2  (red,  Stuiver  et  al.  (1995))  and  GRIP  on  the  ss09sea 
age  model  (orange,  Johnsen  et  al.  (2001)).  Interstadials  (Dansgaard-Oeschger  events)  are 
numbered  in  red  and  orange.  Grey  bars  mark  Heinrich  events,  which  are  numbered  in 
black,  (b)  Hulu  stalagmite  6IS O,  a  proxy  for  East  Asian  monsoon  strength,  where  lower 
<5180  indicates  a  stronger  monsoon  (Wang  et  al.,  2001).  Triangles  are  230Th  dates,  cyan 
data  are  from  speleothem  PD,  blue  from  MSD  (Wang  et  al.,  2001).  (c)-(d)  HLY02-02- 
17JPC7V.  pachyderma  (s.)  (black)  and  Uvigerina  sp.  (green)  5I80.  The  grey  shading  is  the 
95%  confidence  limits  on  the  value  of  the  mean  in  each  sample,  (e)-(f)  N.  pachyderma  (s.) 
(black)  and  Uvigerina  sp .  (green)  <P3C.  (g)  Triangles  mark  the  location  of  radiocarbon 
measurements  on  N.  pachyderma  (s.).  (i)  Fraction  relative  abundance  of  N.  pachyderma  (s.), 
with  95%  confidence  limits. 
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(Figure  4-1)  within  the  sediment  core.  In  Figure  4-2,  I  plot  the  GRIP  5180ice  on  the  ss09sea 
age  model  (Johnsen  et  ah,  2001).  The  authors  consider  ss09sea  more  accurate  than  GISP2 
ice  core  age  model  older  than  40  ky  BP,  but  less  accurate  younger  than  40  ky  BP. 

In  the  Bering  Sea  during  the  Bplling-Allerpd,  there  is  evidence  for  relatively  fresh  and 
warm  surface  water  from  low  N.  pachyderma  (s.)  dl80  and  low  relative  abundance  of 
N.  pachyderma  (s.)  (see  Chapter  3).  Sediments  from  within  the  oxygen  minimum  zone  are 
laminated.  In  contrast,  during  the  Younger  Dryas,  sediments  are  oxic,  there  is  indication 
that  surface  water  was  relatively  salty  and  cold.  The  Younger  Dryas  and  Bplling-Allerpd 
are  unique  events  because  they  occur  during  the  deglaciation,  when  the  global  climate  sys¬ 
tem  was  undergoing  a  large  transition.  But  it  is  not  unreasonable  to  hypothesize  that  the 
previous  stadials  and  interstadials  may  be  analogous  to  these  two  events,  because  the  sense 
of  change  in  the  North  Atlantic  is  the  same.  If  this  is  the  case  for  the  North  Pacific,  then 
I  would  predict  that  warm  interstadial  events  are  associated  with  low  planktonic  6180  and 
low  relative  abundance  of  N.  pachyderma  (s.).  Conversely,  I  would  predict  cold  stadials  and 
Heinrich  events  to  be  associated  with  high  planktonic  £180  and  high  relative  abundance 
N.  pachyderma  (s.).  I  expect  that  the  Heinrich  events  and  interstadials  8  and  12  are  the 
most  likely  of  these  events  to  be  observed  in  this  core  because  of  their  severity  and  longer 
duration,  respectively.  Uvigerina  is  not  the  ideal  benthic  species  to  use  in  gauging  the  <513C 
of  bottom  water,  since  it  lives  within  the  surface  sediments,  so  its  5l3C  reflects  pore-water 
S13 C.  However,  it  was  the  only  benthic  foraminifer  species  that  is  present  throughout  the 
core.  Low  Uvigerina  <513C  would  indicate  low  S13 C  of  porewater  which  could  be  a  function 
of  low  bottom- water  5l3C ,  high  organic  carbon  flux  to  the  seafloor,  or  low  oxygen  levels  in 
bottom  water. 

In  Figure  4-1,  though  each  record  is  variable  through  MIS2  and  3,  there  does  not  seem 
to  be  a  consistent  relationship  between  changes  in  any  of  the  records.  For  example,  the  low 
benthic  £13C  at  around  594  cm  coincides  with  relatively  high  benthic  6180  and  relatively 
low  planktonic  5180.  But  the  planktonic  and  benthic  <5180  records  are  neither  at  maxima 
or  minima  while  benthic  S13 C  is  low  at  492  cm.  It  is  possible  that  the  isotopic  records 
represent  changes  in  the  surface  and  deep  ocean  which  are  independent  of  each  other  and 
not  coherent. 

The  sediment  accumulation  rate  is  1.7  times  greater  in  17JPC  during  MIS3  than  the 
Kiefer  et  al.  (2001)  record  from  the  Emperor  Seamounts  (the  sampling  interval  of  the  cores 
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is  the  same).  If  the  temperature  of  surface  water  at  the  two  sites  changed  in  the  same  way  in 
the  past,  I  would  expect  that  the  signal  observed  by  Kiefer  et  al.  (2001)  would  be  more  well 
resolved  in  17JPC.  The  average  relative  abundance  of  N .  pachyderma  (s.)  is  similar  at  the 
two  sites  (~90%),  indicating  that  the  coldest  sea  surface  temperatures  during  MIS3  were 
similar  at  the  two  sites.  However,  the  17JPC  data  do  not  show  millennial-scale  warmings 
of  surface  water  during  MIS3. 

Analysis  of  upper-ocean  temperature  measurements  from  recent  decades  reveals  a  pat¬ 
tern  of  sea-surface  and  upper-ocean  temperature  changes  during  interdecadal  climate  events 
in  the  North  Pacific  (Deser  et  al.,  1996).  In  this  analysis,  the  temperature  changes  are  spa¬ 
tially  heterogeneous,  where  SST  in  the  central  North  Pacific  is  out  of  phase  with  SST  in  the 
Gulf  of  Alaska,  Bering  Sea,  and  eastern  North  Pacific.  If  a  similar  spatial  pattern  of  SST 
anomalies  occurred  during  climate  events  of  MIS3,  there  could  be  a  significant  temperature 
gradient  between  the  Emperor  Seamounts  (in  the  central  North  Pacific)  and  the  Bowers 
Ridge  so  that  the  SST  proxy  records  at  the  two  sites  would  be  dissimilar. 

4.5  Conclusions 

The  preliminary  age  model  from  this  core  does  not  allow  us  to  correlate  features  in  the 
isotopic  and  faunal  records  to  other  well-dated  climate  records.  There  is  variability  in 
N.  pachyderma  (s.)  and  Uvigerina  sp.  $180  of  an  amplitude  of  ^0.4  %o,  and  variability  of 
Uvigerina  sp.  £13C  of  up  to  1  %o.  There  are  statistically  significiant  changes  in  fraction 
relative  abundance  of  N.  pachyderma  (s.)  of  up  to  15%.  However,  there  does  not  appear  to 
be  any  coherence  between  the  features  in  the  oxygen  and  carbon  isotopic  and  faunal  data. 
These  data  suggest,  contrary  to  a  study  from  the  Emperor  Seamounts,  that  there  was  not 
any  millennial  scale  change  in  SST  that  coincide  with  Dasgaard-Oeschger  or  Heinrich  events 
during  MIS3.  However,  spatial  heterogeneity  in  North  Pacific  SST  could  account  for  the 
dissimilarity  of  the  Bowers  Ridge  and  Emperor  Seamount  records. 
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Chapter  5 


Evidence  of  Methane  in  Bering  Sea 
Sediments  during  Marine  Isotope 
Stage  3 


Abstract 

During  Marine  Isotope  Stage  3  (MIS3,  28-62  ky  BP),  there  are  very  large  negative  excursions 
in  planktonic  and  benthic  <513C  in  a  1467  m  core  from  the  Southeast  Bering  Sea  near  the 
Umnak  Plateau,  a  region  where  methane  hydrates  are  present  today.  Measured  benthic  and 
planktonic  foraminifer  6VS C  during  the  excursions  are  as  low  as  —6  and  — 14  %o ,  respectively, 
and  are  the  probably  the  result  of  overgrowths  of  diagenetic  carbonate  minerals.  These 
minerals  may  be  formed  in  association  with  sulfate  reduction  of  organic  matter  or  anaerobic 
oxidation  of  methane  at  the  base  of  the  sulfate  reduction  zone.  Both  processes  increase  the 
alkalinity  and  supersaturate  the  porewater  with  respect  to  carbonate  minerals.  From  a  linear 
end-member  mixing  model,  I  estimate  that  the  overgrowth  stable  isotopic  composition  is 
S13C  =  —  23.4%o,  5l80  =  6.5  %o,  and  [Mg]/[Ca]  =  72-837 mmol/mol.  This  site  is  deep  and 
cold  enough  to  be  well  within  the  zone  of  methane-hydrate  stability.  Thus,  it  is  unlikely  that 
the  source  of  methane  was  from  local  destabilization  of  methane  hydrates  due  to  changes 
in  pressure  or  temperature  at  the  sediment-water  interface. 


5.1  Introduction 

Episodes  of  low  <513C  measured  in  benthic  foraminifera  (as  low  as  —5  %o)  from  ODP  Site  893 
in  the  Santa  Barbara  Basin  coincide  with  warm  Dansgaard-Oeschger  (DO)  climate  events  of 
the  last  glacial  period  (Kennett  et  al.,  2000).  Kennett  proposed  a  “Clathrate  Gun  Hypoth¬ 
esis”  (Kennett  et  al.,  2003)  in  which  a  warming  of  intermediate  water  destabilizes  methane 
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(CH4)  hydrates  (the  term  hydrates  is  synonymous  with  clathrates)  during  interstadials  (IS). 
According  to  this  controvertial  hypothesis,  methane  (a  strong  greenhouse  gas),  is  released 
into  the  water  column,  and  a  significant  fraction  of  the  gas  reaches  the  atmosphere,  enhanc¬ 
ing  global  warming.  There  is  sedimentary  evidence  of  episodic  flux  of  methane  to  surface 
sediments  at  continental  margin  sites,  not  only  in  the  Santa  Barbara  Basin  (low  benthic 
<513C  (Kennett  et  ah,  2000);  biomarkers  of  aerobic  and  anaerobic  methanotrophy,  (Hinrichs 
et  ah,  2003)),  but  also  in  the  Gulf  of  California  (low  benthic  and  planktonic  613C;  (Keigwin, 
2002),  the  Japan  Margin  (methane  biomarkers;  Uchida  et  ah  (2004)),  the  Papua  Gulf  (low 
planktonic  <513 C;  de  Garidel-Thoron  et  ah  (2004)),  and  the  southwest  Greenland  Sea  (low 
benthic  and  planktonic  <513C;  Millo  et  ah  (2005)). 

In  this  paper,  from  a  core  in  the  southeast  Bering  Sea,  I  identify  three  episodes  during 
MIS3  where  there  is  evidence  of  authigenic  carbonate  minerals.  The  minerals  could  be 
associated  with  high  pore- water  alkalinity  from  the  anaerobic  oxidation  of  methane.  If  this 
was  the  case,  however,  this  site  (1467  m)  is  well  within  the  zone  of  methane  hydrate  stability, 
and  does  not  fit  into  the  mechanism  proposed  by  Kennett  et  ah  (2003),  where  warming  of 
intermediate  water  causes  destabilization  of  methane  hydrates.  If  the  events  in  this  core 
coincide  with  DO  events,  these  findings  support  the  proposition  of  Kennett  et  ah  (2003) 
that  methane  stored  in  continental  shelf  sediments  can  be  mobilized  and  interact  with  the 
climate  system  at  millennial  time  scales. 


Methane  hydrates 

In  organic-carbon-  (C^)  rich  marine  sediments,  when  more  energetically  favorable  oxidants 
are  depleted,  is  reduced  to  CH4  by  anaerobic  methanogenic  archaea  (Berner,  1980). 

2CH20  CH4  +  C02 

The  <$13C  of  this  biogenic  CH4  is  typically  —110  to  —  50  %o  because  of  strong  kinetic  frac¬ 
tionation  during  C02  reduction  (Whiticar,  1999).  In  some  settings,  an  additional  source  of 
CH4  in  marine  sediments  is  thermogenic,  where  CH4  is  one  product  of  thermal  degradation 
of  Corg  at  depth  in  the  sediment  column  where  the  temperature  reaches  80-90°C  (Whiticar, 
1999).  The  613C  of  thermogenic  CH4  is  typically  —50  to  —  20  %o,  and  the  molecular  ratio  of 
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Figure  5-1:  Phase  diagram  for  CH4  hydrate  stability.  The  solid  line  is  for  pure  water,  the 
dashed  line  is  for  seawater.  Above  and  to  the  right  of  these  lines  (lower  pressure  and  higher 
temperature),  methane  hydrate  is  not  stable.  Below  and  to  the  left  of  the  lines,  methane 
hydrate  is  stable,  and  with  sufficient  quantities  of  dissolved  methane,  hydrates  will  form. 
The  HLY02-02-51 JPC  is  at  1467  m  water  depth,  and  today,  the  bottom  water  temperature 
is  2.2°C,  well  within  the  zone  of  hydrate  stability.  This  figure  is  from  Kastner  (2001). 


CH4  to  larger  hydrocarbon  molecules,  ethane  (C2H6)  and  propane  (C3H8), 

[CH4]  _  Cl 

[C2H6]  +  [C3H8]  “  C2  +  C3 

is  less  than  50.  In  contrast,  C2+C3  biogenic  CH4  is  typically  greater  than  300  (Kastner, 
2001).  In  regions  of  the  seafloor  where  thermogenic  hydrocarbons  can  migrate  up  along 
faults,  CH4  within  the  sediment  can  be  a  mixture  from  the  two  sources,  in  which  case  its 
composition  falls  on  a  mixing  line  in  C2+C3~^13C  sPace- 

Methane  hydrates  form  on  Earth  today  primarily  in  high-latitude  permafrost  regions, 
and  in  the  ocean  along  continental  margins,  where  there  is  sufficient  CH4,  and  the  proper 
pressure  and  temperature  conditions  (Figure  5-1)  (Kvenvolden  and  Lorenson,  2001).  Gas 
hydrates  are  isometric  crystalline  solids  which  consist  of  a  lattice  of  water  molecules  with 
cages  that  can  accommodate  gas  molecules  that  are  3.5-7.5x  10~10  m  in  diameter  (Pellen- 
barg  and  Max,  2000).  Methane  hydrate  formation  is  associated  with  a  fractionation  of 
oxygen  isotopes,  on  the  same  order  as  fractionation  associated  with  the  formation  of  sea 
ice,  with  a  fractionation  factor,  cx(hydrate-sw)i  m  the  range  1.0027-1.0031,  so  that  hydrate 
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is  enriched  in  H2180  by  3%o  (Matsumoto  and  Borowski,  2000). 

Gas  hydrates  are  stable  at  higher  temperatures  than  ice,  and  with  increasing  pressure, 
the  maximum  temperature  required  for  stability  increases.  The  lower  boundary  of  methane 
hydrate  stability  in  marine  sediments  is  a  function  of  geothermal  heating  from  below.  Dis¬ 
sociated  gas  at  this  boundary  creates  a  seismic  discontintuity  called  a  “Bottom  Simulating 
Reflector,”  (BSR,  Figure  5-2),  a  term  coined  during  the  Deep  Sea  Drilling  Project  (DSDP) 
Leg  XIX  in  the  Bering  Sea  where  it  was  observed  on  the  Umnak  Plateau  at  sites  184 
and  185  (Scholl  et  ah,  1973).  At  site  185,  in  a  sample  from  within  the  BSR  (665mbsf), 
=  2850,  and  <513C  of  CH4  gas  was  —72.1  %o,  both  measurements  indicating  a  purely 
biogenic  source  for  this  gas  (Claypool  et  al.,  1973). 

Pore  water  containing  methane  can  be  advected  upward  at  convergent  margins,  where 
sediment  is  squeezed,  like  in  an  accretionary  wedge,  around  salt  diapirs,  and  up  along  fault 
planes.  At  passive  margins  as  well,  there  can  be  upward  flux  of  pore  fluids  from  sediment 
compaction. 

Authigenic  carbonates  associated  with  anaerobic  methanotrophy 

Biogenic  methane  will  diffuse  upward  in  the  sediment  column  into  the  overlying  sulfate 
reduction  zone.  In  the  presence  of  sulfate,  CH4  is  oxidized  by  syntropic  anaerobic  sulfate- 
reducing  bacteria  and  archaea  (Alperin  and  Reeburgh,  1984). 

CH4  +  SO^  —  2HCO3  +  HS~  +  H20  (5.1) 

This  reaction  produces  HCO3  ,  which  increases  alkalinity  of  porewater,  and  increases  the 
saturation  of  carbonate  minerals.  Sulfate  reduction  of  organic  matter  also  produces  HCO3 
(Figure  5-2). 

2CH2O  +  SO^-  -►  HCO3  +  H2S  (5.2) 

A  large  enough  flux  of  HCO3  into  the  pore  water  will  cause  it  to  become  supersaturated  with 
respect  to  carbonate  minerals.  The  513 C  of  authigenic  carbonates  associated  with  anaerobic 
oxidation  of  methane  will  be  low,  because  the  dissolved  inorganic  carbon  (DIC)  in  porewater 
from  which  the  carbonate  precipitates  is  a  mixture  of  bottom  water  carbon  (typically, 
513C=  -1  to  1  %o),  carbon  from  the  sulfate  reduction  of  organic  matter  (613C=  —  22  %o, 
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Figure  5-2:  Schematic  of  the  sulfate-methane  interface,  with  simplified  reactions  which  occur 
in  the  upper  sediments.  Both  sulfate  reduction  of  organic  matter  and  anaerobic  oxidation 
of  methane  produce  HCOJ,  which  increases  the  alkalinity  of  porewater.  At  the  Umnak 
Plateau,  the  BSR  is  at  670mbsf  (meters  below  seafloor).  At  the  Blake  Ridge,  another 
passive  margin  site  with  sedimentary  methane  hydrates  off  southeastern  North  America, 
the  sulfate-methane  interface  is  20-30  mbsf  (Hoehler  et  al.,  2000),  and  the  upper  limit  of 
hydrates  is  at  ~250mbsf  (Lee,  2000). 
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Rau  et  al.  (1982)),  and  methane-derived  carbon  (513 C=  -110  to  -50 %o).  Depending  on  the 
seep  rate  of  CH4,  intensity  of  bioturbation,  sediment  accumulation  rate,  and  vertical  fluid 
velocity  within  sediment,  the  zone  of  carbonate  precipitation  can  be  very  near  the  sediment- 
water  interface  (Luff  et  al.,  2004).  At  Eel  River  Basin  (off  Northern  California),  Orphan 
et  al.  (2004)  observed  that  small  dispersed  crystals  of  high-Mg  calcite  form  in  sediments 
with  high  CH4  flux,  and  massive  crusts  composed  of  high-Mg  calcite,  aragonite  (the  high- 
pressure  polymorph  of  CaCOs),  and  dolomite  ([Cao.6,Mgo.4]C03)  form  with  low  CH4  flux. 
Aragonite,  which  has  a  more  compact  crystal  structure,  does  not  easily  incorporate  Mg. 
High-Mg  calcite  has  variable  [Mg]/[Ca]  (4-30 mol%  Mg)  (Berner,  1975). 

Methane  seepage  from  the  seafloor  is  observed  to  be  spatially  and  temporally  hetero¬ 
geneous  (Orphan  et  al.,  2004;  Schmidt  et  al.,  2002).  Torres  et  al.  (2003)  measured  <513C 
in  rose-Bengal-stained  U.  peregrina  at  an  active  seep  site  at  Hydrate  Ridge,  off  Oregon. 
They  also  measured  the  <513C  of  the  pore  water  of  the  sediment  in  which  they  were  found. 
Rose  Bengal  stains  proteins  in  cytoplasm,  and  is  used  to  identify  living  or  recently-living 
foraminifera.  They  discovered  that  though  the  513C  of  pore  water  DIC  was  —48  to  —  6%o, 
the  foraminifer  5i3C  calcite  was  —4  to  0%o.  They  hypothesize  that  the  foraminifera  pre¬ 
fer  to  grow  during  episodes  of  low  CH4  discharge,  and  the  513C  of  their  tests  represents 
near-surface  pore  water  which  is  more  like  the  composition  of  bottom-water.  Samples  with 
low  S13C  also  had  high  Mg/Ca,  which  is  characteristic  of  authigenic  CaC03.  Torres  et  al. 
(2003)  thus  interpreted  the  lowest  613C  measured  in  these  foraminifera  as  the  result  of  an 
overprint  of  diagenetic  carbonates  associated  with  anaerobic  oxidation  of  methane. 

In  the  study  by  Torres  et  al.  (2003),  authigenic  overgrowths  were  not  visible  with  electron 
microscopy,  even  though  their  chemical  signature  was  present.  A  small  amount  of  diagenetic 
carbonate  could  have  a  measurable  effect  on  <513C  and  [Mg]/[Ca]  of  a  foraminifer  sample. 
This  because  authigenic  minerals  have  a  very  different  composition  than  foraminifer  calcite; 
Sl3C  would  be  the  same  as  pore-water  <513 C  (at  this  site,  —6  to  —48  %o)  and  [Mg]/[Ca]  would 
be  1-3  orders  of  magnitude  higher  (further  discussion  of  Mg  in  authigenic  minerals  in  Dis¬ 
cussion).  This  explanation  could  account  for  the  contradictory  conclusions  of  another  study 
conducted  on  Hydrate  Ridge  (Hill  et  al.,  2004),  where  the  absence  of  visible  overgrowths  lead 
the  authors  to  conclude  that  the  very  low  Sl3C  measured  in  benthic  foraminifera  reflected 
the  porewater  <513C. 

One  disadvantage  of  using  rose  Bengal  is  that  cytoplasm  can  persist  in  foraminifer  tests 
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for  weeks  to  months;  in  low-oxygen  environments,  it  can  persist  for  even  longer  (Bernhard, 
1988).  So  it  is  not  capable  of  differentiating  living  from  dead  foraminifera.  In  a  study  of 
surface  sediments  from  cold  seeps  in  Monterey  Bay,  Bernhard  et  al.  (2001)  used  a  technique 
which  marks  adenosine  triphosphate,  which  is  only  present  in  living  cells.  They  found  an  as¬ 
semblage  of  benthic  foraminifera  living  in  the  active  cold  seep  which  included  U.  peregrma , 
implying  that  U.  peregrina  is  capable  of  living  in  porewater  with  high  methane  concentra¬ 
tions.  Whether  the  calcite  of  the  U.  peregrina  test  records  bottom  water  613C  or  depleted 
pore-water  £13C  in  areas  with  active  methane  seepage  remains  an  open  question. 

5.2  Methods 

Core  HLY02-02-51 JPC  was  collected  in  June,  2002,  from  the  eastern  Umnak  Plateau  from 
a  water  depth  of  1467m  (See  Figure  3-1).  The  entire  core  was  sampled  every  48cm,  and 
at  every  8  cm  from  450-750  cm.  Raw  sediment  samples  were  dried  at  50  °C,  then  washed 
with  tap  water  through  a  63  pm  sieve.  Planktonic  foraminifer  Neogloboquadrina  pachyderma 
(sinistral)  and  benthic  foraminifer  Uvigemna  peregrina  were  picked  from  the  150-250  pm  and 
>250  pm  size  fractions,  respectively.  Without  further  cleaning,  8-10  N.  pachyderma  (s.)  and 
a  single  U.  peregrina  were  analyzed  on  a  Finnigan  MAT253  mass  spectrometer  with  a  Kiel 
device  (Figure  5-3d-g).  <5180  and  <513C  are  reported  relative  to  the  PDB  standard  via 
NBS-19  following  the  procedures  of  Ostermann  and  Curry  (2000). 

N.  pachyderma  (s.)  were  submitted  to  the  National  Ocean  Sciences  Accelerator  Mass 
Spectrometer  Facility  at  WHOI  for  AMS  14C  analysis.  They  came  from  samples  where 
measured  planktonic  and  benthic  <P3C  were  not  anomalously  low  (less  than  — 2permil). 
The  14 C  dates  (Table  5.1)  were  calibrated  to  calendar  years  using  Calib5.0html  and  the 
Intcal2004  dataset  (Reimer  et  al.,  2004).  I  assume  that  A R,  the  local  anomaly  of  the  age 
of  surface  water  from  the  global  mean,  is  constant  at  400  y.  This  corresponds  to  a  total 
reservoir  correction,  R ,  of  800  y  (see  Chapter  2  for  a  discussion  of  A R).  I  converted  the 
oldest  measured  14 C  date  to  calendar  years  by  finding  the  intercept  of  the  measured  age 
with  a  curve  of  14 C  age  versus  calendar  age  BP  of  Globigerina  bulloides  in  Cariaco  Basin 
(Hughen  et  al.,  2004a).  This  record  is  tuned  to  GRIP  <5180ice  on  the  ss09sea  age  model 
(Johnsen  et  al.,  2001).  I  construct  a  simple  age  model  for  51JPC  by  assuming  constant 
sediment  accumulation  rates  between  calibrated  14C  ages  (Figure  5-3a-b). 
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depth  (cm) 

14C  y  BP±a 

cal  y  BP  +  cr  t  2a 

134.5-135.5 

10,600±60 

7i  1Qn  TT7T5U  11,110 
ii,lUU  11?230  11,300 

177.5-178.5 

12,500±60 

1  2  c:en  13,460  13,400 

10,00U  ,3  530  13,710 

241.5-242.5 

14,050±85 

1  e;  700  15d80  15,310 

1 0 ,  i  UU  15  890  16,120 

479.5-480.5 

18,200±110 

noi  r oy  20,360  20,230 

/U,OCW  20,690  20,950 

1102.5-1105.5 

33,400±780 

37,980 

Table  5.1:  Uncorrected  AMS  14C  ages  and  corrected,  calibrated  ages  if  N.  pachyderma  (s.)  in 
HLY02-02-51JPC.  Calibrations  performed  with  Calib5.0html  (Stuiver  et  al.,  1998)  using  the 
Marine04  calibration  (Hughen  et  al.,  2004b).  I  use  AR=400y,  corresponding  to  a  reservoir 
correction  of  800  y.  I  report  the  median  calendar  ages  with  the  upper  and  lower  limits  of 
the  a  and  2 o  ranges.  For  the  oldest  measured  14C  ages,  I  report  the  intercept  of  this  age 
with  the  G.  bulloides  record  from  the  Cariaco  Basin  (Hughen  et  al.,  2004a),  which  is  tuned 
to  the  GRIP  ss09sea  age  model. 


[Mg]/[Ca]  was  measured  by  ICP-MS  in  two  samples  of  N.  pachyderma  (s.)  consisting 
of  40  individuals  each.  The  samples  underwent  a  full  oxidative  and  reductive  cleaning 
procedure  in  order  to  remove  clay,  metal  oxides,  and  organic  matter  following  the  procedure 
of  Boyle  and  Keigwin  (1985)  (with  revisions  by  Rosenthal  et  al.  (1995)  and  Boyle  and 
Rosenthal  (1996)).  The  only  deviation  from  this  protocol  is  that  sonication  during  the  final, 
weak-acid  leach  was  skipped  because  of  the  fragility  of  the  samples.  Since  [Ca]  measured  in 
the  samples  was  different  than  that  of  the  standard,  these  two  data  points  were  corrected 
for  the  matrix  effect  of  [Ca]  (K.  Dahl,  R.  Came,  unpublished  data). 


5.3  Results 

We  find  large  negative  excursions  in  J13 C  of  both  N.  pachyderma  (s.)  (as  low  as  —13  %o)  and 
U.  peregrina  (as  low  as  —  6.8  %o).  The  <513C  of  N.  pachyderma  (s.)  exhibit  very  large  negative 
excursions  (<C  —  2%o,  the  thin,  grey,  dashed  line  in  Figure  5-3)  at  520-552  cm,  600-624  cm 
and  672-704  cm.  These  intervals  are  denoted  by  vertical  grey  bars  in  the  Figure.  The  mean 
5lsO  and  5l3C  of  U.  peregrina  in  the  interval  450-750  cm  is  calculated,  excluding  samples 
which  lie  within  the  vertical  grey  bars.  These  values  are  shown  in  orange  in  Figure  5-3,  and 
will  be  regarded  as  “typical”  MIS3  values  to  compare  the  isotopic  excursions  to  (-1.39  %o 
for  <513C  and  4.86  %o  for  <S180). 
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depth  (cm) 

Figure  5-3:  (a)  The  age  model  assumes  linear  sediment  accumulation  rate  between  radiocarbon 
dates  measured  on  N.  pachyderma  (s.),  which  were  calibrated  to  calendar  years  using  A R  —  400, 
corresponding  to  a  reservoir  age  of  surface  water,  R  —  800.  (b)  Sediment  accumulation  rate  in  cm/ky. 
(c)  [Mg]/[Ca]  ratio  measured  in  N.  pachyderma  (s.)  in  two  samples,  (d)-(e)  N.  pachyderma  (s.) 
(black)  and  U.  peregrina  (green)  points  are  individual  <5lsO  measurements.  Lines  connect  the  mean 
of  the  points  measured  in  each  sample.  The  green  squares  on  the  y-axis  are  from  a  rose-Bengal- 
stained  U.  peregrina  from  the  top  of  the  companion  multicore  at  this  coring  site.  The  red  diamond 
is  the  equilibrium  <5180  calculated  from  the  T  —  2.20*C  and  S  =  34.50  %o  using  Equation  3.4.  These 
are  the  values  interpolated  for  the  water  column  at  the  depth  of  this  core,  1467  m,  from  hydrocasts 
conducted  during  the  same  cruise,  (unpublished  data  courtesy  of  D.  McCorkle).  Thus,  U.  peregrina 
appears  to  calcify  at  approximately  the  equilibrium  calcite  value,  (f)-(g)  8 13C  measurements.  Colors 
are  the  same  as  for  <5180.  The  grey  dashed  line  is  <513C=  —2.  The  grey  vertical  bars  denote  the  three 
depth  intervals  where  I  have  high-resolution  samples  in  which  <513C  of  N.  pachyderma  (s.)  falls  below 
this  line.  In  (e)  and  (g),  respectively,  the  orange  lines  are  the  mean  U.  peregrina  6ysO  and  d13C 
outside  these  events,  (h)  Location  of  l4C  dates.  '  ' 


depth 

# 

mass 

mean  mass 

[Mg|/[Ca| 

<51:iC 

O, 

oo 

o 

(cm) 

(Mg) 

(Mg) 

(mmol/mol) 

(%o) 

(%o) 

864 

40 

212 

5.3 

71.7 

-12.39 

5.04 

1104 

40 

97 

2.4 

0.65 

-0.56 

3.34 

Table  5.2:  Two  samples  in  which  N.  pachyderma  (s.)  was  analyzed  for  [Mg]/[Ca]  with  ICP- 
MS.  These  values  are  corrected  for  the  matrix  effect  of  [Ca].  5l80  and  613C  are  the  mean 
values  measured  in  N.  pachyderma  (s.)  in  that  sample. 


Figure  5-4:  a)  and  c)  N.  pachyderma  (s.)  and  U.  peregrina  from  the  sample  at  720  cm 
(without  overgrowth);  b)  and  d)  N.  pachyderma  (s.)  and  U.  peregrina  from  the  sample  at 
816  cm  (with  overgrowth).  The  width  of  the  field  of  view  in  each  panel  is  500  /xm.  Note  that 
the  tests  with  overgrowths  are  yellowish.  The  N.  pachyderma  (s.)  from  the  sample  with 
overgrowth  are  noticably  larger,  and  authigenic  material  fills  the  space  between  the  costae 
of  the  U.  peregrina  tests. 
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Figure  5-5:  HLY02-02-51JPC  N.  pachyderma  (s.)  (black)  and  U.  peregrina  (green)  stable 
isotope  data.  Circled  N.  pachyderma  (s.)  data  have  measured  <513C  less  than  -2.  I  interpret 
these  low  5l3C  values  as  the  signature  of  authigenie  carbonate  overgrowths.  The  dashed 
grey  line  is  a  least-squares  fit  to  the  circled  points.  I  hypothesize  that  this  is  a  mixing  line 
between  100%  glacial  foraminifer  calcite  and  100%  carbonate  overgrowth  CaC03.  The  red 
star  is  an  estimate  of  the  composition  of  authigenie  carbonate  based  on  linear  end-member 
mixing  using  the  samples  circled  in  red  (see  text  for  details).  The  mean  d13C  and  Sl80  of 
samples  without  low  dl3C  from  within  the  high-resolution  sampling  interval  are  marked  by 
stars,  grey  for  N.  pachyderma  (s.)  and  orange  for  U.  peregrina .  The  orange  line  is  the  linear 
mixing  line  between  the  red  and  orange  stars. 
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age  (ky) 

Figure  5-6:  HLY02-02-51JPC  data  plotted  against  age.  Black  trianges  indicate  samples 
where  I  measured  radiocarbon  dates  on  N.  pachy derma  (s.).  Data  in  black  are  of  N.  pachy- 
derma  (s.),  data  in  green  are  U.  peregrina.  The  lines  are  the  mean  value  for  each  sample. 
The  mean-ocean  6l80  change  due  to  changing  global  ice  volume  is  subtracted  from  the 
measured  <5180  record  (Waelbroeck  et  al.,  2002).  618Chce  in  the  GRIP  ice  core  (ss09sea.  age 
model)  is  in  red. 
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5.4  Discussion 


The  composition  of  authigenic  carbonate  overgrowths 

Upon  visual  inspection,  the  planktonic  foraminifera  in  the  samples  with  low  planktonic  613C 
appear  yellowish  (Figure  5-4).  In  the  interval  where  I  sampled  closely  (8cm  resolution), 
I  defined  three  episodes  with  anomalously  low  <513C,  each  of  which  is  less  than  1000  y  in 
duration,  according  to  the  preliminary  age  model.  In  Figure  5-5,  TV.  pachy derma  (s.)  Sl3C 
less  than  —2  %o  are  circled,  and  a  linear  least-squares  fit  is  applied  through  the  circled  points. 

SlsO  =  -0.145(513C  +  3.33  (n  =  47; r2  =  0.826; p  <  10“4)  (5.3) 

This  line  has  a  negative  slope,  indicating  that  low  <513C  is  correlated  with  high  5180. 

We  measured  [Mg]/[Ca]  in  TV.  pachyderma  (s.)  in  a  sample  with  low  5l3C  (864  cm),  and 
a  sample  with  <513C  with  the  range  expected  for  Quaternary  benthic  foraminifera  (1104  cm) 
(see  Figure  5-3c  and  Table  5.2).  In  typical  Quaternary  sediments,  foraminifer  [Mg]/[Ca]  is 
1-10  mmol/mol  (Rosenthal  et  ah,  1999),  and  the  data  from  a  sample  at  1104  cm  falls  close  to 
this  range  (0.7  mmol/mol).  At  864  cm,  a  sample  in  which  mean  <513C  of  TV.  pachyderma  (s.) 
is  —  12.4%o,  measured  [Mg]/[Ca]  is  2  orders  of  magnitude  greater  (71.7 mmol/mol).  The 
mass  of  an  average  TV.  pachyderma  (s.)  in  this  sample  is  5.3 /ig,  whereas  in  the  sample  at 
1104  cm,  the  average  mass  is  2 A  fig. 

Not  all  benthic  foraminifera  in  samples  with  evidence  of  authigenic  carbonates  have 
visual  overgrowth  like  TV.  pachyderma  (s.).  It  may  be  that  the  high  surface  area  of  a 
planktonic  foraminifer  test  is  more  likely  to  accomodate  contamination  by  small  crystals 
of  authigenic  carbonate  than  the  imperforate  test  of  U.  peregrina.  For  the  moment,  I 
assume  that  U.  peregrina  did  not  grow  during  times  of  active  CH4  seeping  and  authigenic 
carbonate  formation  (as  in  the  situation  described  by  Torres  et  al.  (2003)),  and  that  the  low 
<5»13  C  measured  in  U.  peregrina  is  also  the  result  of  contamination  by  authigenic  carbonates 
which  occurred  after  the  foraminifer  calcified  its  test. 

We  consider  a  simple  model  to  explain  the  observed  isotopic  and  chemical  composition  of 
the  foraminfer  tests.  In  this  model,  the  isotopic  or  minor  element  composition  of  foraminifer 
tests  contaminated  with  authigenic  carbonates  falls  on  a  linear  mixing  curve  between  the 
compositions  of  primary  test  calcite  and  authigenic  carbonate  (Schmidt  et  al.,  2002;  Orphan 
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et  ah,  2004) 


a18o 

u  vy meas 

=  S18Oforam4>  +  Sl8Oauth(l  -  </>) 

(5.4) 

513Cmeas 

=  <513C  foram*}*  +  <^13Caut/i(l  —  0) 

(5.5) 

[Mg] 

meas 

[Mg]  ,  ,  [Mg]  (.  ,, 

(5.6) 

where  (p  is  the  fraction  of  the  mass  of  a  sample  that  is  foraminifera  CaC03. 


We  can  estimate  the  composition  of  the  authigenic  carbonate  crust  using  the  S18 O, 
513C,  and  [Mg]/[Ca]  measurements  on  the  two  samples  of  N.  pachyderma  (s.)  in  Table  5.2. 
I  assume  that  the  two  samples,  which  are  each  composed  of  40  individual  N.  pachyderma  (s.) 
have  the  same  mass  of  primary  foraminifer  calcite  with  the  same  <513C,  5180,  and  [Mg]/[Ca], 
i.e.,  that  the  only  difference  between  the  two  samples  is  that  the  sample  from  864  cm  has 
authigenic  calcite  added  to  it.  Then  the  sample  at  1104  cm  is  100%  foram  calcite  and  the 
fraction  foraminifer  calcite  in  the  sample  from  864cm  is  <p  =  =  46%.  Data  from 

the  two  samples  are  circled  in  red  in  Figure  5-5.  With  [Mg]/[Ca]mea5  =  71.7  mmol/mol, 
and  [Mgl/tCaj/orarn  —  0.65  mmol/mol  (Table  5.2),  I  can  then  apply  Equation  5.6  to  find 

[Mg]/[Ca]out*: 


[Mg] 

[^a]  auth 


foram ^ 


(W) 


=  132  mmol/mol. 


Similarly,  with  <513C meas  —  — 12.39  %o,  and  <513C foram  —  ^0.56  %o  (the  mean  of  two  mea¬ 
surements)  (Table  5.2),  I  can  apply  Equation  5.5  to  find  5l8Cauth' 


s13  Cauth 


5l3C 


meas 


<5^  C  f  oramtp 


(W) 


-22.37  %o. 


With  £18Omeas  =  5.04  %o,  and  513C/aram  =  3.34  %o  (the  mean  of  two  measurements)  (Ta¬ 
ble  5.2),  I  can  apply  Equation  5.4  to  find  8l8Oauth'- 

rl8^  ^  Omeas  &  & f oram(P  ct  ac  0/ 

S  O  auth  =  - 77 - rH -  =  6.46  %o. 

(1  -(p) 

These  values  of  S13Cauth  and  5IS Oauth  are  plotted  as  a  red  star  in  Figure  5-5.  The  theoretical 
mixing  line  for  benthic  foraminifer  connects  the  authigenic  carbonate  end  member  to  the 
mean  benthic  <513C  and  <5180  in  samples  from  450-740  cm  excluding  those  pertaining  to  low 
£13C  events  (illustrated  in  orange  in  Figure  5-3e-g).  The  few  benthic  data  scatter  near  this 
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theoretical  mixing  line,  which  is  an  encouraging  sign  that  this  model  is  reasonable. 

In  calculating  these  mixing  lines,  I  assume  that  the  composition  of  the  authigenic- 
carbonate  and  foraminifer-calcite  end  members  during  each  separate  episode  were  identical, 
which  may  not  have  been  the  case.  Variability  in  these  end  members  from  event  to  event 
could  account  for  some  of  the  scatter  in  the  data. 

The  estimate  of  [Mg]/[Ca]  of  overgrowth  could  be  affected  by  the  cleaning  procedure  for 
ICP-MS  analysis  (the  same  chemical  cleaning  procedure  is  not  applied  to  samples  before 
analyzing  for  513C  and  £180).  Fifty  percent  and  95%  of  the  initial  mass  of  the  samples  at 
864  cm  and  1104  cm  were  lost,  respectively,  during  cleaning.  For  this  procedure,  loss  of  50% 
of  the  mass  of  a  sample  is  considered  typical.  If  the  authigenic  carbonate  is  more  susceptible 
to  loss,  then  authigenic  [Mg]/[Ca]  is  underestimated  (it  could  be  as  low  as  72  mmol/mol).  If 
the  foraminifer  calcite  is  more  suceptible  to  loss,  then  authigenic  [Mg]/[Ca]  is  overestimated 
(it  could  be  as  high  as  837  mmol/mol).  High-Mg  calcite  and  aragonite  are  more  soluble  than 
pure  calcite.  However,  in  porewater  that  is  undersaturated  with  respect  to  each  of  these 
minerals,  the  phase  most  suceptible  to  dissolution  may  be  a  function  of  the  kinetics  of 
dissolution  rather  than  thermodynamic  stability;  if  crystals  of  the  less  soluble  mineral  have 
a  higher  ratio  of  surface  area  to  volume,  they  may  actually  dissolve  sooner. 

With  the  data  available,  I  can  not  determine  the  relative  contribution  of  sulfate  reduction 
of  organic  matter  and  anaerobic  methane  oxidation  to  the  formation  of  the  authigenic 
minerals.  The  estimate  of  the  Sl3C  of  carbonate  overgrowths,  — 22.4%o,  can  be  solely  a 
product  of  C org  degradation  by  sulfate  reduction  (Equation  5.2),  which  also  raises  alkalinity 
of  porewater,  without  the  need  to  invoke  methanotrophy.  In  this  case,  nearly  all  the  DIC 
in  the  porewater  could  be  from  reduction  of  marine  organic  matter.  The  porewater  DIC  of 
-22.4  %o  could  also  be  a  combination  of  seawater  DIC  (<513 C=  —1  to  l%o),  sulfate  reduction 
of  organic  matter  (<513C=  —22  %o,  Rau  et  al.  (1982)),  and  methane-derived  carbon  (Sl3C= 
-110  to  -50 %o).  One  test  of  the  relative  contribution  of  these  two  processes  would  be  to 
measure  the  A14C  of  the  overgrowth.  Carbon  from  sedimentary  methane  would  be  14C-free, 
but  fresh  organic  matter  contains  the  same  14C  as  surface  water. 

The  mixing  line  between  benthic  foraminifer  data  points  and  the  estimated  concentration 
of  authigenic  carbonate  (the  orange  dashed  line  in  Figure  5-5)  was  calculated  assuming  that 
the  benthic  foraminifera  did  not  grow  during  active  CH4  seeping  and  carbonate  formation, 
and  that  the  5}3C  of  its  shell  is  that  of  a  typical  benthic  foraminifer.  At  modern  CH4  seep 
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sites,  the  upward  flux  of  CH^-rich  porewater  and  associated  authigenic  carbonate  formation 
is  observed  to  be  very  localized  and  spatially  and  temporally  heterogeneous  (Orphan  et  al., 
2004;  Schmidt  et  al.,  2002),  so  that  the  sedimentary  record  of  a  ~1000y  episode  of  CH4 
degassing  at  one  site  may  record  multiple  episodes  of  active  CH4  venting  and  quiesence.  So 
if  this  were  the  case  at  the  51JPC  core  site,  the  primary  calcite  5U C  may  record  a  Sl3C  of 
pore  water  during  an  episode  without  active  methane  seeping,  when  the  foraminifera  are 
more  likely  to  grow.  And  during  an  episode  of  CH4  venting,  the  authigenic  carbonate  would 
form  with  a  lower  pore  water  5l3C. 

The  theoretical  6IS O  of  authigenic  carbonate  is  1.6%o  higher  than  mean  Sis O  of  U.  pere- 
grina  measured  outside  low  513C  events  in  samples  450-750  cm.  According  to  Equation  3.1 
(Chapter  3),  U.  peregrina  calcifies  near  equilibrium  today.  The  higher  £180  of  the  authi¬ 
genic  carbonate  could  be  from  1)  a  change  in  porewater  £180,  2)  a  thermodynamic  effect, 
such  as  change  in  temperature,  or  a  difference  in  the  mineralogy  between  the  overgrowth 
and  the  foraminifer  calcite,  and/or  3)  a  kinetic  effect,  where  the  [Mg]  of  the  porewater 
and/or  authigenic  calcite  affects  growth  rate  and  incorporation  of  oxygen  isotopes  (Berner, 
1975;  Tarutani  et  al.,  1969). 

The  porewater  <5180  is  conceivably  influenced  by  oceanic  basalt  at  the  bedrock-sediment 
interface  and  the  numerous  ash  layers  in  this  sediment  core  from  nearby  Aleutian  back-arc 
volcanics.  However,  fluids  which  have  circulated  through  these  volcanic  rocks  have  low 
5180,  and  could  not  account  for  the  anomalously  high  <5180  of  authigenic  minerals.  The 
5180  of  water  in  CH4  hydrates  is  enriched  by  ~  3  %o  (Matsumoto  and  Borowski,  2000).  The 
presence  of  water  from  the  decomposition  of  hydrates  could  account  for  some  part  of  the 
enrichment  of  pore  water  5IS O  and  authigenic  minerals.  But  it  probably  does  not  account 
for  all  of  it;  a  1.6  %o  enrichment  requires  the  pore  water  to  be  50%  hydrate  water. 

The  mineralogy  of  the  overgrowths  in  this  core  is  unknown,  but  in  modern  cold  seeps, 
the  authigenic  minerals  can  be  composed  of  high-Mg  calcite,  aragonite  and/or  dolomite 
(Luff  et  al.,  2004).  The  fractionation  of  oxygen  isotopes  as  a  function  of  temperature  has 
been  determined  experimentally  by  Kim  and  O’Neil  (1997)  for  calcite  (at  10-40°C),  by 
Kim  et  al.  (2004)  for  aragonite  (at  5-40°C),  and  by  Fritz  and  Smith  (1970)  for  dolomite 
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Figure  5-7:  In  grey  is  the  equilibrium  <5180  of  calcite  as  a  function  of  temperature  (if 
[Ca]o  =  10 mM,  see  Chapter  2  for  discussion)  (Kim  and  O'Neil,  1997).  In  orange  is  the 
equilibrium  (5180  of  aragonite  (Kim  et  ah,  2004).  At  low  temperatures,  these  two  curves  are 
very  similar.  In  green  is  equilibrium  6l80  of  dolomite  ([Caao^Mg^JCOg)  (Fritz  and  Smith, 
1970). 
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([Ca.6oMg.4o]C03,  at  25-79°C).  The  fractionation  factors,  c*(xx-h20) 


a(XX-H20) 


(ts§)xX 

(t5q)h20 


are 


10001na(calcite_H2o)  =  17.32  x  lO3!^1  -  29.00 
1000  In  C*(aragonite-H2o)  =  19.87  X  103  T  1  —  39.09. 
1000  In  Qi(dolomite— h2o)  =  2.62  X  10 6  T  2  +2.17 


and  are  plotted  in  Figure  5-7,  where  all  three  curves  are  extrapolated  to  temperatures 
down  to  — 2°C.  At  low  temperatures,  the  5lsO  of  aragonite  and  calcite  are  very  similar, 
but  dolomite  6180  is  on  average  2.6  %o  higher.  The  presence  of  authigenic  dolomite  in  the 
overgrowths  could  account  for  its  enriched  5180.  The  [Mg]/[Ca]  of  dolomite  is  667  mmol/mol, 
which  falls  in  the  upper  range  of  [Mg]/[Ca]  I  estimated  possible  for  the  overgrowths  (72- 
837  mmol/mol). 

In  51JPC  (1467 m),  benthic  5180  changes  from  4.9  %o  (mean  glacial  value  for  individuals 
not  affected  by  authigenic  overgrowths)  to  3.5  %o  (modern,  from  a  rose-Bengal  stained  in¬ 
dividual  from  the  coretop)  (see  Figure  5-3).  If  the  glacial-interglacial  change  in  deep-water 
518Oh20  is  — 1.1  %o  (Schrag  et  al.,  2002),  the  additional  0.3  %o  decrease,  if  entirely  due  to 
temperature,  means  that  the  deep  water  at  this  site  was  1.1°C  cooler  (with  Equation  3.2). 
Today,  the  bottom  water  at  this  coring  site  is  2.2 °C  (D.  McCorkle,  unpublished  data).  From 
the  equation  of  state  of  seawater  (FofonofF,  1985)  and  estimates  of  temperature,  salinity  and 
pressure  at  the  coring  site,  the  freezing  point  at  the  sediment- water  interface  at  this  site  was 
— 2.9°C  at  MIS3  (taking  into  account  lower  sea  level).  So  the  maximum  cooling  physically 
possible  during  MIS3  is  [(2.2  —  1.1)  —  (—2.9)]  =  4°C.  If  the  1.6  %o  enrichment  of  5180  was 
entirely  due  to  temperature,  it  would  be  a  5.7°C  cooling  (Equation  3.2).  Therefore,  theoret¬ 
ically,  an  inexplicable  cooling  of  bottom  water  during  these  isotopic  events  could  account 
for  some,  but  not  all  of  the  enriched  5IS O.  Adkins  et  al.  (2002)  estimate  the  glacial  deep 
Pacific  ocean  (at  ODP1123,  at  3290  m  water  depth)  temperature  to  be  2.5  ±  0.5°C  colder 
than  today. 

Finally,  kinetic  factors  affect  the  <5180  of  carbonate  minerals.  Berner  (1975)  found 
that  calcite  which  precipitates  inorganically  from  seawater  contains  2-7  mol%  Mg  (20- 
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Figure  5-8:  The  expected  change  in  5180  of  magnesian  calcite  as  a  function  of  mol%  Mg 
(empirical  relationship  from  Tarutani  et  al.  (1969)). 


75  mmol/mol).  As  Mg  is  incorporated  into  the  calcite  crystal  lattice,  it  causes  the  growth 
surface  to  have  higher  solubility  than  pure  calcite,  and  this  retards  the  growth  rate.  The 
61S0  of  the  Mg-containing  (magnesian)  calcite  is  also  affected,  where 

1000  In  a(MC-H20)  =  0.0776m  +  27.61  (5.7) 

where  MC  is  magnesian  calcite,  and  m  is  the  mol%  Mg  of  the  MC  (Tarutani  et  al.  (1969), 
experiment  conducted  at  25°C).  This  relationship  is  nearly  linear,  where  5180  increases 
0.078  %o  for  every  1  mol%  Mg  increases  (Figure  5-8).  If  I  assume  the  overgrowth  is  all 
magnesian  calcite,  and  I  use  the  estimate  of  the  [Mg]/[Ca]  of  the  overgrowth,  the  expected 
enrichment  of  <5180  (using  Equation  5.7)  is  0.9  %o  (0.5-3. 5  %o).  Enrichment  of  £180  associ¬ 
ated  with  magnesian  calcite  could  also  account  for  some,  but  not  all  of  the  high  <5180  of  the 
overgrowths. 

In  summary,  the  high  <5180  of  the  authigenic  minerals  could  be  accounted  for  if  the 
minerals  are  composed  of  high-Mg  calcite  and/or  dolomite.  The  presence  of  hydrate-derived 
water  would  also  increase  the  6180. 
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The  souce  of  inethane 


If  the  authigenic  minerals  are  associated  with  a  greater  flux  of  methane  to  the  surface 
sediments,  this  raises  a  quandary.  The  coring  site  of  51JPC  is  so  cold  and  deep,  that 
destabilization  of  methane  hydrates  at  this  site  would  require  an  unrealistic  warming  (15°C) 
or  sea-level  drop  (—1100 m).  However,  methane  can  be  transported  within  the  hydrate 
stability  zone  if  (1)  the  concentration  is  not  high  enough  to  form  hydrate  or  (2)  there  is  a 
kinetic  barrier  to  the  formation  of  hydrate.  Brewer  et  al.  (1998)  showed  that  methane  gas 
released  into  seawater  (at  pressure  and  temperature  conditions  within  the  hydrate  stability 
zone)  formed  a  crust  of  hydrate  around  the  periphery  of  the  bubble,  isolating  the  gas  in  the 
interior  of  the  bubble  from  the  seawater,  and  slowing  the  hydrate  formation.  The  positive 
buoyancy  of  the  bubble  then  would  move  it  upward.  In  the  Black  Sea,  there  are  active 
methane  seeps  observed  within  the  hydrate  stability  zone  (Klaucke  et  al.,  2005).  There  may 
be  processes  active  at  the  Umnak  Plateau  in  the  past  which,  like  in  the  Black  Sea,  mobilize 
methane  and  transported  to  the  sediment-water  interface,  but  at  this  time,  this  is  highly 
speculative,  and  if  the  minerals  are  indeed  associated  with  methane,  it  is  a  mystery  what 
the  processes  was. 

Climate  implications 

According  to  the  preliminary  age  model,  the  duration  of  each  of  the  three  episodes  of  low 
<513 C  was  less  than  1000  y,  and  the  events  appeared  to  be  spaced  greater  than  1000  y  apart. 
Though  the  age  model  is  too  preliminary  to  show  if  the  events  in  the  core  coincide  chrono¬ 
logically  with  Dansgaard-Oeschger  (DO)  events,  they  appear  to  have  similar  duration  and 
timing.  If  these  episodes  coincide  with  DO  events,  then  it  supports  the  part  of  the  hypoth¬ 
esis  of  Kennett  et  al.  (2003),  that  sedimentary  methane  hydrates  interact  with  the  climate 
system  on  short  timescales.  However,  the  mechanism  proposed  by  those  investigators  can 
not  explain  the  presence  of  methane  in  a  coring  site  that  is  as  deep  and  cold  as  51JPC. 

5.5  Conclusions  and  Future  Work 

We  find  evidence  of  authigenic  carbonates  with  low  <513C,  high  <5180,  and  high  [Mg]/[Ca] 
which  occur  during  MIS3  in  a  core  from  1500  m  from  the  Umnak  Plateau  region.  These 
overgrowths  may  be  associated  with  anaerobic  oxidation  of  methane.  Future  work  includes 
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sampling  the  rest  of  this  core  at  high  resolution  and  improving  the  age  model.  I  will  also 
measure  bulk  513 C  and  the  concentration  and  <513 C  of  biomarkers  of  anaerobic  and 
aeorobic  methanotrophy,  if  present  during  these  events.  These  measurements  will  help  to 
elucidate  whether  methane  was  involved  in  the  formation  of  the  minerals,  and  if  so,  what 
the  origin  of  the  methane  might  have  been.  I  will  also  make  the  same  measurements  in  a 
shallower  core,  from  within  the  water  depths  where  fluctuations  of  water  temperature  have 
the  potential  to  affect  the  stability  of  sedimentary  gas  hydrates. 
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Chapter  6 


Summary 


The  main  contribution  of  this  thesis  is  that  it  expands  the  geographic  range  of  high- 
resolution  paleoceanographic  observations  to  the  far-North  Pacific.  In  Cook  et  ah  (2005) 
(Chapter  2  of  the  thesis),  I  reported  the  preliminary  stratigraphic  data  from  new  sediment 
cores  collected  during  HLY02-02  in  the  Bering  Sea.  The  most  important  findings  reported 
in  this  paper  are  (1)  that  there  are  sites  in  the  Bering  Sea  with  high  sediment  accumulation 
rate  and  good  calcium  carbonate  preservation,  (2)  that  laminated  and  dysoxic  sediments 
accumulated  at  intermediate  depths,  and  (3)  that  there  are  rapid  (occurring  in  less  than 
220  y)  0.7-0.8%c  decreases  in  N.  pachyderma  (s.)  5iS O  during  the  deglaciation. 

In  Chapter  3,  I  present  high-resolution  stratigraphic  and  faunal  data  of  the  deglaciation 
from  two  Bering  Sea  cores,  one  from  the  Bering  Slope  and  one  from  Bowers  Ridge  (Figure  6- 
1).  I  also  present  a  compilation  of  published  and  unpublished  14 C  data  from  planktonic 
(surface-ocean  dwelling)  and  benthic  (surface-sediment  dwelling)  foraminifera.  I  found  that 
the  vertical  profile  of  14 C  does  not  appear  to  have  been  significantly  different  than  today 
at  the  Last  Glacial  Maximum  (  21  k  BP),  the  Bplling-Allerpd  (BA,  a  brief,  warm  climate 
event  during  the  deglaciation,  14.7-12.9  k  BP)  and  Younger  Dryas  (YD,  a  cold  climate  event 
which  followed  the  BA,  12.9-11.7  k  BP)  (See  Figure  A-l  for  core  locations,  and  Figure  3- 
13  for  data).  The  14C  calculations  are  described  in  detail  in  Appendix  1.  Interpretation 
of  this  dataset  is  limited  by  (1)  the  small  number  of  North  Pacific  data  points  and  (2) 
the  fact  that  many  of  the  benthic  and  planktonic  foraminifer  14 C  measurements  from  this 
dataset  are  probably  affected  by  post-depositional  processes  which  are  known  to  bias  14 C 
measurements.  In  addition  the  simplified  version  of  the  behavior  of  oxygen  and  14 C  in  the 
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Figure  6-1:  The  Bering  Sea.  The  scale  for  color  bathymetry  and  topography  is  in  meters. 
Coring  sites  from  HLY02-02  are  labeled  with  black  dots. 

ocean  as  reviewed  here  is  only  accurate  to  the  first  order,  as  long  as  the  two  mechanisms  I 
described  are  the  only  processes  that  are  changing.  If  I  had  observed  a  significant  change 
in  the  14 C  profile  of  the  North  Pacific,  it  would  be  consistent  with  a  change  in  circulation, 
but  it  would  not  require  it. 

In  the  core  from  the  Bering  Slope,  I  find  that  there  are  very  rapid  (occurring  in  less  than 
220  years)  decreases  in  planktonic  <5180  at  two  times  during  the  deglaciation  corresponding 
to  the  beginning  of  the  BA  warm  period,  and  the  end  of  the  YD  cold  period.  These 
isotopic  shifts  are  each  accompanied  by  a  transient  decline  in  the  relative  abundance  of  the 
cold-water  planktonic  foraminifer  species,  N.  pachyderma  (s.).  Since  the  two  <5180  changes 
are  steps,  and  the  warmings  implied  by  the  decreases  in  abundance  of  N.  pachy derma  (s.) 
are  transient,  I  conclude  that  these  isotopic  events  each  represent  a  combination  of  two 
processes:  a  sudden  and  transient  warming,  and  a  sudden  and  sustained  freshening  of  surface 
water  in  the  Bering  Sea.  I  made  one  simple  test  of  whether  the  faunal  record  was  affected 
by  differential  dissolution  among  the  planktonic  species.  The  test  result  was  negative,  but 
the  inferred  warming  would  be  more  robust  with  another  proxy  for  sea  surface  temperature. 
The  Bering  Slope  core  is  intermittently  laminated,  and  the  sediment  accumulation  rate  is 
very  high  (200cm/ky).  These  two  physical  characteristics  make  the  records  from  this  core 
minimally  affected  by  post-depositional  processes  which  jumble  the  records  from  cores  with 
lower  sediment  accumulation  rates  (like  the  one  from  Bowers  Ridge),  and  therefore  it  is 
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much  more  reliable. 


The  final  important  observation  made  in  Chapter  3  is  that  the  ages  of  the  transitions 
between  laminated  and  not-laminated  sediment  in  the  Bering  Slope  core  correspond  to 
the  ages  of  the  transitions  between  warm  and  cold  events  with  duration  of  centuries  to  a 
millennium  as  they  are  defined  in  an  ice  core  from  the  summit  of  Greenland  (GISP2)  to  well 
within  the  uncertainties  of  the  independent  age  models  of  these  two  records  (see  Figure  3- 
12).  This  correlation  is  robust,  even  considering  the  uncertainty  in  the  apparent  age  of 
surface  water  in  the  Bering  Sea  during  the  deglaciation.  At  this  time,  I  can  not  do  more 
than  speculate  about  a  mechanism  that  connects  the  surface  water  temperature  and  salinity 
in  the  Bering  Sea,  the  oxygen  content  of  1130  m  depth  water  in  the  Bering  Sea,  and  air 
temperature  over  Greenland.  However,  these  observations  will  provide  a  strong  constraint 
to  any  proposed  mechanism  that  connects  the  North  Pacific  and  the  North  Atlantic  on 
centennial  to  millennial  time  scales. 

In  Chapter  4,  I  present  high-resolution  records  of  planktonic  and  benthic  £180  and  613 C, 
as  well  as  relative  abundance  of  N.  pachyderma  (s.)  from  the  core  from  the  Bowers  Ridge 
(2200  m  water  depth)  which  spans  the  period  from  about  50,000  years  ago  to  the  present.  I 
find  no  evidence  of  coherent  millennial-scale  events  in  benthic  <513C  or  N.  pachyderma  (s.) 
relative  abundance  and  6IS O  (Figure  4-2).  These  data  suggest  that  there  was  not  significant 
variability  in  the  surface  water  or  2200  m  water  temperature,  salinity  or  nutrient  content 
from  50,000  to  20,000  years  ago.  Several  investigators  have  suggested  previously  that  there 
were  statistically  significant  changes  in  the  613C  of  intermediate  and  deep  water  coinciding 
with  the  millennial-scale  climate  events  which  are  observed  to  occur  in  the  North  Atlantic 
region  ( i.e .  Lund  and  Mix  (1998);  Mix  et  al.  (1999)),  which  they  interpret  as  a  change  in 
the  rate  of  formation  of  these  water  masses.  In  a  record  of  far-North  Pacific  surface  water 
variability  at  the  Emperor  Seamounts,  Kiefer  et  al.  (2001)  conclude  that  there  were  large 
and  rapid  changes  in  sea  surface  temperatures  which  coincide  with  the  Atlantic  millennial 
climate  events.  The  Bowers  Ridge  record  has  comparable  or  higher  temporal  resolution 
than  these  published  records.  Though  the  age  model  for  the  core  is  very  preliminary,  if  the 
sediment  accumulated  continuously  (and  there  is  no  seismic  or  lithological  evidence  that 
this  was  not  the  case),  events  with  duration  and  spacing  on  the  order  of  a  thousand  or  more 
years  ought  to  be  detectable.  None  are  seen,  which  suggests  that  either  surface  water  at 
the  Emperor  Seamounts  and  Bowers  Ridge  are  not  changing  the  same  way  in  the  past. 
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In  Chapter  5,  I  present  another  set  of  planktonic  and  benthic  5l30  and  S13 C  records 
from  the  last  glacial  period,  this  one  from  the  Umnak  Plateau  region,  in  the  southeast 
Bering  Sea  (1570m  water  depth).  These  records  differ  from  the  ones  presented  in  Chapter 
4  in  that  there  are  extremely  large,  episodic,  negative  excursions  in  benthic  and  planktonic 
<513C.  These  S13 C  excursions  are  not  related  to  the  use  of  <513C  as  a  nutrient  proxy.  The  6l3C 
values  measured  (as  low  as  —  14%oPDB)  are  far  outside  the  range  of  <513C  normally  found  in 
Pleistocene  foraminifera  (-1  to  1%oPDB),  and  are  associated  with  anomalously  high  5180 
(Figure  5-5).  I  hypothesize  that  these  isotopic  features  are  the  result  of  episodic  occurrence 
of  overgrowths  of  carbonate  minerals  on  the  foraminifer  shells.  This  is  supported  by  the 
physical  appearance  of  foraminifera  from  samples  with  anomalous  isotopes  (Figure  5-4), 
which  often  have  visible  yellowish  crusts  on  them,  and  are  significantly  greater  in  mass 
(Table  5.2).  I  constructed  a  simple  end-member  mixing  model  to  estimate  the  isotopic 
composition  and  [Mg]/[Ca]  of  the  overgrowths,  and  find  that  they  are  consistent  with  the 
chemical  signature  of  carbonate  minerals  that  form  inorganically  in  seawater. 

Why  these  overgrowths  occur  is  not  clear.  The  overgrowths  could  conceivably  be  the 
result  of  saturation  of  porewater  with  respect  to  carbonate  minerals  due  to  decomposition 
of  organic  matter  via  sulfate  reduction.  However,  I  speculate  that  the  overgrowths  are 
associated  with  episodic  flux  of  methane-rich  porewater  to  the  surface  sediments.  Inorganic 
carbonate  minerals  frequently  form  at  active  methane  seeps  in  the  modern  ocean.  Also, 
there  is  seismic  evidence  of  methane  deep  in  sediments  (hundreds  of  meters  below  the 
seafloor)  at  this  site.  Based  on  the  preliminary  age  model  for  this  core,  I  can  not  tell  yet 
whether  these  episodes  coincide  with  the  millennial-scale  climate  events  of  the  last  glacial 
period,  but  they  do  have  similar  duration  and  timing  as  those  events.  If  they  are  indeed 
synchronous  with  events,  this  may  imply  that  sedimentary  methane  interacts  with  the 
climate  system  on  short  time  scales. 

Methane  is  a  very  powerful  greenhouse  gas  (more  potent  than  carbon  dioxide),  and 
records  of  atmospheric  methane  concentration  reconstructed  from  ice  cores  show  that  dur¬ 
ing  the  warm  millennial-scale  climate  events,  there  is  more  methane  in  the  atmosphere 
(Chappellaz  et  al.,  ????).  Kennett  et  al.  (2003)  proposed  a  fascinating  but  very  controversial 
hypothesis  dubbed  the  “The  Clathrate  Gun  Hypothesis,”  in  which  warming  of  intermediate- 
depth  water  during  millennial-scale  warm  climate  events  caused  world-wide  destabilization 
of  sedimentary  deposits  of  methane  (in  the  form  of  methane  clathrates).  In  this  model,  some 


94 


of  the  methane  escapes  oxidation  in  the  water  column  and  enters  the  atmosphere,  where 
it  enhances  global  warming.  However,  the  Umnak  Plateau  coring  site  is  cold  and  deep 
enough  that  thermodynamic  destabilization  of  methane  clathrates  is  unlikely  to  occur.  If 
the  authigenic  minerals  are  related  to  methane  in  the  sediments,  it  is  an  open  question  how 
the  methane  was  transported  to  the  surface  sediments. 

Future  directions 

The  work  on  these  Bering  Sea  cores  is  still  in  the  early  stages.  I  have  established  stratigra¬ 
phies  for  the  best  cores  at  each  coring  site,  and  tested  some  of  the  interpretations  of  previous 
investigators  who  have  worked  in  the  North  Pacific.  But  there  are  still  many  open  questions 
that  can  be  addressed  in  the  future:  (1)  What  mechanism  is  responsible  for  the  lithological 
changes  during  the  deglaciation?  There  are  several  other  geochemical  proxies  (for  surface 
ocean  primary  productivity,  porewater  oxygen  content,  and  the  redox  state  of  surface  sedi¬ 
ments)  which  could  shed  light  on  this  question.  As  a  postdoctoral  investigator,  I  will  work 
more  on  this  problem  using  14C  and  613C  of  a  variety  of  benthic  foraminifer  species.  (2) 
On  a  related  note,  what  ties  subsurface  oxygen  in  the  Bering  Sea  to  Greenland  surface  air 
temperature?  Whatever  mechanism  ties  the  two  together  must  be  unaffected  by  the  large 
changes  that  occur  during  the  deglaciation,  including  the  melting  of  the  Laurentide  ice  sheet 
on  North  America,  the  shift  of  the  shoreline  from  the  edge  of  the  continental  slope  to  its 
present  position  hundreds  of  kilometers  away  with  rising  sea  level,  and  opening  of  the  Bering 
Strait.  This  question  could  be  further  explored  in  a  coupled  ocean-atmosphere  model.  (3) 
Where  does  the  freshwater  come  from  that  we  see  evidence  of  during  the  deglaciation?  Con¬ 
ventional  wisdom  says  that  Alaska  was  not  highly  glaciated  at  the  Last  Glacial  Maximum. 
(4)  If  there  was  indeed  methane  in  surface  sediments  of  the  Umnak  Plateau  during  warm 
millennial-scale  climate  events,  how  did  it  get  there?  Could  methane-rich  fluids  migrate 
up  a  fault  from  the  relict  subduction  zone  during  an  earthquake?  To  further  investigate 
this  question  I  will  look  in  this  core  and  at  least  one  other  shallower  core  from  the  Um¬ 
nak  Plateau  and  try  to  establish  whether  these  episodes  were  indeed  related  to  methane 
by  looking  for  organic  biomarkers  which  are  diagnostic  of  methanotrophic  organisms,  and 
whether  they  seem  to  occur  synchronously  in  multiple  cores. 
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Appendix  A 


Apparent  ventilation  age  of  the 
North  Pacific  during  the  last 
deglacial  transition 


Abstract 

I  present  three  ways  to  quantify  the  depletion  of  14 C  in  subsurface  water  relative  to  surface 
water  in  the  past  from  measurements  of  A14C  in  benthic  and  planktonic  foraminifera. 
Each  method  produces  a  subtly  different  quantity.  I  compile  published  and  unpublished 
pairs  of  benthic  and  planktonic  foraminifer  14C  measurements  from  the  North  Pacific  from 
during  the  last  glacial  Termination.  The  data  do  not  show  a  significant,  coherent  change  in 
subsurface  14C  during  the  Bplling-Allerpd  (14.7-12.9  ky  BP)  and  the  Younger  Dryas  (12.9- 
11.7  ky  BP)  and  the  Last  Glacial  Maximum  (23-21  ky  BP)  compared  to  today.  However, 
because  of  the  stratigraphy  and  faunal  abundances  in  North  Pacific  cores,  finding  benthic 
and  planktonic  foraminifera  from  the  same  sample  with  reliable  14 C  is  very  difficult.  This 
null  result  could  be  a  consequence  of  the  limitations  of  this  dataset. 


A.l  Benthic-planktonic  14C  age  difference 

A  conventional  14C  age  is  calculated  from  A14C  measured  on  a  sample  and  the  Libby  14C 
decay  constant,  Xl  —  gpjg - , 


(14Cage)L  =  -ALUn 


A14C  A 
1000  +  )’ 
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where  A14C  is  calculated  relative  to  a  “pre-industrial,  pre-nuclear”  atmosphere  (PIPN): 

Al*c=M4Cb-1) 

Measured  14C/12C  is  corrected  for  (1)  the  effects  of  isotopic  fractionation  and  (2)  so  that 
the  conventional  14 C  age  is  in  years  before  “present”  (BP),  where  the  present  is  defined  as 
the  year  1950  AD.  A  14C  age  calculated  with  the  true  decay  constant,  A  =  8286 - ,  is  1.029 
times  the  magnitude  of  (14Cage)L,  the  ratio  of  A  to  \i 

(“Cage)  =  +  (A.1) 

=  ^(14Cage)L  (A. 2) 

The  difference  between  a  benthic  and  planktonic  14C  age  using  the  true  half  life,  A is 
larger  than  the  difference  using  the  Libby  half  life,  A b-pL,  by  the  same  ratio.  Though  a 
conventional  14 C  age  is  calculated  with  A/,,  I  use  A b~p  because  this  quantity  can  be  directly 
compared  to  other  estimates  of  the  apparent  ventilation  age  of  subsurface  water  in  units  of 
true  calendar  years. 

In  a  compilation  of  pairs  of  benthic  and  planktonic  14C  ages  from  North  Pacific  cores 
(Figure  A-l),  I  plot  A^_p  from  three  time  slices  in  Figure  A-2,  the  Younger  Dryas  (YD, 
11,650-12,890  y  BP),  the  B0lling-Aller0d  (BA,  12,890-14,670  y  BP),  and  the  Last  Glacial 
Maximum  (LGM,  19,000-23,000  y  BP)  (see  Table  A.l).  The  time  periods  for  BA  and  YD 
are  defined  as  the  midpoint  of  the  rapid  transitions  in  5180iCe  from  GISP2  (Stuiver  et  ah, 
1998).  This  record  has  an  independent,  layer-counted  age  model  during  the  deglaciation. 
Closed  symbols  are  data  points  where  the  median  of  the  calibrated  planktonic  calendar  age 
falls  within  the  given  time  slices  with  95%  (2a)  certainty.  The  median  age  of  data  plotted 
with  open  symbols  fall  within  the  time  slice,  but  the  entire  2a  range  do  not.  The  standard 
deviation  of  A^_p  is  the  square  root  of  the  sum  of  squares  of  standard  deviations  of  the 
benthic  and  planktonic  ages  (a  =  (a2  +  a2)1/2),  assuming  that  the  errors  in  these  two  ages 
are  uncorrelated. 

For  reference,  I  plot  A b~p  with  a  profile  of  mean  North  Pacific  A14C  from  the  WOCE 
dataset  (Key  et  ah,  2002).  This  profile  includes  stations  from  0°N  to  55°N.  The  modern  A14C 
profile  is  converted  to  apparent  age  using  Equation  A.l.  This  profile  is  not  corrected  for  the 
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effects  of  fossil  fuel  burning  or  nuclear  tests  on  14C.  In  the  North  Pacific,  the  penetration 
of  CO2  from  these  two  phenomena  is  largely  confined  to  the  uppermost  1  km  and  has  the 
net  effect  of  elevating  A14C  (Key,  2001;  Sarmiento  and  Gruber,  2002). 

The  quantity  that  is  directly  comparable  to  A^_p  is  the  difference  between  the  profile 
of  A14C  converted  to  apparent  age  and  the  apparent  age  of  the  surface  water.  The  global 
mean  apparent  age  of  ocean  surface  water  (R)  is  400  y,  but  it  is  greater  in  areas  with  greater 
upwelling  and  sea  ice  cover  (Bard,  1988).  In  Figure  A-2,  I  subtract  R  from  the  profile  of 
apparent  age  as  the  upper  bound  of  the  shaded  area.  I  subtract  800  y,  a  typical  surface 
water  age  in  an  upwelling  area,  from  the  profile  as  the  lower  bound  of  the  shaded  area. 
This  is  to  illustrate  the  expected  variability  in  A since  the  planktonic  14 C  ages  are  not 
corrected  for  the  deviation  of  local  surface  water  apparent  age  from  the  global  mean  value, 
A  R. 

The  calendar  ages  of  two  new  data  points  from  this  study,  from  HLY02-02-3JPC  and 
HLY02-02-17JPC,  do  not  fall  within  the  BA,  but  just  before.  I  plot  these  data  points  on 
the  BA  time  slices  in  Figures  A-2,  A-4,  and  A-6  because  within  the  sediment  cores,  the 
samples  are  within  the  dysoxic/laminated  stratigraphic  unit  which  I  believe  was  deposited 
during  the  BA.  To  distinguish  them  from  the  other  data,  I  plot  these  points  with  star-shaped 
symbols. 


A. 2  A14C  of  subsurface  water 


Using  the  well-known  rate  of  decay  of  14C,  the  trajectory  of  A14C  of  a  sample  of  benthic 
foraminifera  can  be  “undecayed”  backward  in  time 


A14C(y)  =  1000 


A14Cm 

1000 


eXy  —  1 


(A.3) 


where  A14Cm  is  the  value  measured  in  a  benthic  foraminifer  sample  today,  and  y  is  age  of 
the  sample  in  calendar  years  BP.  This  decay  trajectory  is  plotted  as  the  exponential  black 
line  in  Figure  A-3,  where  the  1  a  and  2 a  measurement  uncertainty  is  propagated  through 
time  and  contoured  in  dark  and  light  grey,  respectively. 

The  A14C  of  the  sample  at  the  time  the  benthic  foraminifera  calcified  their  shells,  A14Co, 
is  estimated  by  finding  the  intersection  of  the  decay  trajectory  with  the  estimated  calendar 
age  of  the  sediment  sample.  In  Figure  A-3,  the  la  and  2 a  contours  of  the  resulting  joint 
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Figure  A-l:  Locations  of  sediment  cores  in  a  compilation  of  published  and  unpublished 
pairs  of  benthic  and  planktonic  14C  ages.  Red  circles  are  coring  sites  shallower  than  1  km 
water  depth,  orange  squares  are  coring  sites  1  -2km,  and  black  diamonds  are  coring  sites 
deeper  than  2  km.  Coring  site  depths  are  not  corrected  for  deglacial  sea  level  change.  Each 
site  is  annotated  with  the  acronyms  used  in  this  chapter.  See  Table  A.l  for  references  and 
details  of  the  cores. 
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Figure  A-2:  Compilation  of  pairs  of  benthic  and  planktonic  foraminifer  14C  ages  from  the 
North  Pacific.  Plotted  here  is  A the  difference  between  the  14 C  age  in  ky  (using  the 
true  half  life)  of  benthics  and  planktonics.  The  average  North  Pacific  A14C  profile  from 
the  WOCE  dataset  (Key  et  ah,  2002)  is  converted  to  apparent  age  with  Equation  A.l.  The 
right-hand  edge  of  the  shaded  region  is  the  apparent  age  profile  minus  400  y  (the  mean 
apparent  age  of  surface  water),  and  the  left-hand  edge  of  the  shaded  region  is  the  apparent 
age  profile  minus  800  y  (a  hypothetical  apparent  age  of  surface  water  in  an  upwelling  region). 
The  time  periods  of  the  YD  and  BA  are  as  defined  in  Stuiver  et  al.  (1998).  For  each  time 
slice,  I  plot  data  for  which  the  median  calibrated  planktonic  calendar  age  falls  within  the 
defined  time  periods.  The  filled  symbols  are  data  points  where  the  calibrated  planktonic 
calendar  age  falls  within  the  time  interval  with  95%  (2a)  certainty.  Table  A.l  contains  the 
data,  references,  and  explanation  of  the  acronyms  used  here. 


probability  distribution  (PDF)  are  drawn  in  black.  The  joint  PDF  is  integrated  along  A14C 
to  find  the  PDF  of  A14Co-  If  A14C  of  benthic  foraminifera  is  the  same  as  that  of  the  bottom 
water  at  the  core  site,  then  A 14 Co  is  an  estimate  of  the  14 C  activity  of  this  water  in  the 
past.  The  PDF  of  A14Co  is  drawn  in  black  on  the  left-hand  y-axis  with  the  la  and  2 a 
ranges. 

In  order  to  compare  the  reconstructed  A14C  activity  of  subsurface  water  to  a  modern 
profile  of  mean  North  Pacific  A14C,  I  subtract  the  estimate  of  the  mean  surface  ocean 
(AMCSW)  (Hughen  et  ah,  2004b)  from  A14C0.  In  Figure  A-4,  I  plot  A14C0-A14Cstt,  for 
the  three  time  slices.  For  reference,  I  plot  in  the  solid  grey  curve  the  Key  et  al.  (2002) 
profile  of  A14C,  where  the  R  (400 y)  is  subtracted  so  that  the  curve  is  relative  to  the  mean 
pre-industrial,  pre-nuclear  surface  ocean. 

The  other  reference,  the  dashed  grey  line  is  what  the  modern  profile  would  look  like 
with  higher  A14C  everywhere.  During  a  time  when  the  atmospheric  and  surface  ocean 
A]4C  were  higher  than  today,  if  the  ocean  circulation  was  precisely  the  same  as  it  is  today, 
the  A  14Cq— A14Csu,  would  be  larger  in  magnitude.  This  is  because  of  the  exponential  decay 
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Figure  A-3:  Schematic  of  calculation  of  Al4Co  and  A14Cq— A14C5lt;.  The  intersection  of 
the  decay  trajectory  of  the  benthic  foraminifer  sample  with  the  calendar  age  of  the  sample 
is  drawn  with  black  contours  at  the  lcr  and  2 o  levels-  The  integral  of  this  joint  probability 
distribution  (PDF)  along  the  y-axis  is  the  PDF  of  A14Co,  which  is  drawn  with  the  la  and 
2a  ranges  on  the  left-hand  axis.  For  comparison  to  a  modern  mean  North  Pacific  A14C 
profile,  the  mean  surface  ocean  A14C  (Hughen  et  ah,  2004b)  is  subtracted  from  this  estimate 
of  A14Co  to  find  A14Co— A14CSUM  the  depletion  of  14C  in  the  benthic  sample  relative  to  the 
global  mean  surface  ocean  value.  This  joint  PDF  is  plotted  in  green,  and  the  integrated 
PDF  is  on  the  right-hand  axis.  The  medians  of  the  joint  PDFs  are  marked  with  stars. 
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Figure  A- 4:  Top  panel:  Mean  surface  ocean  A14C  (Hughen  et  al.,  2004a).  The  mean  A14C 
during  the  YD,  BA,  and  LGM  are  drawn  with  dashed  lines,  and  annotated.  Bottom  panels: 
Compilation  of  pairs  of  benthic  and  planktonic  foraminifer  14 C  ages  from  the  North  Pacific. 
Plotted  here  is  A14Co  —  A14C5UM  the  difference  between  the  estimated  A14C  of  water  in 
which  the  benthic  foraminifera  grew  and  the  mean  surface  ocean  A14C  at  that  time.  The 
solid  grey  line  is  the  mean  North  Pacific  A*4C  profile  from  the  WOCE  dataset  minus  the 
A14C  of  a  mean  PIPN  surface  ocean  where  R  =  400 y  (A14CSU?  =  — 49.6%o)  (Key  et  al., 
2002).  The  dashed  grey  line  is  the  same  profile,  but  where  the  14 C  age  everywhere  in  the 
ocean  was  decreased  by  exactly  the  apparent  age  of  the  mean  surface  water  for  that  time 
period.  If  the  ocean  circulation  was  exactly  the  same,  but  the  atmospheric  and  surface 
ocean  l4C  was  higher,  the  subsurface  A14C  would  be  the  dashed  grey  profile. 
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of  14C.  Given  two  samples  that  are  1000  y  apart  in  age,  the  A14C  difference  between  the 
two  is  greater  the  younger  the  samples  are,  and  the  A14C  difference  is  smaller  the  older  the 
samples  are.  The  dashed  grey  curve  is  the  modern  A14C  profile,  but  where  for  each  time 
period,  the  apparent  age  of  the  mean  surface  water  A14C  is  subtracted. 

A .3  Apparent  ventilation  age 

There  is  a  long-term  decrease  in  14C  production  in  the  atmosphere  through  the  late  Pleis¬ 
tocene,  which  is  thought  to  be  the  consequence  of  a  strengthening  geomagnetic  field  and 
increase  in  shielding  of  Earth  from  cosmic  rays  (Bard,  1998)  (see  Figure  A-5).  This  has 
the  effect  that  a  sample  generally  appears  relatively  younger  them  its  calendar  age,  since 
atmospheric  14C  was  higher  in  the  past  than  today.  The  decrease  in  atmospheric  14C  was 
not  monotonic,  there  are  superimposed  shorter-scale  fluctuations  in  14 C  production  owing 
to  variations  in  solar  activity  and  the  geomagnetic  field.  In  addition,  changes  in  the  global 
carbon  cycle  and  ocean  circulation,  both  important  reservoirs  of  14C,  affect  the  14C/12C 
ratio  in  the  atmosphere. 

Since  a  benthic  foraminifer  grows  in  water  composed  of  constituents  which  were  last  in 
contact  with  the  atmosphere  at  some  point  in  the  past,  the  difference  between  a  benthic  and 
planktonic  14C  age  is  likely  to  be  somewhat  smaller  than  it  would  have  been  if  atmospheric 
14C  was  constant.  Adkins  and  Boyle  (1997)  proposed  a  method  of  interpreting  the  14C 
ages  of  benthic  and  planktonic  foraminifera  that  takes  into  account  changes  in  atmospheric 
A14C. 

The  Adkins  and  Boyle  (1997)  method  uses  the  simplifying  assumption  that  the  parcel 
of  water  in  which  the  benthic  foraminifer  grew  descended  from  the  sea  surface  as  a  closed 
system.  The  authors  acknowlege  that  this  assumption  is  never  true  for  a  parcel  of  water  in 
the  real  ocean,  but  that  calculating  the  14C  content  of  this  theoretical  water  mass  allows 
them  to  take  into  account  changing  atmospheric  A14C.  The  trajectory  of  the  A14C  of  the 
water  parcel,  once  it  is  removed  from  contact  with  the  atmosphere,  is 

A14C(f)  =  1000  (A1-^q0)  +  l)e~At  ~  1  , 

where  A14C(to)  is  the  A14C  of  the  parcel  at  the  moment  it  leaves  the  surface  ocean,  and  t 
is  time. 
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Figure  A-5:  Schematic  of  calculation  of  tv  and  tVAB-  In  the  upper  panel  in  grey  is  atmo¬ 
spheric  A14C  (Reimer  et  ah,  2004),  in  blue  is  the  mean  surface  ocean  Al4C  (Hughen  et  al., 
2004b).  In  red  is  the  trajectory  of  14C  decay  backward  in  time  of  a  planktonic  sample,  in 
green  is  the  decay  trajectory  for  a  benthic  sample.  For  each  of  these  4  curves,  the  la  and 
2 a  uncertainties  are  contoured.  In  the  inset  panel,  the  orange  curve  is  the  trajectory  of  a 
hypothetical  water  parcel  in  which  the  benthic  foraminifera  grew.  The  probability  distri¬ 
butions  (PDFs)  of  the  intersection  of  the  decay  trajectories  with  the  A14C  time  series  are 
drawn,  where  the  median  of  the  PDFs  are  indicated  with  dashed  lines.  From  left  to  right, 
the  PDFs  are  the  intersection  of  the  planktonic  decay  trajectory  with  the  surface  ocean 
curve  (taking  into  account  AT?,  the  local  anomaly  from  the  global  mean  age  of  surface  wa¬ 
ter),  the  intersection  of  the  water  parcel  decay  trajectory  with  the  mean  surface  ocean,  and 
the  intersection  of  the  water  parcel  decay  trajectory  with  the  atmosphere.  The  apparent 
ventilation  age  calculated  as  in  Adkins  and  Boyle  (1997)  is  tVAB •  In  this  study.  I  use  the 
quantity  tu. 
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Figure  A-6:  Compilation  of  pairs  of  benthic  and  planktonic  foraminifer  14C  age  from  the 
North  Pacific.  Plotted  here  is  tv,  the  difference  between  the  14C  age  in  ky  (using  the  true 
half  life)  of  benthics  and  planktonics  taking  into  account  change  in  atmospheric  A14C, 
assuming  that  the  water  in  which  the  benthic  foraminifera  grew  descended  from  the  surface 
ocean  as  a  closed  system.  The  grey  shaded  region,  the  colors,  symbols,  and  acronyms  are 
the  same  as  in  Figure  A-2. 


Assuming  that  the  A14C  of  benthic  foraminifera  is  the  same  as  A14C  of  the  water  in 
which  they  live,  and  that  the  parcel  is  a  closed  system,  the  trajectory  of  the  parcel  of 
water  is  the  same  as  the  “undecay”  trajectory  of  the  benthic  sample  (Equation  A. 3).  This 
is  illustrated  as  the  green  trajectory  in  Figure  A-5,  where  the  orange  portion  beyond  the 
calendar  age  of  the  sample  is  the  trajectory  of  the  hypothetical  water  parcel.  The  la  and 
2a  uncertainty  from  measurement  error  are  shown  as  shaded  contours. 

The  apparent  ventilation  age  of  the  parcel  of  water  with  respect  to  the  atmosphere  is 
the  difference  between  the  calendar  age  of  the  sample  and  the  intersection  of  this  trajectory 
with  the  Intcal2004  record  of  atmospheric  Al4C  (Reimer  et  ah,  2004).  Then,  the  estimate 
of  the  local  apparent  age  of  surface  water  (R  +  A R)  is  subtracted  to  find  the  apparent  age 
with  respect  to  the  local  surface  ocean,  tVAB'  This  calculation  is  functionally  the  same  as 
finding  the  difference  between  the  calibrated  age  of  the  uncorrected  benthic  sample  with 
respect  to  the  atmosphere  and  the  calibrated  age  of  the  planktonic  sample  (corrected  for 
R  +  A R)  with  respect  to  the  surface  ocean,  then  subtracting  R  T  A R.  Since  tv  takes  into 
account  A /?,  like  A^_p,  it  is  the  apparent  age  of  subsurface  water  relative  to  local  surface 
water,  rather  than  from  a  mean  surface  ocean. 

This  hypothetical  parcel  of  water  is  not  one  that  I  expect  to  have  ever  existed  in  the  real 
ocean,  where  it  is  highly  unlikely  that  a  water  parcel  would  descend  from  the  sea  surface 
without  being  transformed  by  mixing  with  surrounding  water  masses.  If  we  consider  the 
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slightly  more  realistic  scenario  in  which  the  parcel  is  a  mixture  of  several  end  members. 


A14C  =  aA14Ci  +6A14C2  +  cA14C3 
1  =  a  +  b  +  c 

The  A14C  of  the  water  parcel  is  the  sum  of  the  A14C  of  each  constituent  water  mass,  each 
weighted  by  the  fraction  of  the  parcel  it  composes.  In  this  case,  change  in  the  A14C  of 
the  parcel  could  be  the  result  of  a  change  in  the  fraction  of  a  constituent,  or  a  change  in 
the  A14C  of  a  constituent.  So  though  the  water  parcel  itself  is  fictional,  it  is  useful  in 
that  it  includes  a  correction  which  takes  into  account  a  known  phenomenon,  the  variable 
atmospheric  A14C. 

I  modified  the  Adkins  and  Boyle  (1997)  method  slightly  by  finding  the  intersection  of 
the  parcel  of  water  with  a  mean  surface  ocean  curve  (Hughen  et  ah,  2004b),  so  that  tv  is 
relative  to  the  mean  surface  ocean  at  the  time  the  hypothetical  water  parcel  would  have 
been  at  the  sea  surface,  rather  than  relative  to  local  surface  water.  This  is  functionally  the 
same  as  finding  the  difference  between  the  calendar  age  of  the  uncorrected  benthic  sample 
with  respect  to  the  mean  surface  ocean,  and  the  calendar  age  of  the  planktonic  sample 
with  respect  to  the  surface  ocean  taking  into  account  A R.  I  report  the  uncertainty  of  this 
difference  as  the  la  and  2a  ranges  of  the  PDF  of  the  difference  between  the  PDFs  of  the 
calibrated  calendar  ages.  These  PDFs  are  typically  non-Gaussian  and  can  be  obtained  from 
the  downloadable  versions  of  Calib  (http://www.calib.org).  I  plot  time  slices  of  tv  in 
Figure  A-6  with  the  modern  profile  of  apparent  age,  where  R  is  subtracted.  I  do  not  take 
into  account  the  local  A R  when  finding  the  intersection  of  the  benthic  trajectory  with  the 
surface  ocean  because  this  would  imply  that  the  local  A R  is  that  same  as  that  in  the  region 
where  subsurface  water  is  formed.  By  making  each  tv  relative  to  a  mean  surface  ocean,  I 
make  no  assumptions  about  the  source  region  of  the  subsurface  water. 


A. 4  Discussion 

The  greatest  advantage  of  A 5_p  is  its  simplicity — it  is  the  difference  between  two  14C  mea- 
surements.  However,  it  does  not  take  into  account  either  variability  in  the  apparent  age  of 
the  surface-ocean,  or  the  known  decrease  of  atmospheric  A14C  through  the  Pleistocene.  In 
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Figure  A-7:  from  a  compilation  of  pairs  of  benthic  and  planktonic  foraminifer  14 C  age 

from  the  North  Pacific,  la  uncertainties  are  plotted  at  each  data  point.  The  three  time 
slices  in  Figure  A-2  are  shaded  in  blue,  pink,  and  green,  and  are  the  YD,  BA,  and  LGM. 
respectively.  The  colors  for  each  time  series  is  the  same  as  in  Figure  A-2. 
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a  region  with  positive  A/?,  A^_p  appears  smaller  than  in  a  region  where  A R  —  0,  even  if 
the  14  C  content  of  the  subsurface  water  is  the  same  at  both  sites.  In  Figure  A-2,  I  expect 
that  there  is  some  spread  of  data  points  to  the  left  (to  smaller  values  of  A 5_p)  from  the 
effect  of  a  positive  AR  in  some  regions.  Since  A R  is  very  poorly  constrained  in  the  past, 
it  is  conservative  to  compare  A^_p  to  a  modern  profile  which  shows  a  likely  range  as  a 
function  of  A R.  However,  since  this  range  is  on  the  same  order  as  the  magnitude  of  change 
in  A during  the  time  periods  of  interest,  the  three  time  slices  do  not  look  significantly 
different  from  one  another  (see  the  individual  time  series  in  Figure  A-7).  Just  regarding  the 
data  points  that  fall  within  the  time  slices  with  95%  certainty,  during  the  YD  two  out  of 
three  points  fall  below  the  modern  range,  during  the  BA  one  out  of  nine  points  falls  below 
the  modern  range.  It  would  be  difficult  to  argue  that  the  data  for  these  time  intervals  are 
significantly  different  than  the  modern  profile.  During  the  LGM,  one  point  falls  below  and 
three  fall  above  the  modern  range  out  of  14  points,  where  there  appears  to  be  the  same  or 
slightly  greater  A &_p  at  2-3  km  depth. 

For  most  sites,  the  time  series  of  A 6_p  is  variable  within  the  time  period  of  interest. 
Though  some  of  the  individual  records  may  contain  statistically  significant  changes  during 
the  deglaciation,  when  amassed  (Figure  A-7),  it  is  evident  that  all  of  the  records  together 
do  not  show  the  same  story.  And  in  fact,  some  of  the  relative  maxima  and  minima  fail 
within  the  same  time  period.  The  uncertainty  in  the  calibration  of  planktonic  i4C  ages  is 
great  enough  that  it  would  not  be  possible  to  tell  whether  changes  occurring  in  each  core 
are  synchronous.  It  is  also  possible  that  there  were  indeed  significant  changes  in  A14C  that 
were  too  brief  to  resolve  in  sediment  cores. 

Calculation  of  A14Co  depends  on  one  additional  piece  of  information,  A R.  A R  is  used  to 
estimate  the  calendar  age  of  the  sample  from  the  planktonic  A14C  measurement.  Knowing 
the  calendar  age  of  the  sample  allows  one  to  transform  the  difference  in  14C  content  of  the 
benthic  and  planktonic  samples  (or  the  mean  surface  ocean)  from  14  C  years  to  A14C,  and 
estimate  the  profile  of  14C  activity  in  the  subsurface  water.  The  profiles  of  A14Co— A14Csti, 
are  different  from  the  profiles  of  A&_p  because  A14C  changes  exponentially  with  time. 

During  the  YD,  A14Co— A14Csiy  appears  significantly  higer  than  the  expected  profile  if 
the  circulation  was  the  same  as  today  (dashed  grey  line  in  Figure  A-4a).  During  the  BA  and 
LGM,  the  data  points  fall  near  the  expected  profile  and  at  values  slightly  higher.  During 
the  LGM,  the  exception  is  from  the  deepest  core  in  the  profile  (ES),  which  is  significantly 
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Figure  A-8:  A14Co— A14CSU,  from  a  compilation  of  pairs  of  benthic  and  planktonic 
foraminifer  14C  age  from  the  North  Pacific,  la  uncertainties  are  plotted  at  each  data,  point. 
The  three  time  slices  in  Figure  A-4  are  shaded  in  blue,  pink,  and  green,  and  are  the  YD, 
BA,  and  LGM,  respectively.  The  colors  for  each  time  series  is  the  same  as  in  Figure  A-4. 
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lower  by  150  %o.  From  these  data,  it  appears  that  during  the  YD  and  BA,  subsurface  Pacific 
water  above  2700  m  water  depth  may  have  been  slightly  enriched  in  14C  relative  to  today, 
and  during  the  LGM,  they  were  the  same  or  slightly  enriched  relative  to  today.  Deep  Pacific 
water  deeper  than  2700  m,  from  one  data  point  from  the  Emperor  Seamounts,  may  have 
been  significantly  depleted  in  14C. 

Compared  to  A^_p,  the  amplitude  of  changes  of  tv  through  time  are  larger  (compare 
Figure  A-9  and  Figure  A-7).  Relative  to  the  modern  profile,  tv  may  be  lower  during  the  YD, 
though  again,  this  is  hard  to  tell,  given  the  small  number  and  large  scatter  of  data  points. 
The  majority  of  data  points  from  the  BA  and  LGM  time  slices  are  higher  than  the  modern 
profile,  which  suggests  that,  subsurface  North  Pacific  water  was  out  of  contact  with  the 
atmosphere  for  longer  during  these  time  periods,  compared  to  today.  This  interpretation 
of  the  tv  data  is  contrary  to  that  of  A14Cq— A14Csli,  data. 


Ill 


YD  BA 
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Figure  A-9:  tv  from  a  compilation  of  pairs  of  benthic  and  planktonic  foraminifer  14C  age 
from  the  North  Pacific.  1<t  uncertainties  are  plotted  at  each  data  point.  The  three  time 
slices  in  Figure  A-4  are  shaded  in  blue,  pink,  and  green,  and  are  the  YD,  BA,  and  LGM, 
respectively.  The  colors  for  each  time  series  is  the  same  as  in  Figure  A-6. 
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Table  A.l:  Compilation  of  pairs  of  benthic  and  planktonic  14C  dates.  I  report  the  maxima 
and  minima  of  the  la  and  2 a  ranges  of  each  quantity.  When  the  calendar  age  is  during  the 
YD,  BA  and  LGM  (as  defined  in  Stuiver  et  al.  (1998))  with  95%  (2a)  certainty,  the  time 
period  is  bold.  I  use  the  same  A R  as  in  the  original  publications. 
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